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)
                  ataCTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTa                                                                                                                         >  3:112059/1‑149 (MQ=255)
                   taCTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTAt                                                                                                                        >  5:138621/1‑149 (MQ=255)
                     cTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCg                                                                                                                      <  6:129883/149‑1 (MQ=255)
                      tCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTa                                                                                                                         <  8:55212/145‑1 (MQ=255)
                      tCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTa                                                                                                                         >  7:55212/1‑145 (MQ=255)
                      tCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGt                                                                                                                     >  3:36361/1‑149 (MQ=255)
                            gCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATctct                                                                                                                                     <  6:46368/127‑1 (MQ=255)
                            gCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATctct                                                                                                                                     >  5:46368/1‑127 (MQ=255)
                            gCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGaaaa                                                                                                               <  1:142844/149‑1 (MQ=255)
                              aTGTTTCTGGCCGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAg                                                                                                              >  6:52886/1‑148 (MQ=255)
                                   tCTGGCGGTCATAGATCTATCGGCAAAGCGATATTTGCCATGGCGCGTCGCAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGGAACACGTGAGGGGGAGGAAGAATATG‑AAAGAGCTGTGCAAGGAAGGTAGGGTTGGAAAggcaggt                                                                                                        >  8:105969/1‑143 (MQ=255)
                                     tGGCGGTCATAGCTCTATCGGCAAAGCGATAATTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGa                                                                                                      >  8:39971/1‑149 (MQ=255)
                                       tcgatCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGa                                                                                                      >  4:22088/5‑147 (MQ=255)
                                       gCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGAc                                                                                                     <  4:58041/148‑1 (MQ=255)
                                                    tATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATc                                                                                         <  3:22088/147‑1 (MQ=255)
                                                      tCGGCAAAGCGATAATTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCt                                                                                        <  7:39971/146‑1 (MQ=255)
                                                        ggCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCAATCTtata                                                                                    >  3:107507/1‑148 (MQ=255)
                                                                 atatTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGTGCTGGAGGGTAACAACAGATGGGGAGGAATAATAAG‑AAAAATCTATGACAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATATGATAAAATGAATTa                                                                          >  4:132492/1‑149 (MQ=255)
                                                                      tGCCATGGCGCGCCACAGGGCTGTTGAATGCCTTACTGGGCAGGACGGGAACACCAGATCGGCAGCAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTCCAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGtt                                                                     <  6:133193/149‑1 (MQ=255)
                                                                         cATGGCGCTTCACAGTTCTGTTGAATTGGTTAATGGGCAGGAGGTTAACACCAGATTGGGATTAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAATAATTATGGTTTa                                                                   <  4:107507/148‑1 (MQ=255)
                                                                          aTGGCGCTTCACAGTTCTGTTGAATGGGTTAATGGTCATGAGGGTAACAACAGATGGGGAGGAATAATAA‑TAAAAATCTCTGAAAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATATGATAAAAAGAATGATGGTTTAg                                                                  >  6:114066/1‑148 (MQ=255)
                                                                          aTGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGc                                                                 <  4:7423/149‑1 (MQ=255)
                                                                               gcgTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCt                                                                                        >  3:123793/1‑121 (MQ=255)
                                                                               gcgTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCt                                                                                        <  4:123793/121‑1 (MQ=255)
                                                                               gcgTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAg                                                            >  4:37924/1‑149 (MQ=255)
                                                                                   cacaGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCa                                                        >  5:74255/1‑149 (MQ=255)
                                                                                     caGTGCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAtt                                                      <  5:63673/149‑1 (MQ=255)
                                                                                        ttCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAttgtt                                                   <  5:52886/149‑1 (MQ=255)
                                                                                             ttGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGa                                                 <  1:144260/146‑1 (MQ=255)
                                                                                             ttGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGa                                                 <  4:36361/146‑1 (MQ=255)
                                                                                                            tgCAGGATGGTAACACCATATTGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAATAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGaaa                                <  7:105969/147‑1 (MQ=255)
                                                                                                                ggaggGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAgtg                           <  6:74255/149‑1 (MQ=255)
                                                                                                                 gaggGTAACACCAGATGGGGAGGAATAATAAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAgtgt                          <  4:112059/149‑1 (MQ=255)
                                                                                                                                         ataataAG‑AAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAATTGTGCTAATCTGCCCTCCGTTCGGCTg  >  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: ‑