Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,143 C→T 16.7% intergenic (+541/+73) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,270,1430CT16.7% 29.7 / 2.6 26intergenic (+541/+73)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base C (14/6);  new base T (2/2);  total (18/8)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

GCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTC  >  NC_000913/3270012‑3270289
                                                                                                                                   |                                                                                                                                                  
gCTTTCTTTTGTTTTTTGTGTCGTGGAATGCGTTTTATATTTCTGCCACTCCAAAAGCATAAACCCGCCGTAGCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattca                                                                                                                                  <  1:976248/150‑1 (MQ=11)
          ttttttaGTGTCGTGGACTGTCTTTATTCTAACGGGCACTCAAAAATCAAAAACCCGCCGCATCGGGTTTTTACGTAAACCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGcc                                                                                                                        <  1:1339091/145‑1 (MQ=11)
                               cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGCCAGGCATTCATCGTTAAGCGCCTACAAAAACAAAAACccgc                                                                                                   >  1:334846/1‑150 (MQ=11)
                               cattcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGAAAGTCCTTAATCGTTTAGCCCCGCCAAAACCCAAAAAccgc                                                                                                   >  1:1338272/1‑150 (MQ=11)
                                attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGGTTTCTTTTGGGTATAGCGTCGTGTAAAGTCATTCATATTTAAACACCAACAAAAACAAAAACccgga                                                                                                  >  2:744199/1‑148 (MQ=11)
                                   catCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTAGGGTAAAGAGTAGTGTAAAGTCAATCAAAATTGAACACGAACAAAAGCAAAAACccgcagca                                                                                               >  2:1243835/1‑146 (MQ=11)
                                        ttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGCCCGATAAGCCGCGTTCTTTTGGTTATAGCGTCGTGGCCCGGCCTTCAACGTAGTCCGCCACCAACAGCCAACACCCGCCGaccgcgg                                                                                          >  2:851757/1‑144 (MQ=11)
                                           tGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTGGGGTATAGCGTCGTGGAGAGTCCTTCAAGATTGTGCGCGAGCAAACGCCAAAACCCGCAGAAGggggggt                                                                                       >  1:835093/1‑147 (MQ=11)
                                                         cagaaaCCCGCCCCATCGGGTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAgg                                                                         <  1:724815/147‑1 (MQ=11)
                                                         cAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATATTTCAGCACCGCCAAAAGCAAAAACCCGCCGCAGCGCGGTTTTACGTAAATGAgg                                                                         >  2:385400/1‑150 (MQ=11)
                                                                 cgccgAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGGCAATAATCTTTCTGCCCCTCCAAAAACAAAAACCCGCCGAAGCGGGTTTTTAAGTTAATCATGTGAAACTg                                                                 >  1:998657/1‑150 (MQ=11)
                                                                     gAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGTCAATCATCTTTAAGCCCCTACAAAAACAAAAACCCGCAGAAGAGGGGTTTTAAGTTAATCAGGTGCAACttaacc                                                             >  1:701178/1‑145 (MQ=11)
                                                                                 tACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTAGGGTATAGCGTCGTGGCCAGTCAATCATAATTATCCACCACCAAAAGCAAAAACCCGCCGAAACGCGGTTTTACGTTAATAAGGTGAAACTTCCCCATAAGCCGGGtt                                                 >  2:295461/1‑150 (MQ=14)
                                                                                   cGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTAAACTTTCAGCCCCTCCAAAAGCAAAAACCCGCCGAAGCCGGGTTTTAAGTAAATCAGGGGCAACTGACAGATAAACCGGGTTCt                                               >  2:984514/1‑150 (MQ=14)
                                                                                                             aGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGCTAAGCCGGGTTCTGTCGTGGTCAGTCATTCATCTAGGcc                     >  2:803485/1‑150 (MQ=255)
                                                                                                              gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt                                             >  1:732292/1‑125 (MQ=21)
                                                                                                              gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt                                             <  2:732292/125‑1 (MQ=21)
                                                                                                               ccGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTAATATTTATGCCCCTCCAAAAGCAAAAACCAGCCGCAGCGGGTTTTTACGTAAATAAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAg                   >  1:1130885/1‑150 (MQ=18)
                                                                                                                cGCTTTCTTTTGGGTCTAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGc                  <  1:452792/150‑1 (MQ=255)
                                                                                                                cGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGGTTTTACGGAAATCAGGTGAAACTGACCGAAAAGCCGGGTTCTGTAGTGGACAGTCATTCATCTAGGCCAGc                  >  1:1291514/1‑150 (MQ=255)
                                                                                                                 gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCggg                                                   <  2:572936/116‑1 (MQ=14)
                                                                                                                 gCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCggg                                                   >  1:572936/1‑116 (MQ=9)
                                                                                                                         ttGGGTCTAGCGTCGTGTTAAGTGATTCATCCTTCTGACCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCa         <  1:249640/150‑1 (MQ=21)
                                                                                                                                taGCGTCGTGGACAGTCATTCATCTTTCTGACCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGGGGACAGTGATTCATCTTGGCCAGCAATCGCTCACTGGCTc  >  1:76377/1‑150 (MQ=255)
                                                                                                                                  gCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCg                                                     <  1:73376/97‑1 (MQ=35)
                                                                                                                                  gCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCg                                                     >  2:73376/1‑97 (MQ=35)
                                                                                                                                   |                                                                                                                                                  
GCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTC  >  NC_000913/3270012‑3270289

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