Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,921,544 T→G 17.6% V441G (GTA→GGA)  rsmF → 16S rRNA m(5)C1407 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,921,5440TG17.6% 188.7 / 4.4 85V441G (GTA→GGA) rsmF16S rRNA m(5)C1407 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (18/50);  new base G (15/2);  total (33/52)
Fisher's exact test for biased strand distribution p-value = 4.50e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGCCTCAACTAGGCTGAAAAATGGTGCGATCGGACTGGTCGTACCACAATCGGCAGCTAAATGGA  >  NC_000913/1921320‑1921769
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                        
gttCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTAtt                                                                                                                                                                                                                                       >  1:1653450/1‑228 (MQ=255)
gttCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGAGGTTTAACCGCAAGCCGAGCCAGTGTCGGGAGATGGAtt                                                                                                                                                                                                                                       >  2:860332/1‑228 (MQ=255)
  tCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTa                                                                                                                                                                                                                                         <  2:117651/224‑1 (MQ=255)
   cccGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGt                                                                                                                                                                                                                                    >  1:269479/1‑228 (MQ=255)
     cGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTa                                                                                                                                                                                                                                  <  1:1476015/228‑1 (MQ=255)
           gCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACGAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGGAGGTGGATTGGGTACttttc                                                                                                                                                                                                                            >  1:66998/1‑226 (MQ=255)
           gCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTcc                                                                                                                                                                                                                            <  2:688922/228‑1 (MQ=255)
              ttGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGAGGGTTTGGTTACTTTTCAGc                                                                                                                                                                                                                        >  1:1621308/1‑228 (MQ=255)
                             gTAAAGTCCGATTTGCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGtt                                                                                                                                                                                                         <  1:1907563/228‑1 (MQ=255)
                               aaaGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCCAGCCGCGCCAGTGGCGGATGCTGTATTGGGGACTTTTCAGCATCACCCGACTGGTGTa                                                                                                                                                                                                       >  1:1111526/1‑228 (MQ=255)
                                 atgccGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACCCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGcc                                                                                                                                                                                                    <  2:1242548/225‑1 (MQ=255)
                                   tCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGGTGGTGTATTGGTTACTTTCTCGCGTCAACCGCTTGGtgtttccc                                                                                                                                                                                                   >  2:696077/1‑221 (MQ=255)
                                   tCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGtt                                                                                                                                                                                                                                   >  1:1239680/1‑197 (MQ=255)
                                    ccGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTa                                                                                                                                                                                                                                  <  2:1237466/197‑1 (MQ=255)
                                       aTTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGt                                                                                                                                                                                                                                          >  1:1650274/1‑186 (MQ=255)
                                        ttttCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTa                                                                                                                                                                                                                                         <  2:1647355/186‑1 (MQ=255)
                                             tcGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGGATTGGGTTCTTTTCAGGCTTCAACCACTCGGTTTTCCCaacagg                                                                                                                                                                                              >  2:1189816/1‑220 (MQ=255)
                                              cGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGa                                                                                                                                                                                             <  1:1191950/223‑1 (MQ=255)
                                                tttGGGGATTAAACTTGCCGAGCCGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTGGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGtt                                                                                                                                                                                       <  2:1113095/226‑1 (MQ=255)
                                                  tGGGGATTAAACTTGCCGAACCGCACAACAAAGGTTAACGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTcgc                                                                                                                                                                                    <  1:697246/228‑1 (MQ=255)
                                                     gtatttatcTTGCCGATACGCACAACAAAGGTGATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAAGATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGGGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGAt                                                                                                                                                                                <  1:1947936/220‑1 (MQ=255)
                                                       aTTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTg                                                                                                                                                                              <  1:861837/228‑1 (MQ=255)
                                                         tAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGGATTGGTTACTTTCCAGCAGCAACCGATTGGGTTAGGCCAACGGATTGGTTCGCGAttta                                                                                                                                                                             >  1:1385662/1‑225 (MQ=255)
                                                          aaaCTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGaaa                                                                                                                                                                           <  1:768082/228‑1 (MQ=255)
                                                                 ccGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTCCCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTa                                                                                                                                                                                                       <  1:760424/194‑1 (MQ=255)
                                                                  cGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGAGGATGGATTGGTTACTTTTCAGCATCAACCGATTGGTTTTGCCAAAAGGAGTGGTTTGCGGATTTaaaaaaacca                                                                                                                                                                    >  2:288689/1‑223 (MQ=255)
                                                                   gagaCGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTa                                                                                                                                                                   <  2:1650525/227‑1 (MQ=255)
                                                                   gagaCGAACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTAt                                                                                                                                                                  <  1:1584716/228‑1 (MQ=255)
                                                                     tgagcgcaacaaAGGTTATCGCTGGCAGAATGAAGCGGTTAGTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATcc                                                                                                                                                                <  1:1339421/222‑1 (MQ=255)
                                                                             acaaAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGaa                                                                                                                                                        <  2:269013/227‑1 (MQ=255)
                                                                             acaaAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGGATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGaa                                                                                                                                                        <  1:1832772/227‑1 (MQ=255)
                                                                                        cGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGGATTGGTTACTTTTCCGCATCCACCGAATGGGTTAGCCAAACGGGTTGGGTCGCGATTGCAAAACAGCTATTCCCGGGAACTGGTGCGCGAt                                                                                                                                           >  1:1308519/1‑228 (MQ=255)
                                                                                              cagcaTGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGAGTGTGGATTGGTTAATTTTCa                                                                                                                                                                                                                           >  2:1955387/1‑146 (MQ=255)
                                                                                               agcaTGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAg                                                                                                                                                                                                                         <  1:1958875/146‑1 (MQ=255)
                                                                                                                     tGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATTGTGGATTTGTTTATTTCCCGCATCAACCCCTTGGTTTTGGCAAACCGACTGGGTTGCGATTGAAAAACAGCCATCCGCGCTCACTGGGGGGCCGTGGGGAACTTCTTTCCGTGTACGCCTGGCa                                                                                                           >  2:31596/1‑228 (MQ=255)
                                                                                                                     tGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTGTTGGGTTCTTTCCAGCATCAACCGATTGGTTTTGCCCAAAGGATTTGGTTGCGCTTGTAAAACAGCTATCCGCGTGAACTGGGGCGCGCTGGGGAACCTGTTTACGGTCAAGCCTGAct                                                                                                            >  2:1589550/1‑227 (MQ=255)
                                                                                                                     tGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTTCTTTCCCCCATCAACCCATTGGGTTAGCCAAAAGGATTGGGTCGCGGTTTAAAAAAAGCAATCCCCGCGAAACTGGGTGGGCGTGGGAACATTTTTCCCGTTAAGCCTGGCa                                                                                                           >  2:950779/1‑228 (MQ=255)
                                                                                                                      gCCTCAACCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTa                                                                                                                       <  1:952466/217‑1 (MQ=255)
                                                                                                                        cacACCCGATAATATGAATGCATTCGAGCGGCCACCGCAGGAAGCGGAGGAGTGGTATCGCGGGAGCGATGTTTACCCGCAAGCCGCGCCAGTGGAGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAgcg                                                                                                        <  1:1762830/226‑1 (MQ=255)
                                                                                                                                ataataTGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGGATTGGTTACTTTCCAGCATAAACCGATTGGGTTAGACAAACGGATTGGTTTGCGATTGAAAAACAGCTATCCGCGGGAAATGGGGGGCGATGGGGAAAATTTTTAcg                                                                                                                           >  2:1765688/1‑203 (MQ=255)
                                                                                                                                 taataTTAATGCATTCGAGCTGACGCCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCgg                                                                                                                         <  1:1768812/204‑1 (MQ=255)
                                                                                                                                 taataTGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAg                                                                                                                                                                                                                         >  2:1424050/1‑112 (MQ=255)
                                                                                                                                  aataTGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGc                                                                                                                                                                                                                        <  1:1426606/112‑1 (MQ=255)
                                                                                                                                       gAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTgg                                                                                                                                                                                                                                     >  2:1502330/1‑95 (MQ=255)
                                                                                                                                       gAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGGATTGGTTACTTTTCAGCATCAACCGATTGGGTTATGCAAACGGATTGGTTCTCGATTGAAAAACAGCTa                                                                                                                                                                   >  2:1944735/1‑159 (MQ=255)
                                                                                                                                        aaTGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGt                                                                                                                                                                                                                                    <  1:1505015/95‑1 (MQ=255)
                                                                                                                                        aaTGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTAt                                                                                                                                                                  <  1:1948193/159‑1 (MQ=255)
                                                                                                                                         aTGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGCTGATGGATTGGGTACTTTTCAGCATCAACCCATTGGTTTAGCCCAACGGATTGGTTCGCGATTGGAAAAAAACCATCCCCGGGAACTGGGTGGCCGTGGCGAACTTTTTACCGGTTATGCCTGGCAAGCGACCAAAAAAgcggca                                                                                        >  2:2010955/1‑224 (MQ=255)
                                                                                                                                                  ggCTGACCCCGAAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTg                                                                              <  2:1618445/227‑1 (MQ=255)
                                                                                                                                                       acacCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCCATTGGGTTTGCCCAAAGGGTTGGTTTGCGGTTGGAAAACAAGCATCCCCGTTACCCGGGGGGCGGGTGGAAACCTTTTTCCGGGAAAGCCCGAACAGGCCCAAAAAAAACGCAACTTTTgtattggac                                                                          >  2:17102/1‑219 (MQ=255)
                                                                                                                                                        cacCGGAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGcaca                                                                         <  1:675343/227‑1 (MQ=255)
                                                                                                                                                            ggaGGAAGCGGGGGAGTGGTATCGGGGGCGCGATTTTTACTCGGAAGCCGCGCCAGGGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCCATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCgg                                                                    <  2:1383177/226‑1 (MQ=255)
                                                                                                                                                            ggaGGAAGCGGAGGATTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCgg                                                                    <  1:289182/226‑1 (MQ=255)
                                                                                                                                                            gaaGGAAGCGTAGGGGTGGTATCGCGGGCGCGATGTTTCCCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCgg                                                                    <  2:652655/226‑1 (MQ=255)
                                                                                                                                                            gCAGGGAGGGGAGGAGTGGTAGCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCgg                                                                    <  1:617024/228‑1 (MQ=255)
                                                                                                                                                             cAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGGGCGATGGGGAACTTTTTTCCCGTAACGCCCGACAACGCCCCAAAAAAGGGGACCTTTTGACTTGCACCATCGGc                                                                   >  2:1644947/1‑228 (MQ=255)
                                                                                                                                                                 aaGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGc                                                                <  1:17136/227‑1 (MQ=255)
                                                                                                                                                                     ggaggaGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTgg                                                                                                                                                                                                           >  2:1689379/1‑90 (MQ=255)
                                                                                                                                                                      gaggagTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGt                                                                                                                                                                                                          <  1:1692368/90‑1 (MQ=255)
                                                                                                                                                                            tGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGCCTCAACTAgg                                                      <  2:1113776/226‑1 (MQ=255)
                                                                                                                                                                              gTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCcgc                                                                                                                                                              <  1:1081677/125‑1 (MQ=255)
                                                                                                                                                                                    cgGTCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGCCTCAACTAGGCTGAAAAAt                                             <  1:1277451/227‑1 (MQ=255)
                                                                                                                                                                                    cgGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGGATTGGTTAATTTTCCGCATCAACCGATTGGTTTAGCCCAAAGGATTGGTTTGCGATTTAAAAACAGGTATTCGCGTGAAATGGGGGGGGGTGGGGAAATTTTTAACggg                                                                                                                        >  2:1051050/1‑153 (MQ=255)
                                                                                                                                                                                     gggCGCGATGTTTACCCGCAAGCCGCGACAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTa                                                                                                                       <  1:1052943/154‑1 (MQ=255)
                                                                                                                                                                                      ggCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGGATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAAAGGGTTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGNNNNNCGATGGGAAACTTTTTACCGGTAACGCCTGACAgc                                                                                                         >  2:1836424/1‑165 (MQ=255)
                                                                                                                                                                                      ggCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGAAGGTGTAATGGGTTCTTTCCCGCATCAACCGATTGGGTTAGCCAAAAGGGTTGGTTCGCGATTGAAAAACAGATATCCGCGGGGACTGGGGCGGGATGGGGAAATTTTTTCCGGTAACGGCTGGCAGCGCACCAAAAAAACGCACCTTTTTACTGGCACATTCGTCTGCCTCAACTaaggcggcaaat                                             >  2:701508/1‑214 (MQ=255)
                                                                                                                                                                                       tcagggtGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTNNNCCNGTAACGCCTGACAgcg                                                                                                        <  1:1839683/159‑1 (MQ=255)
                                                                                                                                                                                        cgcgATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTTCAGCAACAACCGAATGGTTTAGGCCAACGGATTGGTTCGCGATTGAAAAACAGCCATCCGCGGGAAATGGGGCGGGAGGGGAAACGTTTTACCGGTAAAGCCTGACAACGCACAAAAAAAGAGCCCtt                                                                                      >  2:419598/1‑182 (MQ=255)
                                                                                                                                                                                         gcgATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACttt                                                                                     <  1:420300/182‑1 (MQ=255)
                                                                                                                                                                                         gcgATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGAACATTCgg                                                                    <  1:1843378/199‑1 (MQ=255)
                                                                                                                                                                                            ttGTTTACCCGAAGGCCGGGACGGGTGCGGAGTATTTATTTGTTACTTTCCCGCACTAACGTATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGCCTCAACTAGGCTGAAAAATGGTGCGAt                                     <  2:1570463/226‑1 (MQ=255)
                                                                                                                                                                                                ttACCCGCAAGCCGTGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGCCTCAACTAGGCTGAAAAATGGTGCGATCGGAc                                <  2:1557658/228‑1 (MQ=255)
                                                                                                                                                                                                ttACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGCCTCAACTAGGCTGAAAAATGGTGCGATCGGa                                 <  1:1301697/227‑1 (MQ=255)
                                                                                                                                                                                                aagcccGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAACCAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGCCTCAACTAGGCTGAAAAATGGTGCGATCGGAc                                <  1:1733452/225‑1 (MQ=255)
                                                                                                                                                                                                       cAAGCCGCGCCAGTGGCGGATGATGTATTGGTTAATTTCCAGCATCAACCGATTGGTTTAGCCAAACGGGTTGGTTTGCGATTGAAAAAAAGCTATCCGCGTGAACTGGTGCGCGGTGGGAAACTTTTTTCCg                                                                                                                          >  2:1212123/1‑133 (MQ=255)
                                                                                                                                                                                                        aaGGCGCGCCAGTGGTGGGTGATGTATTGGTTACTTGCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGGGATGGGAAACTTTTGACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGCCTCAACTAGGCTGAAAAATGGTGCGATCGGACTgg                             <  1:834045/223‑1 (MQ=255)
                                                                                                                                                                                                        aaGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCgg                                                                                                                         <  1:1214291/133‑1 (MQ=255)
                                                                                                                                                                                                           ccGCGCCAGTGGCGGATTATGTATTGGTTACTTTCCAGCAACAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGCCTCAACTAGGCTGAAAAATGGTGCGATCGGACTGGTCGTACca                     <  2:1597586/228‑1 (MQ=255)
                                                                                                                                                                                                              cgcCAGTGGCGGATGATGTATTGGGTACTTTTCAGCATCAACCGATTGGGTTAGCCCAACGGATTTGGTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGGGCGCGGTGGGAAAAttttt                                                                                                                              >  1:218794/1‑122 (MQ=255)
                                                                                                                                                                                                              cgcCAGTGGCGGATGATGGATTGGTTACTTTCCAGCATCAACCCATTGGTTTAGCCAAAAGGATTGGTTCGCGGTTGAAAAAAAGCTATCCGCGTGAACCGGGGCGCGATGGGGAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAACGCACTTTTTGACTGGCACATTCGGCTGCCTCAAATAGGGTGAAAAATGGTGTGGTCGGGCCGGTCggaacccaa                  >  2:571401/1‑220 (MQ=255)
                                                                                                                                                                                                               gcCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTa                                                                                                                             <  2:218421/122‑1 (MQ=255)
                                                                                                                                                                                                                  aGTGGCAGATGATGGATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGCCTCAACTAGGCTGAAAAATGGTGCGATCGGACTGGTCGTACCACAATCgg              <  1:1658797/228‑1 (MQ=255)
                                                                                                                                                                                                                  aGTGGAGGATGACGTATTGGTTACTTACCAGCATCAACCTGTTGGTGTAGCCAAACGGATTGGGTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGCCTCAACTAGGCTGAAAAATGGTGCGATCGGACTGGTCGTACCACAATCgg              <  2:1476444/228‑1 (MQ=255)
                                                                                                                                                                                                                   gTGGGGGATGATGTATGGTTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGCCTCAACTAGGCTGAAAAATGGTGCGATCGGACTGGTCGTACCACAATCgg              <  1:1647861/227‑1 (MQ=255)
                                                                                                                                                                                                                           tgatgTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTAACGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGCCTCAACTAGGCTGAAAAATGGTGCGATCGGACTGGTCGTACCACAATCGGCAGCTAAAt     >  1:1313988/1‑228 (MQ=255)
                                                                                                                                                                                                                              tgTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGCCTCAACTAGCCTGAATACTGGTGTGATCGGACTGGTCGTACCCCAATCGGCAGCTAAATGga  >  2:1987560/1‑228 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                        
GTTCCCGGTGGGCATTGAAGCCCTGATCGGTAAAGTCCGATTTTCTCGGTTGGGGATTAAACTTGCCGAGACGCACAACAAAGGTTATCGCTGGCAGCATGAAGCGGTTATTGCTCTTGCCTCACCCGATAATATGAATGCATTCGAGCTGACACCGCAGGAAGCGGAGGAGTGGTATCGCGGGCGCGATGTTTACCCGCAAGCCGCGCCAGTGGCGGATGATGTATTGGTTACTTTCCAGCATCAACCGATTGGTTTAGCCAAACGGATTGGTTCGCGATTGAAAAACAGCTATCCGCGTGAACTGGTGCGCGATGGGAAACTTTTTACCGGTAACGCCTGACAGCGCACAAAAAAAGCGCACTTTTTGACTGGCACATTCGGCTGCCTCAACTAGGCTGAAAAATGGTGCGATCGGACTGGTCGTACCACAATCGGCAGCTAAATGGA  >  NC_000913/1921320‑1921769

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