Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,270,196 A→G 11.4% intergenic (+594/+20) 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,1960AG11.4% 103.6 / 6.6 44intergenic (+594/+20)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base A (21/17);  new base G (2/3);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

AAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGC  >  NC_000913/3270065‑3270342
                                                                                                                                   |                                                                                                                                                  
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                                                                                                                                               >  8:110574/1‑96 (MQ=33)
                                         tGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCggg                                                                                                                                               <  7:110574/96‑1 (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                                                                                            <  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                                                                                                                       <  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)
                                                                                                                                 gAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGc  <  4:118199/149‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                                  
AAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAATCAGGTGAAACTGACCGATAAGCCGGGTTCTGTCGTGGACAGTCATTCATCTAGGCCAGCAATCGCTCACTGGCTCAAGCAGCCTACCCGGGTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGC  >  NC_000913/3270065‑3270342

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: