Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,606,903 A→G 26.7% S267P (TCC→CCC)  rsmC ← 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,9030AG26.7% 94.4 / 7.1 58S267P (TCC→CCC) rsmC16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (9/32);  new base G (14/3);  total (23/35)
Fisher's exact test for biased strand distribution p-value = 5.16e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTTAACGTC  >  NC_000913/4606665‑4607132
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
gaTAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCCTCGTGGAACGGCGGGTTGga                                                                                                                                                                                                                                                >  2:69892/1‑239 (MQ=255)
  aatttAACCTTTCTTCCCCGGGGGGGGTCTAAGGGCGGGATCCCCCTTGAAGCGGCCTGTTGGCGCGATCCCTTCGTGGAAGCCAAATGTCTCATCCAGCCCGGCCGGATAAGGCAGGAAGGCGGTCGCTAAAAGTCGCAGCTGGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgag                                                                                                                                                                                                                                               <  1:195090/235‑1 (MQ=255)
          cTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTTCAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCCGGCTGGTTTTCCTCCCAACGGGGAACGGGGGCGTGGGgatgtgcctt                                                                                                                                                                                                                                      >  1:898359/1‑233 (MQ=255)
             tcttcGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCGCCATCCCGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGgatgagtcaatca                                                                                                                                                                                                                                   >  2:59990/1‑231 (MQ=255)
              cttcGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGACCGGCGGGTTGGGGAGTATCatattca                                                                                                                                                                                                                                 >  2:608157/1‑236 (MQ=255)
                cccgctgtggggTCATAATGGCGCGATCCCCCTTGAAGCCCCCTGTTTGCCCGACCACTTCGGGGAAGCCAAATGCCTCATCCAGCCCGTCCGGATAAGGCAGGAAGGCGTCCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCaaaa                                                                                                                                                                                                                               <  1:967400/230‑1 (MQ=255)
                             tttAATGGCGCGTTACCCCTTTAAGCGCCCGTTTTGCGGGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCTTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCAcc                                                                                                                                                                                                                 <  1:98782/236‑1 (MQ=255)
                             taTAATGGGGCGATACTCCTTGAAGCGCCCTGTTTGCGCGGTCCCTTCGGGGAAGCCAAAGGTCTAACCCAGCACGTCCGGATAAGGCAGGAAGGCGGTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCAcc                                                                                                                                                                                                                 <  1:370407/237‑1 (MQ=255)
                              aTAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCAACGGGGAACGGCGGGGTGGGGGAGATACTTATCaaaaaaacattcacccc                                                                                                                                                                                                                  >  1:168104/1‑226 (MQ=255)
                                     cgcgATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGAGTTTGCGCCCCATCCAGGCTGGTTTGCATCCCCTCGGGGAACGCCGGGGTGGGGGTGATCAGATCaaaaacacacttcccccccggaaag                                                                                                                                                                                                          >  1:987714/1‑218 (MQ=255)
                                           accccTTGAAGCCCCCTGTTTGCGCGATCGCTTCGGGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTg                                                                                                                                                                                                 <  2:999406/236‑1 (MQ=255)
                                            ccccTTGAAGCGCCCTGTTTGCGCGATCACTTCGGGGAAGCCAAATGTCTCATCCAGCCCGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGc                                                                                                                                                                                                <  2:898080/237‑1 (MQ=255)
                                                 tggaccgcctcttttGCGCCATCCCTTCGTGGAACCCAAATTTCTCATCCAGCACGCCCGGATATGGATGGAAGGCGTTCGCGACAATTCCCAGCTCGCCGGCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCGCCTCGGAAAAGACGTTGCTGGCa                                                                                                                                                                                          <  1:377598/227‑1 (MQ=255)
                                                   aaGCGCCCTGTTTGCGGGTTCACTTCGGGGAACGCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATGCGCAGCTCGCCGCCGCTATTAAGATGACGCCCCGCACCGCGAATCAGGGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCaaa                                                                                                                                                                                        <  1:220485/239‑1 (MQ=255)
                                                        cccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCAACCCATCGTGGAACGGCGGGTTGGAGATGATCACATCACAACGGCCTTTCACCCCGGGAAAGACGGTGTTGGCcaaaaac                                                                                                                                                                                    >  2:953450/1‑234 (MQ=255)
                                                         ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCCCCGCGAATCAGCGTTTGCCCCGCATCCAGGCTGGTTTGCATCCCATCGTGGACCGGCGGGTGTGGGATGATTATTTCAACAAGACCCTTCCCCCCCGGaaaaaaggtgccgggaaaaacat                                                                                                                                                                                      >  1:1016778/1‑219 (MQ=255)
                                                         ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGATCCTATCCAAAAGACCTTTCCCCCCGGGAAAGACGTTGCTGGCAAAAAATTc                                                                                                                                                                                  >  2:254971/1‑239 (MQ=255)
                                                         ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGCTGGAGATGATCATTTCAAAACGACACTTTACCCCGGGAAAAACGATGCTGGCAAAGACCTc                                                                                                                                                                                  >  2:495453/1‑239 (MQ=255)
                                                           tttttgggggttaccttgggggcaGCCAAAGGTCTCATCCAGCACTTCCGGATAAGGCGGGAAGGCGTTCGCTCCGATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCAc                                                                                                                                                                               <  1:528439/216‑1 (MQ=255)
                                                           tttttGGGCGATTCCTTTGGGAAGGCAAAATGTTTCAGCCAGCACGCCCGGAAAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCAc                                                                                                                                                                               <  1:648191/237‑1 (MQ=255)
                                                                     aTCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATTATCATTTCACAACGACCTTTTCCCCCGGAAAAAACGTTGCTGGGAAAGACTTACCCTTCCacac                                                                                                                                                                     >  2:388093/1‑238 (MQ=255)
                                                                       cACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgag                                                                                                                                                                                                                                               >  2:416620/1‑169 (MQ=255)
                                                                       cACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgag                                                                                                                                                                                                                                               <  1:416620/169‑1 (MQ=255)
                                                                           tggTGGAAGCCAAATTTCCAATCCACGCCGTCCGGATAAGGAAGGAAGGGGTTCGCTGCAATTCGCAGTCCCCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGc                                                                                                                                                              <  2:748183/237‑1 (MQ=255)
                                                                              tGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGCTTGCATCCCATCGTGGAACGGCGGGGTGGGGCGGTTACTATCAAAACCACCCTTCACCCCCGACACAACGCTGCGTGGCAAAACCTCACCCCTCAAAACCCTTGGcgc                                                                                                                                                            >  1:956637/1‑237 (MQ=255)
                                                                                      aaaTGTCTCAGCCAGCACGGCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTg                                                                                                                                                   <  2:701680/239‑1 (MQ=255)
                                                                                       aaTGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGGTGATCATACCAAAAACACCCTTCTACTCGGAAAAGAAGTTGATGGGAAAGACTTCAACCTCAACAACggcgaccgccccagtagg                                                                                                                                                  >  1:844909/1‑222 (MQ=255)
                                                                                             tcATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGTGTGGGCGTGAACCTATTAAAAAAACCCTTCCACCCCGAAAAAGAGTGTCGCGCACAAACCTCACCCTCCACACCCATGGCGGCaaaggggggcgggcg                                                                                                                                            >  1:864438/1‑226 (MQ=255)
                                                                                                tCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCt                                                                                                                                                                                                                <  2:980437/172‑1 (MQ=255)
                                                                                                    gcccGTCCGGATAAGGCAGGAATGCGTTCGCTACAATTCCGATCCAGCCGCCGCTATTAAGATGCCCCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCt                                                                                                                                     <  2:671176/236‑1 (MQ=255)
                                                                                                     cACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGATTCTATCAAAACGCCCTTTCACCTCGGGAAAGACGTTGCTGGCAAAAACTTCACCTTCAAAACCCTTGGGCGCAAAGGATGCGCGGCTGGGTTCt                                                                                                                                     >  2:114369/1‑238 (MQ=255)
                                                                                                         ttcggaaagggcgaggaGGCGGTCGCAACAATTCGTACGCCGCCGCCGCTATTAAGATGACGTACCGCCCCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCAATCGAGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGGTGCTGGCAAAGACTTCACCTGCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTAccgc                                                                                                                                <  1:608157/223‑1 (MQ=255)
                                                                                                          ccGGATAAGGCAGGAAGGCTTTCCCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTAccgcc                                                                                                                               <  2:168104/239‑1 (MQ=255)
                                                                                                                 aGGCAGGAAGGGGTTCGCTACAATTCGAAGCCCCCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCag                                                                                                                       <  2:924349/239‑1 (MQ=255)
                                                                                                                  ggCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGATCGATTCAAAACACACTTTCCACCTGGAAAAGAAGGTTGCTGCGAAAACTTCCCCTTCCACACCACTGGCCGCCAGAGTGGGGCgggtggccatcaccccccgggcaag                                                                                                                       >  1:306598/1‑221 (MQ=255)
                                                                                                                       gAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGGTGCTGGCAAAGACTTCACCTTCAAAACCGTTGGCCGCAAGTGTTGCGCGGCTGGGTTTTACCGCCCGGGCCGAGACa                                                                                                                  >  2:207327/1‑238 (MQ=255)
                                                                                                                          ggCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGAGAGGATCCTCTCAAAAAGACCCTTCAACCCCGAAAAAAAGTTGGTGGCAAAGACTTCAACTCCAAAACCCTTGTCCGCACAGGATTCGCGGCGGGGTTCTCCCCCCcgggggacgaccaccc                                                                                                               >  1:668236/1‑227 (MQ=255)
                                                                                                                            cggtCGCTACAATTCGCAGCCCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGca                                                                                                             <  2:440127/235‑1 (MQ=255)
                                                                                                                            cGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTc                                                                                                                                                                                  <  2:713659/172‑1 (MQ=255)
                                                                                                                            cGTTCCCTAAAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGca                                                                                                             <  2:956637/238‑1 (MQ=255)
                                                                                                                              ttCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGATGATCATATCAAAACGACCTTTCACCTCGGCAAAGACGTTGCTGGGAAAGACTTCCACCTCAACACCATTGGCCGCAAGGGTTTGCCGGCTGGCCTCTACCCCCGGGGCAGAGACaacccgaaa                                                                                                          >  1:826954/1‑231 (MQ=255)
                                                                                                                                 gCTACAATTCGCCGCTGCCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGt                                                                                                      <  2:826954/239‑1 (MQ=255)
                                                                                                                                                                  tcccgccccgcccCCCGAAGCACGGTTGGCGCCGCAGCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATAACAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGCCTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAgg                                                                      <  1:114369/235‑1 (MQ=255)
                                                                                                                                                                    aCGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGAGATGATCATATCAAAAAGACCTTTCACCCCGGAAAAGACGTTGCTGGCCAAGACTTCACCTTCAACACCGTTGGCCGCAAGGGTTGCGCGGCGGGCTTCTACCGCCCGGGCAAAGAAACCGCCCAAAGGGAGACGAAGTTTTGGCGAATGGCGGGCGAAGGCaa                                                                   >  1:225453/1‑239 (MQ=255)
                                                                                                                                                                                 cgAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGATGATCATATCCAAACGACCTTTCACCTCCGGAAAGACGTTGCTTGCAAAGAATTCACCTTCACCCCCCTTTGCCCGCAAGGGTTTGCGGCTGGCTTCTAACCACCGGGCAGGGGCAACGCACAAGGGGGAGCGAGATTTTGGCGAACAGCGGGCaaaaaacaccgaaaaaaaccc                                                      >  1:596915/1‑221 (MQ=255)
                                                                                                                                                                                  gAACCAGGGTTGGCGCCGCATCCGGGCTGGTTTGCACCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGcc                                                      <  2:987714/238‑1 (MQ=255)
                                                                                                                                                                                         ggTTTGCGCCTCAACCGGGCTGGGTTGCCTCCCATCGTGGACCGGGGGGTGGGGGATGGACATATCAAAACGGCCTTTCCCCTGGGAAAAGACGTTGCTGGCAAAGACTTCACCTCCAACACCGTGGGCCGCAAGTGTTGCGCGGCTGGCTTCTCCCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCcac                                              <  1:897351/238‑1 (MQ=255)
                                                                                                                                                                                         cGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAAACGTTGCTGGCAAAGACTTCACCTTCAAAACCGTTGGCCGCAAGGGTTGCGCGGCTGGCTTCTACCGCCCGGGCAGAGACATCGCACAAGGGGGGACGAATTTTCGGCGAATGGCGCGCCAAGGCAACTTAAAgaaacccccgccccc                                              >  2:756928/1‑225 (MQ=255)
                                                                                                                                                                                              tgccccgctCCGGGGTGGTTTGCATCCAAGCGTGGAAGGGGGGGTTGGGGATGATAAAATCAAAACGACCTTTCCCCTCGGAAAAGACGTGGCGGGCAAAGCCTTCCCCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATTGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCg                                         <  1:388093/231‑1 (MQ=255)
                                                                                                                                                                                              gcgcCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACTGCCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCg                                         <  2:306598/239‑1 (MQ=255)
                                                                                                                                                                                                   gcATCCAGGCTGGTTTGCATCCCATCGTGGAACGGGGGGTTGGAGATGATCAGATCAAAACGACCTTGCACCTCGGAAAAGACGTTGCTGGCAAAGACTTGCCCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATc                                    <  2:919470/239‑1 (MQ=255)
                                                                                                                                                                                                     aaccAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATc                                    <  2:5702/235‑1 (MQ=255)
                                                                                                                                                                                                     aTCCGGGCTGGTTTGCACCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCa                                  <  2:596915/239‑1 (MQ=255)
                                                                                                                                                                                                            cttgTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATACCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACtt                            <  2:646286/235‑1 (MQ=255)
                                                                                                                                                                                                               gTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGAAGGTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGGTGCGCGGCTGGCTTCTCCCGCCGGGGCAGGAGACTCGCCCAAAGGGGAGCGAAATTTCGGCGACAGGGCCGCAAACGGAACCGAAAACAGCCCCGCCCCCCAccccgacaccacccccttac                         >  1:602453/1‑221 (MQ=255)
                                                                                                                                                                                                                 ttGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACCCCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCt                       <  2:579132/238‑1 (MQ=255)
                                                                                                                                                                                                                        ccATCGTGGAACGGGGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGtgtg               <  2:1016778/239‑1 (MQ=255)
                                                                                                                                                                                                                         catgggggAAGGGGGGGTTGGGAATGAACATATCAAAACGGCCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCCCCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGc              <  2:429674/233‑1 (MQ=255)
                                                                                                                                                                                                                         cATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGc              <  2:225453/239‑1 (MQ=255)
                                                                                                                                                                                                                                      cggGTTGGAGATGATCATAGCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAACACTTCACCTACAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTGTCGTGTGCGGAGTTAACGTc  <  2:237373/238‑1 (MQ=255)
                                                                                                                                                                                                                                        ggTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCag                                                                                                                       <  2:25323/120‑1 (MQ=255)
                                                                                                                                                                                                                                        ggTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCag                                                                                                                       >  1:25323/1‑120 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
GATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTTAACGTC  >  NC_000913/4606665‑4607132

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: