Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,193 C→T 12.6% intergenic (+591/+23) 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,1930CT12.6% 93.5 / 5.0 40intergenic (+591/+23)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base C (20/14);  new base T (2/3);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTG  >  NC_000913/3270048‑3270340
                                                                                                                                                 |                                                                                                                                                   
atCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAACCTGACCGATAAGCCGCTTTCTTTTGGGGAGAGCGTATTGGCAAGTCATTCAAGGGATCGCACCGCCAACAGCAAAAAAccgccg                                                                                                                                                    >  5:241590/1‑147 (MQ=11)
                 aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                                       <  8:6450/143‑1 (MQ=37)
                 aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTa                                                                                                                                       >  7:6450/1‑143 (MQ=37)
                      aaaaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTGCttttggaggg                                                                                                                            >  4:101478‑M1/1‑139 (MQ=255)
                        aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTg                                                                                                                          <  1:3783/149‑1 (MQ=17)
                                   aGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGTTTTTGCttttggaggggcagaaagatgaa                                                                                                               <  2:59756‑M1/149‑24 (MQ=255)
                                    gCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGTTTTTGCttttggaggggcagaaagatgaat                                                                                                              <  3:101478‑M1/149‑25 (MQ=255)
                                                          tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                                                                                                             <  7:110574/96‑1 (MQ=33)
                                                          tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                                                                                                             >  8:110574/1‑96 (MQ=33)
                                                           gAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGCGAGTCATTCATCGTTGGGCCGCGGCAATAGCAAAAACCCGCCGAAGCGGGGGTTTACGTATATGAGGTGAAAGTGGCCGATAGGCCGGGGTGTGTCGTgggctg                                                                                       >  4:19912/1‑145 (MQ=9)
                                                                     gATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCgg                                                                                                       <  2:230671/123‑1 (MQ=11)
                                                                          gCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGc                                                                         >  5:102051/1‑148 (MQ=255)
                                                                                  ctttTGTGCATACCGGCGCTGACAGTACTTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATc                                                                 <  7:188584/148‑1 (MQ=18)
                                                                                     ttGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCt                                                              >  1:246781/1‑148 (MQ=255)
                                                                                         gTATAGCGTCGTGGACAGTCATTCATCTTTCAGCCCCACCAAAAACAAAAACCCCCCGCAGCGGGTTTTTACGTAAATCACGTTAAACTGCCCCATAAGCCGGGTAATGTCCTGGACAGTCATTAATCTTGGCCCGCAAACGCtacaat                                                         >  6:231589/1‑144 (MQ=17)
                                                                                                   gtGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGtt                                                                                                    >  5:5038/1‑96 (MQ=38)
                                                                                                   gtGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGtt                                                                                                    <  6:5038/96‑1 (MQ=38)
                                                                                                    tGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGa                                                                                          >  3:179218/1‑105 (MQ=255)
                                                                                                    tGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGa                                                                                          <  4:179218/105‑1 (MQ=255)
                                                                                                         aGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTAcc                                         >  6:162151/1‑149 (MQ=255)
                                                                                                             attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt                                                                                                >  7:204834/1‑90 (MQ=255)
                                                                                                             attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGt                                                                                                <  8:204834/90‑1 (MQ=255)
                                                                                                             attcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCggg                                     >  2:201927/1‑149 (MQ=255)
                                                                                                              ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGtt                                                                                                    <  7:214051/85‑1 (MQ=39)
                                                                                                              ttcatCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGtt                                                                                                    >  8:214051/1‑85 (MQ=39)
                                                                                                                   cTTTCTGCGCATCCATAAGCATAAACCCGTAGACGAGTGTTTTTACGTAAATGAGGCGAATCTGACCGCTATGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTc                                  <  3:139199/146‑1 (MQ=255)
                                                                                                                    tttCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTa                              >  3:118199/1‑149 (MQ=255)
                                                                                                                         gCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGAGCGATAAGCCGGGTTCTGTCGGGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGc                         >  3:174590/1‑149 (MQ=255)
                                                                                                                           cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGGGTTCTGTAGTgaa                                                                                          >  8:9680/1‑80 (MQ=11)
                                                                                                                             cTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTa                     >  1:162352/1‑149 (MQ=255)
                                                                                                                              tCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTtttgggtatagcgtcgtggacagtcattcatctttctgcccctccaaaagcaaaaacccgcttcggcgggtttttgct                    <  1:59756‑M2/149‑79 (MQ=255)
                                                                                                                                cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACt                                                                                                                     <  2:229237/50‑1 (MQ=38)
                                                                                                                                cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACt                                                                                                                     >  1:229237/1‑50 (MQ=38)
                                                                                                                                  aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCg                                                                                                                <  7:215562/53‑1 (MQ=35)
                                                                                                                                    aGCAAAAACCCGCCGAAGCGGGTTTTTGCttttggaggggcagaaagatgaatgactgtccacgacgctatacccaaaagaaagcggcttatcggtcagtttcacctggtttacgtaaaaacccgcttcggcgggtttttgcttttgga              <  1:59756‑M1/149‑121 (MQ=255)
                                                                                                                                      taaaaaCCCGCTTCGGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtca                                                                                    <  5:73216/76‑1 (MQ=18)
                                                                                                                                      taaaaaCCCGCTTCGGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtca                                                                                    >  6:73216/2‑77 (MQ=18)
                                                                                                                                      cAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGaa            >  5:164251/1‑149 (MQ=255)
                                                                                                                                       aaaaaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAAc           >  6:188274/1‑149 (MQ=255)
                                                                                                                                       aaaaaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAAc           >  2:87269/1‑149 (MQ=255)
                                                                                                                                        aaaaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGa                                                                                                               <  7:111397/48‑1 (MQ=38)
                                                                                                                                              cgccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAtt    <  1:194694/149‑1 (MQ=255)
                                                                                                                                               gccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTAttt   >  4:146165/1‑149 (MQ=255)
                                                                                                                                                ccgAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTg  <  1:117495/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTG  >  NC_000913/3270048‑3270340

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: