Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,214 T→C 9.0% intergenic (+612/+2) 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,2140TC 9.0% 112.2 / 4.0 45intergenic (+612/+2)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base T (20/20);  new base C (2/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

AAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTT  >  NC_000913/3270072‑3270361
                                                                                                                                              |                                                                                                                                                   
aaaCCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTg                                                                                                                                               <  1:3783/149‑1 (MQ=17)
                                   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                                                                                                               <  4:179218/105‑1 (MQ=255)
                                                                            tGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGa                                                                                                               >  3:179218/1‑105 (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                                                                                                                                          >  1:229237/1‑50 (MQ=38)
                                                                                                        cAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACt                                                                                                                                          <  2:229237/50‑1 (MQ=38)
                                                                                                          aaaGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCg                                                                                                                                     <  7:215562/53‑1 (MQ=35)
                                                                                                              taaaaaCCCGCTTCGGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtca                                                                                                         >  6:73216/2‑77 (MQ=18)
                                                                                                              taaaaaCCCGCTTCGGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGtca                                                                                                         <  5:73216/76‑1 (MQ=18)
                                                                                                              cAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGaa                                 >  5:164251/1‑149 (MQ=255)
                                                                                                               aaaaaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAAc                                >  2:87269/1‑149 (MQ=255)
                                                                                                               aaaaaCCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAAc                                >  6:188274/1‑149 (MQ=255)
                                                                                                               aaaaaCCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTtttgggtatagcgtcgtggacagtcattcatctttctgcccctccaaaagtaaaaacccgcttcggcgggtttttgcttttggaggg                                <  4:101478‑M2/149‑88 (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)
                                                                                                                          gAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGc                     <  4:118199/149‑1 (MQ=255)
                                                                                                                            aGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTtttgggtatagcgtcgtggacagtcattcatctttctgcccctccaaaagcaaaaacccgcttcggcgggtttttgcttttggaggggcagaaagatgaa                   >  2:59756‑M2/1‑49 (MQ=255)
                                                                                                                             gCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTtttgggtatagcgtcgtggacagtcattcatctttctgcccctccaaaagtaaaaacccgcttcggcgggtttttgcttttggaggggcagaaagatgaat                  >  3:101478‑M2/1‑48 (MQ=255)
                                                                                                                               gggTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGc                <  4:140763/149‑1 (MQ=255)
                                                                                                                                        cGTAAATCAGGTGAAACTGACCGCTAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTg        <  6:164251/148‑1 (MQ=255)
                                                                                                                                         gTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGa      <  5:162151/149‑1 (MQ=255)
                                                                                                                                          tAAATCAGGTGAAACTGACCGATAAGCCGCGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAg     <  2:162352/149‑1 (MQ=255)
                                                                                                                                            aaTCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGtt   >  7:141620/1‑149 (MQ=255)
                                                                                                                                             aTCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGttt  <  6:241590/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
AAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTT  >  NC_000913/3270072‑3270361

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: