Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,221,684 C→A 19.6% intergenic (+99/‑111) rpoZ → / → rpoZ DNA‑directed RNA polymerase subunit omega/bifunctional phosphopantothenoylcysteine decarboxylase/phosphopantothenate‑‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,221,6840CA19.6% 50.3 / 7.2 22intergenic (+99/‑111)rpoZ/rpoZDNA‑directed RNA polymerase subunit omega/bifunctional phosphopantothenoylcysteine decarboxylase/phosphopantothenate‑‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (10/6);  new base A (2/2);  total (14/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

GCCCTGTTATTTCAAGTGATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATATGATGATGTAAAGCGTAAATGTCATAGTGATATTTAAGGGAGATTGAATATGAAGAA  >  CP000730/1221548‑1221808
                                                                                                                                        |                                                                                                                            
gCCCTGTTATTTCAAGTGATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAt                                                                                                                          <  2:137461/141‑1 (MQ=255)
gCCCTGTTATTTCAAGTGATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAAt                                                                                                                          <  2:227804/141‑1 (MQ=255)
                        atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATACAATGATTGATGTTATAATAATACTaaa                                                                                                  >  2:128442/1‑141 (MQ=255)
                        atGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATACAATGATTGATGTTATAATAATACTaaa                                                                                                  >  2:73106/1‑141 (MQ=255)
                                    aGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCtata                                                                                      >  2:67668/1‑141 (MQ=255)
                                    aGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCtata                                                                                      >  2:21135/1‑141 (MQ=255)
                                                                 aGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATatg                                                         <  2:194077/141‑1 (MQ=255)
                                                                   aCGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATatgat                                                       >  1:169242/1‑141 (MQ=255)
                                                                                aCGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTTGTTTATTTAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATATGATGATGTAAAGCGTa                                          >  2:127125/1‑141 (MQ=255)
                                                                                aCGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATATGATGATGTAAAGCGTa                                          >  2:157961/1‑141 (MQ=255)
                                                                                      aGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATATGATGATGTAAAGCGTAAATGTc                                    >  1:179834/1‑141 (MQ=255)
                                                                                       ggCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATATGATGATGTAAAGCGTAAATGTCa                                   >  1:262263/1‑141 (MQ=255)
                                                                                       ggCATATCCAAATCAAACGCACCTATGTAGGCGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATATGATGATGTAAAGCGTAAATGTCa                                   >  2:105766/1‑141 (MQ=255)
                                                                                              ccAAATCAAACGCACCTATGTAGGTGTGTTTTTTTTTTTTTTTTTATCATTGATGTTATAATAATACTATATTTGTATTTATAAAAAATTAGTAAGCCTTTGTGCGCATTTGATGCTGTAACGCGTAGATGACATAGCGat                            >  2:162031/1‑141 (MQ=255)
                                                                                                   tCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATATGATGATGTAAAGCGTAAATGTCATAGTGATATTTa                       <  2:126615/141‑1 (MQ=255)
                                                                                                      aaCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATATGATGATGTAAAGCGTAAATGTCATAGTGATATTTAAgg                    <  1:218186/141‑1 (MQ=255)
                                                                                                      aaCGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATATGATGATGTAAAGCGTAAATGTCATAGTGATATTTAAgg                    <  1:137525/141‑1 (MQ=255)
                                                                                                      aaCGCACCTATGTAGGTGTGTTTTTTTAGTTTATACAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATATGATGATGTAAAGCGTAAATGTCATAGTGATATTTAAgg                    <  1:73106/141‑1 (MQ=255)
                                                                                                      aaCGCACCTATGTAGGTGTGTTTTTTTAGTTTATACAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATATGATGATGTAAAGCGTAAATGTCATAGTGATATTTAAgg                    <  1:128442/141‑1 (MQ=255)
                                                                                                        cGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATATGATGATGTAAAGCGTAAATGTCATAGTGATATTTAAGGga                  >  1:99964/1‑141 (MQ=255)
                                                                                                                        tgtTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATATGATGATGTAAAGCGTAAATGTCATAGTGATATTTAAGGGAGATTGAATATgaagaa  >  1:199859/1‑141 (MQ=255)
                                                                                                                        tgtTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATATGATGATGTAAAGCGTAAATGTCATAGTGATATTTAAGGGAGATTGAATATgaagaa  >  1:12534/1‑141 (MQ=255)
                                                                                                                                        |                                                                                                                            
GCCCTGTTATTTCAAGTGATTATTATGGTAAAGAATAGTTTCGACCATTAAAAATATGTGAAGCTAGACGTACATCAAATACGTCAAGGCATATCCAAATCAAACGCACCTATGTAGGTGTGTTTTTTTAGTTTATCCAATGATTGATGTTATAATAATACTAAATTTGTATCTATAAAAAAGTAATGAGCATTTGTGCGCATATGATGATGTAAAGCGTAAATGTCATAGTGATATTTAAGGGAGATTGAATATGAAGAA  >  CP000730/1221548‑1221808

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: