Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,030 T→C 13.8% intergenic (+428/+186) yhaC → / ← rnpB uncharacterized protein YhaC/RNase P catalytic RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,270,0300TC13.8% 59.6 / 4.3 34intergenic (+428/+186)yhaC/rnpBuncharacterized protein YhaC/RNase P catalytic RNA component
Reads supporting (aligned to +/- strand):  ref base T (6/21);  new base C (2/2);  total (11/23)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

AACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGAT  >  NC_000913/3269883‑3270119
                                                                                                                                                   |                                                                                         
aaCTGACCGATAAGCCGATTTCTTTTTCCTATAGTGTCGGGGAAAGTGATTCATCTTTCTGCCCCTCCAAAATTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAgtgt                                                                                         <  1:538847/150‑1 (MQ=11)
                    tctttTGGGTATTGTGTCGTGGACTGTCTTTCATCTTTCTGGCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTAACCGATAAGCCGCTTTCTTTTGGGTACAGTGTCGTGGACAGTCATTCATCtt                                                                     <  1:2461600/150‑1 (MQ=11)
                                     cGTGGACTGTCTTTACTCTTAATCCCACTCAAAAAGATAAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGc                                                    <  1:2110783/150‑1 (MQ=11)
                                                   catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTg                                                                                     >  2:770499/1‑103 (MQ=25)
                                                    taatttATGCCACTCAAAAATAACAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAg                                     <  1:990779/147‑1 (MQ=11)
                                                                 tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTgg                                                                                    >  1:715106/1‑90 (MQ=18)
                                                                 tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTgg                                                                                    <  2:715106/90‑1 (MQ=18)
                                                                 tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattcat                                                                        <  1:1786489/102‑1 (MQ=17)
                                                                 tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGtcattcat                                                                        >  2:1786489/1‑102 (MQ=25)
                                                                       aggtcaaaCCCGCCGTAGCGGGTTTTTACGTAAAACAGGTCAAACTGACCGATCAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCa                  <  1:1835349/145‑1 (MQ=255)
                                                                        gTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTgg                                                                                    >  2:158402/1‑83 (MQ=17)
                                                                             aaCCCGCCGTATCGGTTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGaa            <  1:766049/150‑1 (MQ=11)
                                                                               cccGCCGTATCTGTTTTTTACGTAAATCAGGTTAAACTGACTGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCAAAAAGCATAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAc          <  2:1741421/150‑1 (MQ=11)
                                                                               cccGCCGTATCGGTTTTTTACGTAAATCAGGTTTAACTGACCGATAAGCCGCTTTCTTTTGTGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAATCTAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAc          <  2:407769/150‑1 (MQ=11)
                                                                               cccGCCGTATCGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAc          <  2:1444630/150‑1 (MQ=255)
                                                                               cccGCCGCACCGGTTTTTTACGTAAAACAGGTCAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAc          <  1:125876/150‑1 (MQ=255)
                                                                               cccGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgc                                          >  1:519478/1‑118 (MQ=35)
                                                                                  gccgAATCGTTTTTTTACGTTAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGa       <  1:1972357/150‑1 (MQ=11)
                                                                                  gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgcc                                         >  1:1025286/1‑116 (MQ=35)
                                                                                  gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccgcc                                         <  2:1025286/116‑1 (MQ=35)
                                                                                    cgtatcGGTTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGTTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc     <  2:571167/145‑1 (MQ=255)
                                                                                    cgtatcGGTTTTTTACGTAAAACAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc     <  2:1744328/145‑1 (MQ=255)
                                                                                       atcgtgTTTTTACGTAAATCAGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGAt  <  1:511976/145‑1 (MQ=255)
                                                                                                   cGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGAGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttttgc                          >  1:1464792/1‑112 (MQ=11)
                                                                                                           aGGTTAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                           <  2:1339651/105‑1 (MQ=17)
                                                                                                                 aaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccg                                           >  1:1161930/1‑83 (MQ=33)
                                                                                                                 aaCTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACccg                                           <  2:1161930/83‑1 (MQ=33)
                                                                                                                     gACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAAcc                   <  2:1120482/103‑1 (MQ=35)
                                                                                                                      aCCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAc          <  2:208105/111‑1 (MQ=35)
                                                                                                                      aCCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAAc          >  1:208105/1‑111 (MQ=35)
                                                                                                                         gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCAGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGAcc     <  2:2024103/113‑1 (MQ=25)
                                                                                                                          aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                                   <  2:2060946/82‑1 (MQ=33)
                                                                                                                          aTAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCg                                   >  1:2060946/1‑82 (MQ=33)
                                                                                                                              gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                  <  2:1474527/79‑1 (MQ=34)
                                                                                                                              gCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                  >  1:1474527/1‑79 (MQ=34)
                                                                                                                                                   |                                                                                         
AACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGAT  >  NC_000913/3269883‑3270119

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: