Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,165,051 A→C intergenic (‑97/‑265) coaA ← / → BW25113_RS20630 type I pantothenate kinase/tRNA‑Thr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,165,0510AC93.1% 98.7 / ‑5.7 29intergenic (‑97/‑265)coaA/BW25113_RS20630type I pantothenate kinase/tRNA‑Thr
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (16/11);  total (17/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCGTCATTAACGTTTGCTCTTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAA  >  NZ_CP009273/4164922‑4165184
                                                                                                                                 |                                                                                                                                     
aGGCGTCATTAACGTTTGCTCTTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATggg                                                                                                                              <  2:23431/139‑1 (MQ=255)
                tGCTCTTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGaa                                                                                                               >  1:392988/1‑138 (MQ=255)
                                 aaGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAg                                                                                             >  1:278048/1‑139 (MQ=255)
                                              gCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGaa                                                                                  <  2:344191/137‑1 (MQ=255)
                                                           tATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACa                                                                     <  1:40721/137‑1 (MQ=255)
                                                            aTCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAg                                                                                             <  2:115869/112‑1 (MQ=255)
                                                            aTCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAg                                                                                             >  1:115869/1‑112 (MQ=255)
                                                            aTCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGa                                                                   >  1:226124/1‑138 (MQ=255)
                                                                   aaGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCAt                                                           <  1:376659/139‑1 (MQ=255)
                                                                     gCGATATTTGCCATGGCGCGACACAGTTACGTTGAATGGGTGAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGAAATGTGACAGAAATTCCATCt                                                         >  2:513628/1‑139 (MQ=255)
                                                                           tttGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATaaa                                                   <  2:427529/139‑1 (MQ=255)
                                                                                   ggCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCa                                                            <  2:91953/122‑1 (MQ=255)
                                                                                   ggCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCa                                                            >  1:91953/1‑122 (MQ=255)
                                                                                               ttCTGTTGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAg                                >  1:523330/1‑138 (MQ=255)
                                                                                                    ttGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAgcg                           >  1:226496/1‑138 (MQ=255)
                                                                                                    ttGAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAgcg                           >  1:85128/1‑138 (MQ=255)
                                                                                                    ttGAATGGGTTAATGGGCAGGAGGGTAACCACAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAgcg                           >  2:497886/1‑138 (MQ=255)
                                                                                                      gAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCa                         <  2:523330/138‑1 (MQ=255)
                                                                                                      gAATGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCa                         <  1:175285/138‑1 (MQ=255)
                                                                                                      gAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCa                         <  2:409349/138‑1 (MQ=255)
                                                                                                         tGGGTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAttgt                     >  2:341448/1‑139 (MQ=255)
                                                                                                         tGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGATTAGCAGGAGCGCAttgt                     >  2:163958/1‑139 (MQ=255)
                                                                                                            gTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTAt                                          >  2:29592/1‑115 (MQ=255)
                                                                                                            gTTAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTAt                                          <  1:29592/115‑1 (MQ=255)
                                                                                                              tAATGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAttgtttat                 >  2:150152/1‑135 (MQ=255)
                                                                                                                aTGGGCAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGca               <  2:85128/138‑1 (MQ=255)
                                                                                                                     cAGGAGGGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAatga         >  1:243549/1‑139 (MQ=255)
                                                                                                                      aggaggGTAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAatgat        >  2:393456/1‑139 (MQ=255)
                                                                                                                             tAACCCCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGaa  >  2:457075/1‑138 (MQ=255)
                                                                                                                                 |                                                                                                                                     
AGGCGTCATTAACGTTTGCTCTTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAA  >  NZ_CP009273/4164922‑4165184

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: