Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 558,477 A→C 5.9% E21D (GAA→GAC)  rplR → 50S ribosomal protein L18

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947558,4770AC 5.9% 196.4 / 4.1 68E21D (GAA→GAC) rplR50S ribosomal protein L18
Reads supporting (aligned to +/- strand):  ref base A (29/35);  new base C (2/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

CCAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGC  >  NC_002947/558330‑558620
                                                                                                                                                   |                                                                                                                                               
ccAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTCCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAAc                                                                                                                                                >  8:340634/1‑149 (MQ=255)
     gCCTTACAAAGGTAAAGGTGTGCGTTACGGGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCTAATCACGTCTCAAGATGCACGAACTCGaa                                                                                                                                           <  7:340634/149‑1 (MQ=255)
     gCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGaa                                                                                                                                           <  2:55153/149‑1 (MQ=255)
          aCAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCGCAAGATGCACGAACTCGAAGTCGt                                                                                                                                      >  8:97702/1‑149 (MQ=255)
            aaaggtaaaggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTgc                                                                                                                                    >  8:309303/1‑149 (MQ=255)
                gtaaaggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGACTTCGAGTTCGTgcatct                                                                                                                                >  5:94311/1‑145 (MQ=255)
                 taaaggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGACTTCGAGTTCGTgcatctt                                                                                                                               <  6:94311/149‑6 (MQ=255)
                  aaaggtGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCtgt                                                                                                                              >  6:101426/1‑149 (MQ=255)
                      gtgtGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCgt                                                                                                                          >  5:218096/1‑149 (MQ=255)
                           cGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTc                                                                                                                      <  5:103166/148‑1 (MQ=255)
                             ttACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGc                                                                                                                   >  8:75836/1‑149 (MQ=255)
                                cgcgGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGctc                                                                                                                 <  5:50781/148‑1 (MQ=255)
                                 gcgGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGCCTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACGCGAAGTCGTTCGCCTGTGCGTGGTcccag                                                                                                                  >  6:259883/1‑143 (MQ=255)
                                 gcgGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGctcct                                                                                                               <  4:16089/149‑1 (MQ=255)
                                    gACGAAGTAGTCCGTCGTAAAGAAGCGAAGAAGAAGTAGTGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCgc                                                                                                            <  3:352525/149‑1 (MQ=255)
                                       gAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCgcag                                                                                                          >  3:153870/1‑148 (MQ=255)
                                         agtagtCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGcac                                                                                                       <  5:101426/149‑1 (MQ=255)
                                            agtCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACAt                                                                                                     <  1:271563/148‑1 (MQ=255)
                                             gtCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCt                                                                                                   <  3:327052/149‑1 (MQ=255)
                                                cgtcgtAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTAc                                                                                                 <  6:379831/148‑1 (MQ=255)
                                                     tAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTACCGCTCCTCGCAGCACATCTACGCCCAg                                                                                           <  1:208933/149‑1 (MQ=255)
                                                           agCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCAt                                                                                      <  4:59766/148‑1 (MQ=255)
                                                            gCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCAttt                                                                                    <  2:347993/149‑1 (MQ=255)
                                                            gCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCAtt                                                                                     >  1:347993/1‑148 (MQ=255)
                                                             ccAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTc                                                                                   >  3:220096/1‑149 (MQ=255)
                                                             ccAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAG‑TATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCg                                                                                  >  2:318289/1‑149 (MQ=255)
                                                               aagaagaagTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCgg                                                                                 >  1:277552/1‑149 (MQ=255)
                                                               aagaagaagTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCgg                                                                                 >  5:202908/1‑149 (MQ=255)
                                                               aagaagaagTAGGGCCTAGTAAATGACCCACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCgg                                                                                 >  1:208486/1‑149 (MQ=255)
                                                                agaagaagTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGc                                                                                <  5:259883/149‑1 (MQ=255)
                                                                agaagaagTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGc                                                                                <  7:75836/149‑1 (MQ=255)
                                                                       gTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTgc                                                                                                                                    >  6:306078/1‑90 (MQ=255)
                                                                       gTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTgc                                                                                                                                    <  5:306078/90‑1 (MQ=255)
                                                                       gTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGgc                                                                         <  7:309303/149‑1 (MQ=255)
                                                                        tAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGCGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGgca                                                                        <  5:99511/149‑1 (MQ=255)
                                                                           ggCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTAGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGgcagc                                                                      >  5:16221/1‑148 (MQ=255)
                                                                            gCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGgcagc                                                                      >  2:359536/1‑147 (MQ=255)
                                                                             ccTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAg                                                                   >  3:183249/1‑149 (MQ=255)
                                                                                 gTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGtt                                                                <  6:99289/148‑1 (MQ=255)
                                                                                 gTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTc                                                               >  7:304183/1‑149 (MQ=255)
                                                                                    aaTGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTg                                                             <  4:8631/148‑1 (MQ=255)
                                                                                        aCCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCaa                                                         <  1:149666/148‑1 (MQ=255)
                                                                                          cGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAgcg                                                      >  7:46100/1‑149 (MQ=255)
                                                                                                aaaaGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGCTCTGGCATGCGCCTCGa                                                >  1:157383/1‑149 (MQ=255)
                                                                                                 aaaGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTGCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGAc                                               <  6:357496/149‑1 (MQ=255)
                                                                                                   aGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCGTCGACCt                                             >  5:300918/1‑149 (MQ=255)
                                                                                                     ttATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTg                                           <  4:183249/149‑1 (MQ=255)
                                                                                                     ttATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTAGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTg                                           <  6:16221/149‑1 (MQ=255)
                                                                                                        ttCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGAc                                        <  1:318289/149‑1 (MQ=255)
                                                                                                                  tCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGc                              <  6:218096/149‑1 (MQ=255)
                                                                                                                     cgcTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTg                           <  4:220096/149‑1 (MQ=255)
                                                                                                                      gcTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGa                          <  4:153870/149‑1 (MQ=255)
                                                                                                                       cTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCAGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGAt                         >  5:158812/1‑149 (MQ=255)
                                                                                                                        tCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGtcgat                        >  4:370491/1‑145 (MQ=255)
                                                                                                                          gCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATg                        <  8:46100/147‑1 (MQ=255)
                                                                                                                          gCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGgc                      <  6:202908/149‑1 (MQ=255)
                                                                                                                          gCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGgc                      <  2:277552/149‑1 (MQ=255)
                                                                                                                            aaaGCACGTCTCAAGATGCACGAACTCGAAGTCGTGAGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGgcgc                    >  6:180533/1‑149 (MQ=255)
                                                                                                                             aaGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATgttgc                    >  3:351066/1‑144 (MQ=255)
                                                                                                                               gCACGTATCAAGATGCACGAACTCGAAGTCGTGCGCCTCTGCGTGTTACGCTCCTCGCAGCACATCTACGCACAGGTCATTTCGGAAGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACATTGGACAAAGACCTGCGTGATGGCGCCAc                 >  2:102609/1‑149 (MQ=255)
                                                                                                                                           gcgCACGACTTCGAGTTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGAc                                  >  1:167554/3‑120 (MQ=255)
                                                                                                                                           gcgCACGACTTCGAGTTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGAc                                  <  2:167554/118‑1 (MQ=255)
                                                                                                                                            tGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGAc    <  2:208486/149‑1 (MQ=255)
                                                                                                                                             gCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGAcg   >  7:42808/1‑149 (MQ=255)
                                                                                                                                              cACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGTTGCCAGTGGCGCCATcacgc   <  4:351066/148‑5 (MQ=255)
                                                                                                                                              cACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGAcgc  >  6:106275/1‑149 (MQ=255)
                                                                                                                                              cACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGAcg   <  7:97702/148‑1 (MQ=255)
                                                                                                                                              cACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGGCATTTCGGCCGACGGCAGCAAGGCTCTGGCATGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGAcgc  <  2:157383/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CCAGAGCCTTACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGC  >  NC_002947/558330‑558620

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑