Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,391,247 A→C 5.2% E48A (GAG→GCG) ‡ rlmKL → 23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,391,2470AC 5.2% 521.2 / 13.0 153E48A (GAG→GCG) ‡rlmKL23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (80/65);  new base C (2/6);  total (82/71)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 155 total aligned reads displayed.

CGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGGTGTGCACGCGGTC  >  NC_002947/2391098‑2391395
                                                                                                                                                     |                                                                                                                                                      
cgGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGa                                                                                                                                                        <  2:3252157/150‑1 (MQ=255)
 ggggCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAg                                                                                                                                                       <  1:818178/150‑1 (MQ=255)
      gATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCacac                                                                                                                                                  >  1:1799583/1‑150 (MQ=255)
           gTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCgg                                                                                                                                             <  1:1850119/150‑1 (MQ=255)
            tCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGCCGAGGTGTGCTCGCGCACCTCgtc                                                                                                                                            <  2:3994313/150‑3 (MQ=255)
              gCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGGAGAGGCCAAGGGCCTCGGCCTTGCCGAGGTGTGCTCGCGCACCTCgtcaa                                                                                                                                          <  2:2394562/150‑5 (MQ=255)
              gCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGCCGAGGTGTGCTCGCGCACCTCgtcaa                                                                                                                                          <  1:2732126/150‑5 (MQ=255)
                   gtcagtcaAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTcgc                                                                                                                                   >  1:913786/1‑150 (MQ=255)
                      agtcaAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCgcggc                                                                                                                                <  2:4614198/150‑1 (MQ=255)
                               gCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGAc                                                                                                                       <  2:494527/150‑1 (MQ=255)
                                cccAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACa                                                                                                                      >  2:1508364/1‑150 (MQ=255)
                                 ccAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGTCGAGGTGCGCGAGCACACCTCGGCCATGTCGGCGGcgccgcgaat                                                                                                                       <  1:4150757/150‑5 (MQ=255)
                                  cAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATg                                                                                                                    >  2:5358187/1‑150 (MQ=255)
                                  cAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATg                                                                                                                    >  1:5188795/1‑150 (MQ=255)
                                    aGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGa                                                                                                                  <  2:1692637/150‑1 (MQ=255)
                                       cgcTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAAc                                                                                                               <  1:4067426/150‑1 (MQ=255)
                                        gcTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACt                                                                                                              <  2:3137538/150‑1 (MQ=255)
                                         cTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTg                                                                                                             <  2:3140637/150‑1 (MQ=255)
                                             aaTACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCCCACCTCGGCCCTTCGCGGCGCCGCCGACATGGAAACTGCTTa                                                                                                         >  1:4795219/1‑150 (MQ=255)
                                                  tCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAATAGGACAAGGGACTCGGCCTTGACGAGGTGCGCGAGCACACATCGGCCATTCGCGGCGCCGCCGACATTCAAACAGCTTACCGcc                                                                                                    >  2:2407342/1‑150 (MQ=255)
                                                  tCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGcc                                                                                                    >  2:1236239/1‑150 (MQ=255)
                                                  tCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGcc                                                                                                    >  2:740120/1‑150 (MQ=255)
                                                  tCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGcc                                                                                                    >  2:2438578/1‑150 (MQ=255)
                                                  tCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCACAAAGGCCTTGAAAGCCTGCTCGACGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACAAGGAAACTGCTTACCGcc                                                                                                    >  2:4787146/1‑150 (MQ=255)
                                                       tCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCacac                                                                                                                                                  <  2:91787/101‑1 (MQ=255)
                                                         ggACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCtgtgtgg                                                                                             >  1:4905927/1‑149 (MQ=255)
                                                         ggACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTgtgtgt                                                                                             >  2:3687270/1‑150 (MQ=255)
                                                            ccGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTg                                                                                          <  1:969804/150‑1 (MQ=255)
                                                            ccGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTg                                                                                          >  2:4296958/1‑150 (MQ=255)
                                                               tttCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTc                                                                                       <  2:1850910/150‑1 (MQ=255)
                                                               tttCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTc                                                                                       >  1:1850910/1‑150 (MQ=255)
                                                               tttCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTc                                                                                       >  1:4592070/1‑150 (MQ=255)
                                                               tttCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTc                                                                                       <  2:4592070/150‑1 (MQ=255)
                                                                ttCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCa                                                                                      >  1:286483/1‑150 (MQ=255)
                                                                ttCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCa                                                                                      >  2:3028966/1‑150 (MQ=255)
                                                                 tCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCAc                                                                                     <  2:2116245/150‑1 (MQ=255)
                                                                  cGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACg                                                                                    >  1:2562812/1‑150 (MQ=255)
                                                                   gAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGTGCTCGCGCACCTCgtca                                                                                                                                           <  2:2878386/96‑4 (MQ=38)
                                                                   gAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGt                                                                                   >  1:1056843/1‑150 (MQ=255)
                                                                    aaCTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTc                                                                                  >  1:3884121/1‑150 (MQ=255)
                                                                    aaCTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTc                                                                                  >  1:3884179/1‑150 (MQ=255)
                                                                      cTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAATGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTg                                                                                >  2:2223712/1‑150 (MQ=255)
                                                                       tttACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTgg                                                                               >  2:681248/1‑150 (MQ=255)
                                                                       tttACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTgg                                                                               >  1:4172652/1‑150 (MQ=255)
                                                                       tttACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTgg                                                                               >  1:1590426/1‑150 (MQ=255)
                                                                        ttACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGc                                                                              <  1:195584/150‑1 (MQ=255)
                                                                             tCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAAcc                                                                         >  1:3335415/1‑150 (MQ=255)
                                                                               aCCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACcgc                                                                       >  2:1736772/1‑150 (MQ=255)
                                                                                 cTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCt                                                                                                   >  1:4617691/1‑120 (MQ=255)
                                                                                       cAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTcgcg                                                                                                                                  >  2:242111/1‑83 (MQ=255)
                                                                                       cAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTcgcg                                                                                                                                  <  1:242111/83‑1 (MQ=255)
                                                                                       cAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTcgcg                                                                                                                                  >  1:2411964/1‑83 (MQ=255)
                                                                                       cAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCt                                                               >  2:322038/1‑150 (MQ=255)
                                                                                        aaaGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTg                                                              >  1:1303214/1‑150 (MQ=255)
                                                                                        aaaGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTg                                                              <  1:4296958/150‑1 (MQ=255)
                                                                                        aaaGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTg                                                              <  1:245123/150‑1 (MQ=255)
                                                                                        aaaGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATAGCTg                                                              >  2:2265942/1‑150 (MQ=255)
                                                                                          aGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGt                                                            <  1:217553/150‑1 (MQ=255)
                                                                                            gCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGc                                                          <  2:1863044/150‑1 (MQ=255)
                                                                                              cTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTc                                                        >  2:1468886/1‑150 (MQ=255)
                                                                                              cTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTc                                                        >  2:1468902/1‑150 (MQ=255)
                                                                                              cTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTc                                                        >  2:4665850/1‑150 (MQ=255)
                                                                                              cTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTc                                                        >  1:1204/1‑150 (MQ=255)
                                                                                               ttGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCa                                                       >  2:3561589/1‑150 (MQ=255)
                                                                                               ttGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCa                                                       >  2:3561497/1‑150 (MQ=255)
                                                                                                 gAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTg                                                                                          <  2:5278565/113‑1 (MQ=255)
                                                                                                 gAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTg                                                                                          >  1:5278565/1‑113 (MQ=255)
                                                                                                 gAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAg                                                     >  1:1049849/1‑150 (MQ=255)
                                                                                                 gAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACAGGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAg                                                     >  2:384938/1‑150 (MQ=255)
                                                                                                       cTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTtc                                               >  2:4015287/1‑150 (MQ=255)
                                                                                                        tGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTtct                                              <  2:4816942/150‑1 (MQ=255)
                                                                                                         gCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTtctc                                             >  1:59017/1‑150 (MQ=255)
                                                                                                           tCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGa                                           >  2:104317/1‑150 (MQ=255)
                                                                                                           tCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGa                                           <  2:117904/150‑1 (MQ=255)
                                                                                                           tCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACATCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGa                                           <  2:1116352/150‑1 (MQ=255)
                                                                                                           tCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCAGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGa                                           <  2:2495696/150‑1 (MQ=255)
                                                                                                                 aaGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATgaaga                                     >  1:3192686/1‑150 (MQ=255)
                                                                                                                 aaGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATgaaga                                     >  1:4302809/1‑150 (MQ=255)
                                                                                                                    agGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACg                                  >  1:539632/1‑150 (MQ=255)
                                                                                                                     ggCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGc                                 <  2:5122852/150‑1 (MQ=255)
                                                                                                                      gCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGcc                                >  1:1402934/1‑150 (MQ=255)
                                                                                                                      gCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGcc                                <  2:1517988/150‑1 (MQ=255)
                                                                                                                      gCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGcc                                >  1:353178/1‑150 (MQ=255)
                                                                                                                      gCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGcc                                >  2:214941/1‑150 (MQ=255)
                                                                                                                               cTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATgaa                                       >  2:1443666/1‑134 (MQ=255)
                                                                                                                               cTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCtt                       <  1:1474552/150‑1 (MQ=255)
                                                                                                                                tCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCATACCGCCAGTGTGTCTGGTCACGTCTGGCCAACCCCGTATTGCTGGTGCTCAAGCGCTACTCGATGAAGAACGCCGACGAccaat                      >  2:1154655/1‑147 (MQ=255)
                                                                                                                                 cGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTa                     <  1:1625465/150‑1 (MQ=255)
                                                                                                                                    ccTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTacga                  <  1:1485317/150‑1 (MQ=255)
                                                                                                                                      ttGACGTGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATTAAGAACGCCGACGACCTTTacgacg                <  1:2971042/150‑1 (MQ=255)
                                                                                                                                          cGAGGTGTGCTCGCGCACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTgg                                                             <  1:2462957/101‑1 (MQ=255)
                                                                                                                                          cGAGGTGTGCTCGCGCACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTgg                                                             >  2:2462957/1‑101 (MQ=255)
                                                                                                                                            aGGTGCGCGAGCACACCTCGGCCATTCGCGGTGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGGTGTGc          <  2:2796043/150‑1 (MQ=255)
                                                                                                                                                gcgcGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGGTGTGCAcgc      <  2:286483/150‑1 (MQ=255)
                                                                                                                                                gcgcGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGGTGTGCAcgc      <  1:1212976/150‑1 (MQ=255)
                                                                                                                                                  gcgAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGGTGTGCACGCgg    <  2:1799583/150‑1 (MQ=255)
                                                                                                                                                  gcgAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGATTCTCGATGAAGAACGCCGACGACCTTTACGACGGTGTGCACGCgg    <  1:369996/150‑1 (MQ=255)
                                                                                                                                                   cgAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGGTGTGCACGCGGt   <  1:4850373/150‑1 (MQ=255)
                                                                                                                                                    gAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGGTGTGCACGCGGTc  <  2:4302809/150‑1 (MQ=255)
                                                                                                                                                    gAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGGTGTGCACGCGGTc  >  1:3758855/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGGTGTGCACGCGGTC  >  NC_002947/2391098‑2391395

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