Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,913,204 G→A 10.8% intergenic (‑186/‑96) PP_4317 ← / → PP_4318 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,913,2040GA10.8% 95.0 / 3.7 37intergenic (‑186/‑96)PP_4317/PP_4318hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base G (18/15);  new base A (2/2);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

GAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTG  >  NC_002947/4913079‑4913312
                                                                                                                             |                                                                                                            
gAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTat                                                                                                    >  1:64606/1‑135 (MQ=21)
gAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAg                                                                                                    >  1:147554/1‑136 (MQ=35)
aaGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTat                                                                                                    >  1:15779/2‑135 (MQ=18)
             ggggACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAAt                                                                                       >  1:277099/1‑136 (MQ=40)
             ggggACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAAt                                                                                       >  1:514752/1‑136 (MQ=40)
                gACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCa                                                                                    >  1:222673/1‑136 (MQ=38)
                     cGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAg                                                                               >  2:506176/1‑136 (MQ=37)
                     cGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCAAGTCAg                                                                               >  2:55137/1‑136 (MQ=37)
                       gAGCGCAGCGCTGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGa                                                                             <  2:358098/136‑1 (MQ=37)
                         gcgcagcgcaGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAg                                                                           >  1:457104/1‑136 (MQ=37)
                           gcagcgcagGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGAc                                                                         >  1:379505/1‑136 (MQ=37)
                             agcgcagGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGAccc                                                                       >  1:374651/1‑136 (MQ=34)
                                gcagGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAg                                                                    >  2:491821/1‑136 (MQ=34)
                                   ggggCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGa                                                                 <  2:346570/136‑1 (MQ=34)
                                     ggCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGAtg                                                               >  2:361681/1‑136 (MQ=34)
                                       ccGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGaa                                                                                         >  1:536658/1‑108 (MQ=35)
                                       ccGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGaa                                                                                         <  2:536658/108‑1 (MQ=35)
                                       ccGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGAtgtg                                                             <  1:498962/136‑1 (MQ=34)
                                       ccGGATGATGGGAGCGCAGCGTTTTTTGGTTACTCTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGAtgtg                                                             >  1:383134/1‑136 (MQ=21)
                                         ggATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTggg                                                           <  2:141389/136‑1 (MQ=34)
                                          gatgatGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGATCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTggg                                                           <  1:442724/135‑1 (MQ=21)
                                          gatgatGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGATCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTggg                                                           <  1:440208/135‑1 (MQ=21)
                                          gatgatGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGATCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTggg                                                           >  2:442724/1‑135 (MQ=21)
                                          gatgatGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGATCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTggg                                                           >  2:440208/1‑135 (MQ=21)
                                                         gCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGcgcg                                           >  1:40795/1‑136 (MQ=32)
                                                         gCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGcgcg                                           >  1:386975/1‑136 (MQ=32)
                                                          cGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGa                                          <  2:188878/136‑1 (MQ=32)
                                                          cGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGa                                          >  1:401757/1‑136 (MQ=32)
                                                             tttttGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTc                                       <  1:361681/136‑1 (MQ=32)
                                                                      aCTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGaaaaaa                              <  2:200203/136‑1 (MQ=32)
                                                                      aCTTTTTGTCGCGGTTGACAAAAAGTGACTGGCCTTAAGGGCCATACCTTTCTTTGTCTTCTTTAGATTTCGCTGAACTCCAGTCAGGTTGACCGGAGCTATGTGGGCTGATGTACTAGCGCGATTCGATaaaaaa                              <  2:249109/136‑1 (MQ=11)
                                                                                            aaGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCa        <  1:140980/136‑1 (MQ=32)
                                                                                             tgTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCaa       <  2:194838/135‑1 (MQ=17)
                                                                                             aGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCaa       <  2:109063/136‑1 (MQ=32)
                                                                                                gACTCGCCGTAATGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACt    <  2:207269/136‑1 (MQ=17)
                                                                                                  cTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTg  <  2:350545/136‑1 (MQ=31)
                                                                                                  cTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTg  <  2:420440/136‑1 (MQ=31)
                                                                                                                             |                                                                                                            
GAGGCCGGAGTGAGGGGACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCCATCCCTTTCAAGGTCTTCGGTAGAATTCGCTGAAATCCAGTCAGGAAGACCCGAGCGATGTGGGCTGATGTACTAGCGCGATTCGAGAAAAAAGCACCCGCCAGCGTAATGGCCAAACTTG  >  NC_002947/4913079‑4913312

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