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.97.6% 181.6 / ‑1.8 42intergenic (‑120/‑242)coaA/thrUpantothenate kinase/tRNA‑Thr
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (20/21);  minor base C (1/0);  total (21/21)
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

CTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCT  >  NC_000913/4175022‑4175271
                                                                                                                            |                                                                                                                             
cTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCa                                                                                                                 <  6:152978/138‑1 (MQ=255)
cTCATAAGAATGATTCTGGCCGTCATAGCTCTATCGCCAACGCGATATTTGCAATGCCGCGTCACAGTACTGTTGAAAGGGTTAAAGGGCAGGAGGGTAACACCAGATGGGGAGGAATAAAAACAAAAATCTCAGCCAGGAAGCTATCg                                                                                                      >  6:87955/1‑149 (MQ=255)
   aTAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGtt                                                                                                    >  4:61642/1‑148 (MQ=255)
       gCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGaa                                                                                                 >  2:130208/1‑147 (MQ=255)
            tttCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCaa                                                                                           <  5:26574/148‑1 (MQ=255)
              tCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAAtg                                                                                         <  3:61642/148‑1 (MQ=255)
                       cATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGatat                                                                               >  1:75124/1‑149 (MQ=255)
                               tATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCt                                                                        >  6:4611/1‑148 (MQ=255)
                                 tCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCt                                                                        <  2:75124/146‑1 (MQ=255)
                                      aaaGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATaaaa                                                                 >  6:48017/1‑148 (MQ=255)
                                         gCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGTGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGaa                                                             <  1:130208/149‑1 (MQ=255)
                                          cGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAAt                                                            <  5:4611/149‑1 (MQ=255)
                                                  gCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACa                                                                                    <  3:52343/117‑1 (MQ=255)
                                                  gCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACa                                                                                    >  4:52343/1‑117 (MQ=255)
                                                    cATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAg                                                  >  1:12466/1‑149 (MQ=255)
                                                             tCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAgcg                                          >  1:81080/1‑148 (MQ=255)
                                                               acaGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAt                                       <  2:12466/149‑1 (MQ=255)
                                                               acaGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAt                                       >  5:135558/1‑149 (MQ=255)
                                                                 aGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAttg                                     <  2:81080/149‑1 (MQ=255)
                                                                       gTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGc                                                                                                                   <  2:173357/65‑1 (MQ=255)
                                                                       gTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGc                                                                                                                   >  1:173357/1‑65 (MQ=255)
                                                                        ttGTATGGGTTAATGGGCTGGTGGGTAACACCAGATGGGGAGGATTAATAAGAAAAATCTCTGCCTGGAAGGTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAATAGAATTATGGTTTAGCAGGAGCGCATTGTTGa                                 <  5:87955/146‑1 (MQ=255)
                                                                        ttGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGca                              >  6:16578/1‑149 (MQ=255)
                                                                          gAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGa                                 >  3:55187/1‑144 (MQ=255)
                                                                          gAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGa                                 <  4:55187/144‑1 (MQ=255)
                                                                          gAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGcac                             <  5:16578/148‑1 (MQ=255)
                                                                           aaTGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGaaa                                                                                                >  2:138588/1‑80 (MQ=255)
                                                                           aaTGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGaaa                                                                                                <  1:138588/80‑1 (MQ=255)
                                                                            aTGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAgc                                           <  1:85119/132‑1 (MQ=255)
                                                                            aTGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAgc                                           >  2:85119/1‑132 (MQ=255)
                                                                                gTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAatgatg                      >  6:100787/1‑149 (MQ=255)
                                                                                 ttAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGt                     <  5:48017/149‑1 (MQ=255)
                                                                                          aggaggGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGa                                                                                      >  1:163322/1‑75 (MQ=255)
                                                                                          aggaggGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGa                                                                                      <  2:163322/75‑1 (MQ=255)
                                                                                          aggaggGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCAt                                                                          <  1:166322/87‑1 (MQ=255)
                                                                                          aggaggGTAACACCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCAt                                                                          >  2:166322/1‑87 (MQ=255)
                                                                                                    cacCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCt  <  5:100787/149‑1 (MQ=255)
                                                                                                     acCAGATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAg                                                  >  2:94161/1‑100 (MQ=255)
                                                                                                         gATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATg                                                        >  6:164647/1‑90 (MQ=255)
                                                                                                         gATGGGGAGGAATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATg                                                        <  5:164647/90‑1 (MQ=255)
                                                                                                               gaggaATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAttg                                     <  7:96488/103‑1 (MQ=255)
                                                                                                               gaggaATAATAAGAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCAttg                                     >  8:96488/1‑103 (MQ=255)
                                                                                                                            |                                                                                                                             
CTCATAAGCATGTTTCTGGCGGTCATAGCTCTATCGGCAAAGCGATATTTGCCATGGCGCGTCACAGTTCTGTTGAATGGGTTAATGGGCAGGAGGGTAACACCAGATGGGGAGGAATAATAAGAAAAAATCTCTGCCAGGAAGCTATCGTTGAAAAGCAATGTGACAGATATTCCATCTGATAAAAAGAATTATGGTTTAGCAGGAGCGCATTGTTGAGCACAATGATGTTGAAAAAGTGTGCTAATCT  >  NC_000913/4175022‑4175271

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: