Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 558,484 G→T 5.6% V24F (GTC→TTC)  rplR → 50S ribosomal protein L18

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947558,4840GT 5.6% 221.7 / 3.5 71V24F (GTC→TTC) rplR50S ribosomal protein L18
Reads supporting (aligned to +/- strand):  ref base G (33/34);  new base T (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

ACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGCCACTAAG  >  NC_002947/558340‑558630
                                                                                                                                                |                                                                                                                                                  
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                                                                                           >  5:202908/1‑149 (MQ=255)
                                                     aagaagaagTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCgg                                                                                           >  1:277552/1‑149 (MQ=255)
                                                     aagaagaagTAGGGCCTAGTAAATGACCCACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCgg                                                                                           >  1:208486/1‑149 (MQ=255)
                                                      agaagaagTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGc                                                                                          <  7:75836/149‑1 (MQ=255)
                                                      agaagaagTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGc                                                                                          <  5:259883/149‑1 (MQ=255)
                                                             gTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTgc                                                                                                                                              <  5:306078/90‑1 (MQ=255)
                                                             gTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTgc                                                                                                                                              >  6:306078/1‑90 (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                                <  2:277552/149‑1 (MQ=255)
                                                                                                                gCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGgc                                <  6:202908/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                                            <  2:167554/118‑1 (MQ=255)
                                                                                                                                 gcgCACGACTTCGAGTTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGAc                                            >  1:167554/3‑120 (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)
                                                                                                                                           tCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGCCa       <  6:291141/147‑1 (MQ=255)
                                                                                                                                           tCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGCCa       >  5:291141/1‑147 (MQ=255)
                                                                                                                                           tCGAAGTCGAGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGCCatg     >  4:26728/1‑147 (MQ=255)
                                                                                                                                             gAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGCCCTGCGTGATGGCGCCACTGGCAACATCGACGCGGCCACt     >  7:342540/1‑147 (MQ=255)
                                                                                                                                             gAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGCCACt     >  4:313910/1‑147 (MQ=255)
                                                                                                                                              aaGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGCCACTAAg  >  2:106049/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
ACAAAGGTAAAGGTGTGCGTTACGCGGACGAAGTAGTCCGTCGTAAAGAAGCCAAGAAGAAGTAGGGCCTAGTAAATGACCGACAAAAAAGTTATTCGACTGCGTCGCGCTCGCAAAGCACGTCTCAAGATGCACGAACTCGAAGTCGTGCGCCTGTGCGTGTTCCGCTCCTCGCAGCACATCTACGCCCAGGTCATTTCGGCCGACGGCAGCAAGGTTCTGGCAAGCGCCTCGACCTTGGACAAAGACCTGCGTGATGGCGCCACTGGCAACATCGACGCGGCCACTAAG  >  NC_002947/558340‑558630

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