Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,119,900 G→A 6.7% intergenic (‑140/+39) rlmF ← / ← valS ribosomal RNA large subunit methyltransferase F/valine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,119,9000GA 6.7% 693.1 / 26.3 224intergenic (‑140/+39)rlmF/valSribosomal RNA large subunit methyltransferase F/valine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (121/88);  new base A (11/4);  total (132/92)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 230 total aligned reads displayed.

TCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCACCGCAGGCGGGGCCTTGTCGAC  >  NC_002947/1119751‑1120049
                                                                                                                                                      |                                                                                                                                                     
tCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTcgcg                                                                                                                                                       >  2:3136501/1‑150 (MQ=255)
      tCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGTTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAggg                                                                                                                                                 <  2:2934126/150‑1 (MQ=255)
       cATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGGGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGc                                                                                                                                                <  2:5040147/150‑1 (MQ=255)
          ggCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGt                                                                                                                                             >  1:1558811/1‑150 (MQ=255)
           gCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTa                                                                                                                                            >  1:2498819/1‑150 (MQ=255)
            cAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTAc                                                                                                                                           >  1:3676000/1‑150 (MQ=255)
             aGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACa                                                                                                                                          <  2:6593432/150‑1 (MQ=255)
              gtgtCCAGCCCTCTCCGAGCGGTGGCTTGAAGCCCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACCGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAg                                                                                                                                         <  2:1316843/150‑1 (MQ=255)
               ggtCCAGCCCGCGCCGCGCGGGGGCGTGCCGACCCTGGGGCGGCGGCGCAGCCCTTTTGCGACACCCGGCCGCGCCCCCCGGCGGGTGGCGCCTTGAAGGTTGGCGCCGGCCTGTGGGAGCGGCCTTGGGTCGCGACAGGGCCGTACAGc                                                                                                                                        <  2:6335463/149‑1 (MQ=255)
                   cAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACTGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGcc                                                                                                                                    <  1:4848704/150‑1 (MQ=255)
                           taCGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACCGGAGGGTGGCGCCGGGAAGGTTGGCGCCGTCCGGTGGGCGCGGCGTGGTGTCGCGGAAGGGCCGGACAGCGGCCCCAACGGt                                                                                                                            >  1:1132647/1‑150 (MQ=38)
                                gCGGTGGCTTGATGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACTGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGtctctc                                                                                                                       <  1:4788152/150‑1 (MQ=255)
                                 cGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTctctct                                                                                                                      <  2:4111664/150‑1 (MQ=255)
                                      ctagagaGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGa                                                                                                                 <  2:3595485/144‑1 (MQ=255)
                                       cTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGAt                                                                                                                <  1:3852469/150‑1 (MQ=255)
                                        ttGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGAtt                                                                                                               >  2:1895646/1‑150 (MQ=255)
                                          gAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTAc                                                                                                             >  1:3908753/1‑150 (MQ=255)
                                           ccgACCCTGGGGCTGCGGCGCAGCCCTTTTGCGACACACGGCCGCTCCCACCGGCGGGTGGCGCCTTGACGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCACCGGTCTCTCGAGTGATTACa                                                                                                            <  2:145022/148‑1 (MQ=38)
                                           aaGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACa                                                                                                            <  2:2784284/150‑1 (MQ=255)
                                           aaGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACa                                                                                                            >  2:1488325/1‑150 (MQ=255)
                                           aaGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCGTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACa                                                                                                            >  2:3350997/1‑150 (MQ=255)
                                              aCCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAgcg                                                                                                         >  2:5646935/1‑150 (MQ=255)
                                                ccTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGTAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCt                                                                                                       >  1:4017509/1‑150 (MQ=255)
                                                ccTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACCGCGGCCCCAACGGTCTCTCTCGTGATTACAGCGCt                                                                                                       >  2:3536404/1‑150 (MQ=255)
                                                ccTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCt                                                                                                       >  2:6106403/1‑150 (MQ=255)
                                                ccTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCt                                                                                                       >  1:4484785/1‑150 (MQ=255)
                                                ccTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCt                                                                                                       >  2:2640543/1‑150 (MQ=255)
                                                ccTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCt                                                                                                       >  1:2051581/1‑150 (MQ=255)
                                                aaaagggCTGCGCAGCAGCCCCTTTGCGACACAAGGCCGCTACCACAGGCGGGTGGCGCCGTGCAGGGTGGCGCCGTCCGGGGGGAGCGGCCTTGGGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGAGTACAGCGCt                                                                                                       >  1:5872526/5‑150 (MQ=255)
                                                 cTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTg                                                                                                      <  2:2156922/150‑1 (MQ=255)
                                                       ctgctgCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATcc                                                                                                <  2:5446555/150‑1 (MQ=255)
                                                         gctgcGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCGGGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCGAGTGATTACAGCGCTGCAATCcgc                                                                                              >  1:1328018/1‑150 (MQ=255)
                                                         gcggcgcAGCCCTCTGGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTGGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCTCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCcgc                                                                                              <  1:1200029/147‑1 (MQ=38)
                                                           tgcGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCgcgc                                                                                            <  2:5972337/150‑1 (MQ=255)
                                                                   ccctttcgcgacacaaggccgctcccacaggacggcgccaaccttcaaggcgccaccctccggtgggagcggccgtgtgtcgcAAAAGGGCTGCGCAGCAGCCCCAGGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATgctgct                                                                                    >  2:328227‑M2/84‑150 (MQ=255)
                                                                     ttttGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCg                                                                                  >  1:311454/1‑150 (MQ=255)
                                                                     ttttGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCGTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCAGGCTGCTCg                                                                                  >  1:2732605/1‑150 (MQ=255)
                                                                      tttGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGTTCgg                                                                                 >  2:3208675/1‑150 (MQ=255)
                                                                         GCGACACAAGGCCGCTCCCACAGGACGGCGCCAACCTTCAAGGCGCCACCCTCCTGTGGGAGCGGCCTTGTGTCGCAAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGa                                                                              <  2:1733808‑M1/150‑1 (MQ=255)
                                                                           GACACAAGGCCGCTCCCACAGGATGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAg                                                                            <  1:5013679‑M1/150‑1 (MQ=255)
                                                                           GACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAg                                                                            <  2:1671898‑M1/150‑1 (MQ=255)
                                                                            ACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGt                                                                           <  2:1867907‑M1/150‑1 (MQ=255)
                                                                               CAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTgg                                                                        <  1:2769022‑M1/150‑1 (MQ=255)
                                                                                AAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAgcgc                                                                                                        <  2:205318‑M1/117‑1 (MQ=255)
                                                                                AAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGc                                                                       <  1:1426400‑M1/150‑1 (MQ=255)
                                                                                AAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGc                                                                       <  2:4824568‑M1/150‑1 (MQ=255)
                                                                                AAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGc                                                                       <  1:5778162‑M1/150‑1 (MQ=255)
                                                                                AAGGCCGCTCCCACAGGAGGGTGGCACCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGc                                                                       <  1:6339744‑M1/150‑1 (MQ=255)
                                                                                     CGCTCCCACAGGCGGGTGGCGCCTTGATGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAggg                                                                  >  1:5321652‑M1/1‑150 (MQ=255)
                                                                                     CGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAggg                                                                  >  2:3707899‑M1/1‑150 (MQ=255)
                                                                                        TCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGccc                                                                                                                                   <  1:721130‑M1/82‑1 (MQ=255)
                                                                                        TCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCt                                                               <  2:6794005‑M1/150‑1 (MQ=255)
                                                                                        TCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCt                                                               <  2:1001891‑M1/150‑1 (MQ=255)
                                                                                        TCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCt                                                               <  1:2626639‑M1/150‑1 (MQ=255)
                                                                                         tcccacaggacggcgccaaccttcaaggcgccaccctcctgtgggagcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCt                                                               >  1:169726‑M2/62‑150 (MQ=255)
                                                                                         tcccacaggacggcgccaaccttcaaggcgccaccctcctgtgggagcggccttgtgtcgcACAAGGGCTGCGCAGCAGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCt                                                               >  2:459487‑M2/62‑150 (MQ=255)
                                                                                         tcccacaggacggcgccaaccttcaaggcgccaccctcctgtgggagcggccttgtgtcgcAAAATGGCTGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCt                                                               >  1:6284273‑M2/62‑150 (MQ=255)
                                                                                         tcccacaggacggcgccaaccttcaaggcgccaccctccggtgggagcggccttgtgtcgcGAAAGGGCCGTACCGCGGCCCCAACGGTCTCTCGAGTGAGTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGACGTTGGCCAGGGCCt                                                               >  2:5579730‑M2/62‑150 (MQ=255)
                                                                                          cccacaggacggcgccaaccttcaaggcgccaccctcctgtgggagcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTg                                                              >  1:3134239‑M2/61‑150 (MQ=255)
                                                                                          CCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGt                                                             <  1:1900852‑M1/150‑1 (MQ=255)
                                                                                           ccacaggacggcgccaaccttcaaggcgccaccctcctgtgggagcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGt                                                             >  2:2343644‑M2/60‑150 (MQ=255)
                                                                                           ccacaggacggcgccaaccttcaaggcgccaccctcctgtgggagcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGt                                                             >  1:5126750‑M2/60‑150 (MQ=255)
                                                                                            cacaggacggcgccaaccttcaaggcgccaccctcctgtgggagcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGtt                                                            >  2:797795‑M2/59‑150 (MQ=255)
                                                                                              AGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCgg                                                         <  1:4192645‑M1/150‑1 (MQ=255)
                                                                                               GGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGa                                                        <  1:2652535‑M1/150‑1 (MQ=255)
                                                                                                GAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGAc                                                       <  2:2110354‑M1/150‑1 (MQ=255)
                                                                                                GAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGAc                                                       <  1:2985313‑M1/150‑1 (MQ=255)
                                                                                                 gacggcgccaaccttcaaggcgccaccctcctgtgggagcggccttgtgtcgcAAAAGGGCTGCGCAGCAGCCCCAGGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACt                                                      >  1:654234‑M2/54‑150 (MQ=255)
                                                                                                   ggTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCg                                                    <  2:4664674/150‑1 (MQ=255)
                                                                                                       ggcgccaaccttcaaggcgccgtcctgtgggagcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGcc                                                 >  1:157937‑M2/48‑150 (MQ=255)
                                                                                                       gcgcCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCa                                                >  1:6566803/1‑150 (MQ=255)
                                                                                                        cgcCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAg                                               <  2:4914225/150‑1 (MQ=255)
                                                                                                            cTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCGTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCtt                                            >  2:1414121/1‑150 (MQ=255)
                                                                                                              ttcaaggcgccaccctcctgtgggagcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTgg                                          >  2:4497191‑M2/41‑150 (MQ=255)
                                                                                                              tGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTgg                                          >  1:3155605/1‑150 (MQ=255)
                                                                                                                 aaggcgccaccctcctgtgggagcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGgcgc                                       >  2:1884898‑M2/38‑150 (MQ=255)
                                                                                                                     tGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTc                                   >  1:759079/1‑150 (MQ=255)
                                                                                                                     tGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTc                                   >  2:5355798/1‑150 (MQ=255)
                                                                                                                          ccGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCcgc                                                                                              <  1:5455171/86‑1 (MQ=255)
                                                                                                                          ccGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCcgc                                                                                              >  2:5455171/1‑86 (MQ=255)
                                                                                                                           cgacctgggggtgcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTtctc                             <  1:3475198‑M2/123‑1 (MQ=255)
                                                                                                                           cctcctgtgggagcggccttgtgtcgcAAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTtctc                             <  2:669201‑M2/123‑1 (MQ=255)
                                                                                                                           cGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTtctc                             >  2:4766914/1‑150 (MQ=255)
                                                                                                                             tcctgtgggagcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGa                           >  1:5556227‑M2/26‑150 (MQ=255)
                                                                                                                                tgtgggagcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTtctc                             >  1:4596035‑M2/23‑145 (MQ=255)
                                                                                                                                tgtgggagcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCa                        <  2:5126750‑M2/128‑1 (MQ=255)
                                                                                                                                tgtgggagcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCa                        >  2:6744346‑M2/23‑150 (MQ=255)
                                                                                                                                  tgggagcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCAcc                      >  1:6039014‑M2/21‑150 (MQ=255)
                                                                                                                                  tgggagcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCAcc                      >  1:2572812‑M2/21‑150 (MQ=255)
                                                                                                                                    ggagcggccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCACCGc                    >  1:586059‑M2/19‑150 (MQ=255)
                                                                                                                                          gccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCACCGCAGGCgg              <  1:2426667‑M2/138‑1 (MQ=255)
                                                                                                                                          gccttgtgtcgcAAAAGGGCTGCGCAGCAGCCCCAGGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCACCGCAGGCggg             <  1:5286180‑M2/138‑1 (MQ=255)
                                                                                                                                           ccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCACCGCAGGCggg             <  1:1895646‑M2/139‑1 (MQ=255)
                                                                                                                                           ccttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCACCGCAGGCggg             >  2:37094‑M2/12‑150 (MQ=255)
                                                                                                                                           ccttgtggcgcGAAAGGGCCGTACAGCGTCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCCTGGCGCGTTCCTGCTCGATCCCCGCAGGCggg             <  2:5206022‑M2/139‑1 (MQ=255)
                                                                                                                                             ttgtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCACCGCAGGCGGGGc           >  1:915701‑M2/10‑150 (MQ=255)
                                                                                                                                             ttgtgtcgcGAAAGGGCCGTACAGCGGCCACAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCACCGCAGGCGGGGc           <  2:4529780‑M2/141‑1 (MQ=255)
                                                                                                                                               gtgtcgcGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCACCGCAGGCGGGGCCt         >  2:4240808‑M2/8‑150 (MQ=255)
                                                                                                                                                     cGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCACCGCAGGCGGGGCCTTGTCGa   >  1:4910159‑M2/2‑150 (MQ=255)
                                                                                                                                                      gAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCACCGCAGGCGGGGCCTTGTcgac  <  1:5355798/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TCGGGGTCATGGCAGTGTCCAGCCCTATACGAGCGGTGGCTTGAAGACCCTGGGGCTGCTGCGCAGCCCTTTTGCGACACAAGGCCGCTCCCACAGGAGGGTGGCGCCTTGAAGGTTGGCGCCGTCCTGTGGGAGCGGCCTTGTGTCGCGAAAGGGCCGTACAGCGGCCCCAACGGTCTCTCTAGTGATTACAGCGCTGCAATCCGCGCATGCTGCTCGGTGAAGTTGGCCAGGGCCTGTTCGGACTCGGCCAGCTTGGCGCGTTCCTTCTCGATCACCGCAGGCGGGGCCTTGTCGAC  >  NC_002947/1119751‑1120049

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