Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,175,151 (A)6→5 intergenic (‑125/‑237) coaA ← / → thrU pantothenate kinase/tRNA‑Thr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,175,1460A.94.3% 129.2 / 6.4 36intergenic (‑120/‑242)coaA/thrUpantothenate kinase/tRNA‑Thr
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (15/18);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

ATTAACGTTTGCTCTTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTG  >  NC_000913/4175001‑4175287
                                                                                                                                                 |                                                                                                                                             
aTTAACGGTTGCTCTTTTATACTCATAAGCATGTTTCTGGCTGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACATCAGATGGGGAGGAATAATAAGaaaa                                                                                                                                            <  6:2191/149‑1 (MQ=255)
                  ataCTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTa                                                                                                                         >  3:112059/1‑149 (MQ=255)
                   taCTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTAt                                                                                                                        >  5:138621/1‑149 (MQ=255)
                     cTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCg                                                                                                                      <  6:129883/149‑1 (MQ=255)
                      tCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTa                                                                                                                         <  8:55212/145‑1 (MQ=255)
                      tCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTa                                                                                                                         >  7:55212/1‑145 (MQ=255)
                      tCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGt                                                                                                                     >  3:36361/1‑149 (MQ=255)
                            gCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATctct                                                                                                                                     <  6:46368/127‑1 (MQ=255)
                            gCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATctct                                                                                                                                     >  5:46368/1‑127 (MQ=255)
                            gCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGaaaa                                                                                                               <  1:142844/149‑1 (MQ=255)
                              aTGTTTCTGGCCGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAg                                                                                                              >  6:52886/1‑148 (MQ=255)
                                   tCTGGCGGTCATAGATCTATCGGCAAAGCGATATTTGCCATGGCGCGTCGCAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGGAACACGTGAGGGGGAGGAAGAATATGAAAGAGCTGTGCAAGGAAGGTAGGGTTGGAAAggcaggt                                                                                                        >  8:105969/1‑143 (MQ=255)
                                     tGGCGGTCATAGCTCTATCGGCAAAGCGATAATTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGa                                                                                                      >  8:39971/1‑149 (MQ=255)
                                       tcgatCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGa                                                                                                      >  4:22088/5‑147 (MQ=255)
                                       gCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGAc                                                                                                     <  4:58041/148‑1 (MQ=255)
                                                    tATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATc                                                                                         <  3:22088/147‑1 (MQ=255)
                                                      tCGGCAAAGCGATAATTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCt                                                                                        <  7:39971/146‑1 (MQ=255)
                                                        ggCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCAATCTtata                                                                                    >  3:107507/1‑148 (MQ=255)
                                                                 atatTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGTGCTGGAGGGTAACAACAGATGGGGAGGAATAATAAGAAAAATCTATGACAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATATGATAAAATGAATTa                                                                          >  4:132492/1‑149 (MQ=255)
                                                                      tGCCATGGCGCGCCACAGGGCTGTTGAATGCCTTACTGGGCAGGACGGGAACACCAGATCGGCAGCAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTCCAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGtt                                                                     <  6:133193/149‑1 (MQ=255)
                                                                         cATGGCGCTTCACAGTTCTGTTGAATTGGTTAATGGGCAGGAGGTTAACACCAGATTGGGATTAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAATAATTATGGTTTa                                                                   <  4:107507/148‑1 (MQ=255)
                                                                          aTGGCGCTTCACAGTTCTGTTGAATGGGTTAATGGTCATGAGGGTAACAACAGATGGGGAGGAATAATAATAAAAATCTCTGAAAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATATGATAAAAAGAATGATGGTTTAg                                                                  >  6:114066/1‑148 (MQ=255)
                                                                          aTGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGc                                                                 <  4:7423/149‑1 (MQ=255)
                                                                               gcgTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCt                                                                                        >  3:123793/1‑121 (MQ=255)
                                                                               gcgTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCt                                                                                        <  4:123793/121‑1 (MQ=255)
                                                                               gcgTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAg                                                            >  4:37924/1‑149 (MQ=255)
                                                                                   cacaGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCa                                                        >  5:74255/1‑149 (MQ=255)
                                                                                     caGTGCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAtt                                                      <  5:63673/149‑1 (MQ=255)
                                                                                        ttCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAttgtt                                                   <  5:52886/149‑1 (MQ=255)
                                                                                             ttGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGa                                                 <  1:144260/146‑1 (MQ=255)
                                                                                             ttGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGa                                                 <  4:36361/146‑1 (MQ=255)
                                                                                                            tgCAGGATGGTAACACCATATTGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAATAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGaaa                                <  7:105969/147‑1 (MQ=255)
                                                                                                                ggaggGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAgtg                           <  6:74255/149‑1 (MQ=255)
                                                                                                                 gaggGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAgtgt                          <  4:112059/149‑1 (MQ=255)
                                                                                                                                         ataataAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAATTGTGCTAATCTGCCCTCCGTTCGGCTg  >  7:103628/1‑149 (MQ=255)
                                                                                                                                           aataaGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAgt                            >  4:86115/1‑122 (MQ=255)
                                                                                                                                           aataaGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTACCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAgt                            <  3:86115/122‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
ATTAACGTTTGCTCTTTTATACTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCTGCCCTCCGTTCGGCTG  >  NC_000913/4175001‑4175287

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