Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,893,530 T→C 6.6% F284L (TTC→CTC)  rarD → chloramphenicol‑sensitive protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,893,5300TC 6.6% 311.7 / 8.2 105F284L (TTC→CTC) rarDchloramphenicol‑sensitive protein
Reads supporting (aligned to +/- strand):  ref base T (45/53);  new base C (2/5);  total (47/58)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 110 total aligned reads displayed.

TGCCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAGCTCATCACC  >  NC_002947/3893381‑3893678
                                                                                                                                                     |                                                                                                                                                    
tgcCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGTGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCt                                                                                                                                                      >  1:299986/1‑150 (MQ=255)
tgcCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCt                                                                                                                                                      >  1:1755543/1‑150 (MQ=255)
tgcCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCt                                                                                                                                                      >  2:652371/1‑150 (MQ=255)
    aTTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCaa                                                                                                                                                  <  1:30931/150‑1 (MQ=255)
       cAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCTTTGCTGCTGGGAGAGACCGGCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAAGCCCTCAAGCa                                                                                                                                               <  2:565570/150‑1 (MQ=255)
        aGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAc                                                                                                                                              >  2:1524533/1‑150 (MQ=255)
        aGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAc                                                                                                                                              <  2:2346229/150‑1 (MQ=255)
        aGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAc                                                                                                                                              <  2:1688890/150‑1 (MQ=255)
        aGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAc                                                                                                                                              <  1:1144309/150‑1 (MQ=255)
        aGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAAGCCCTCAAGCAc                                                                                                                                              <  2:1449901/150‑1 (MQ=255)
        aGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAAGCCCTCAAGCAc                                                                                                                                              <  2:1449431/150‑1 (MQ=255)
         gCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCAcc                                                                                                                                             >  2:305540/1‑150 (MQ=255)
         gCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGAACTGGTGCTGGTGCTTGAGGGCTTCAAGCAcc                                                                                                                                             <  2:1347828/150‑1 (MQ=255)
           cTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTg                                                                                                                                           >  2:106623/1‑150 (MQ=255)
                   ccTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCa                                                                                                                                   <  2:741431/150‑1 (MQ=255)
                    cTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAg                                                                                                                                  >  2:1490710/1‑150 (MQ=255)
                       cTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGt                                                                                                                               >  1:1350357/1‑150 (MQ=255)
                        tCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTc                                                                                                                              >  2:681511/1‑150 (MQ=255)
                               cGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTc                                                                                                                       <  2:2360704/150‑1 (MQ=255)
                               cGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGACGCAGTACTGGTGCTGGTTCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTc                                                                                                                       <  2:1902359/150‑1 (MQ=255)
                                   gAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCgggg                                                                                                                   >  1:1808155/1‑148 (MQ=255)
                                    attgcGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGgtgt                                                                                                                  >  2:576876/5‑150 (MQ=255)
                                     aCCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTa                                                                                                                 >  2:1893041/1‑150 (MQ=255)
                                        cGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGc                                                                                                              >  2:887637/1‑150 (MQ=255)
                                        cGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGc                                                                                                              <  2:1755543/150‑1 (MQ=255)
                                          tgctgcTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCg                                                                                                            <  2:195359/150‑1 (MQ=255)
                                          tgctgcTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAAGCCCTCAAGCAc                                                                                                                                              >  1:85331/1‑116 (MQ=255)
                                          tgctgcTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAAGCCCTCAAGCAc                                                                                                                                              <  2:85331/116‑1 (MQ=255)
                                             tgctgGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAAc                                                                                                         >  2:28119/1‑150 (MQ=255)
                                             tgctgGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGATCGGTGTAAGCCGAAc                                                                                                         >  2:1223996/1‑150 (MQ=255)
                                               ctgGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCt                                                                                                       <  2:321572/150‑1 (MQ=255)
                                                tgGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTg                                                                                                      <  1:305540/150‑1 (MQ=255)
                                                  gTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGc                                                                                                    <  2:2204038/150‑1 (MQ=255)
                                                       cgtggcgtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGAAATACCTTGTTGCCAATCTGGGCCGCAGCAGATTTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTAcc                                                                                               <  1:1116037/150‑1 (MQ=255)
                                                        gtggcgtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTgg                                                                                                     >  1:31059/1‑143 (MQ=255)
                                                         tggcgtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTc                                                                                                                                                    >  1:464486/1‑95 (MQ=255)
                                                         tggcgtTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTc                                                                                                                                                    <  2:464486/95‑1 (MQ=255)
                                                                gctgctGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCAGGTGCTTGAAGCCCTCAAGCAccagc                                                                                                                                          >  2:171335/1‑95 (MQ=37)
                                                                gctgctGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCAGGTGCTTGAAGCCCTCAAGCACCAGc                                                                                                                                          <  1:171335/98‑1 (MQ=37)
                                                                  tgctgGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTc                                                                                    <  1:1893041/150‑1 (MQ=255)
                                                                        gCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGc                                                                              >  2:466037/1‑150 (MQ=255)
                                                                        gCGAAACCATCGGCCCAGATCAGTGGCTGACCTGCCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCTCCCTTCGGTGTAAGCCGAACCTGGCGTACCCCCCGCTCTTCAACCGc                                                                              >  2:489898/1‑150 (MQ=255)
                                                                         cGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGCAGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGcc                                                                             <  1:158136/150‑1 (MQ=255)
                                                                          gAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCt                                                                                                                                         >  2:184154/1‑89 (MQ=255)
                                                                          gAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCt                                                                            >  1:2449204/1‑150 (MQ=255)
                                                                               cATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGgtgt                                                                                                                  <  1:1580019/107‑1 (MQ=255)
                                                                               cATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGgtgt                                                                                                                  >  2:1580019/1‑107 (MQ=255)
                                                                                aTCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACgg                                                                      >  2:2350998/1‑150 (MQ=255)
                                                                                aTCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTCGCGTACCCGTCGCTCTTCAACCGCCTCCACgg                                                                      >  2:456846/1‑150 (MQ=255)
                                                                                  cGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTACGCCGAACCTGGCGTACCCGTCGCCATTCAACCGCCTCCACGGTc                                                                    <  2:1166428/150‑1 (MQ=255)
                                                                                    gCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAg                                                                  <  1:1898725/150‑1 (MQ=255)
                                                                                      ccGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCa                                                                >  2:185291/1‑150 (MQ=255)
                                                                                      ccGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCa                                                                >  2:2419673/1‑150 (MQ=255)
                                                                                         gATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATcc                                                             >  1:426044/1‑150 (MQ=255)
                                                                                          aTCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGcc                                                                                                             >  1:1612499/1‑101 (MQ=255)
                                                                                            cAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGc                                                          <  2:1433284/150‑1 (MQ=255)
                                                                                            cAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGc                                                          <  1:1524533/150‑1 (MQ=255)
                                                                                               tGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCAt                                                       <  2:299986/150‑1 (MQ=255)
                                                                                               tGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCCGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCAt                                                       <  2:1643199/150‑1 (MQ=255)
                                                                                                   tGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCa                                                   >  2:181498/1‑150 (MQ=255)
                                                                                                    gACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAg                                                  >  1:1847234/1‑150 (MQ=255)
                                                                                                    gACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAg                                                  <  1:2008594/150‑1 (MQ=255)
                                                                                                       ctacctGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTAcc                                                                                               <  1:550840/102‑1 (MQ=255)
                                                                                                       ctacctGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAg                                                  <  1:955144/147‑1 (MQ=255)
                                                                                                       ctacctGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCt                                               >  2:2391482/1‑150 (MQ=255)
                                                                                                       ctacctGCCAATCTGGGCCGCAGTACTGGTGCTGATGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCCACCGCCTCCACGGTCAGCATCCCGCCATTCCAGGCt                                               >  2:698595/1‑150 (MQ=255)
                                                                                                        tacctGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAGGc                                                                                                              >  2:1364163/1‑86 (MQ=255)
                                                                                                           ctGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGcc                                           >  2:1994430/1‑150 (MQ=255)
                                                                                                            tGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCgg                                                                                                                     >  2:1400958/1‑75 (MQ=255)
                                                                                                            tGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCgg                                                                                                                     <  1:1400958/75‑1 (MQ=255)
                                                                                                            tGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCg                                          <  1:1041304/150‑1 (MQ=255)
                                                                                                            tGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCg                                          <  2:660455/150‑1 (MQ=255)
                                                                                                             gCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCGCGCTATCCgg                                         >  1:701378/1‑150 (MQ=255)
                                                                                                             gCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCgg                                         >  1:1007795/1‑150 (MQ=255)
                                                                                                               cAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTc                                       <  1:132335/150‑1 (MQ=255)
                                                                                                               cAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTc                                       <  1:197883/150‑1 (MQ=255)
                                                                                                                aaTCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCa                                      >  2:101209/1‑150 (MQ=255)
                                                                                                                 aTCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCAc                                     <  1:2350998/150‑1 (MQ=255)
                                                                                                                    tGGGCCGCAGTACTGGTGCTGGTGCTTGAGGTCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGa                                  <  2:388088/150‑1 (MQ=255)
                                                                                                                       gCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCa                                                                   >  1:841470/1‑114 (MQ=255)
                                                                                                                       gCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCa                                                                   <  2:841470/114‑1 (MQ=255)
                                                                                                                        ccGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCa                                                            <  2:1416071/120‑1 (MQ=255)
                                                                                                                        ccGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCa                                                            >  1:1416071/1‑120 (MQ=255)
                                                                                                                         cGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATc                             >  2:325249/1‑150 (MQ=255)
                                                                                                                              tACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAg                        <  2:410516/150‑1 (MQ=255)
                                                                                                                                 tggtgctggtgcTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCg                     <  1:185291/150‑1 (MQ=255)
                                                                                                                                  ggtgctggtgctTGAGGGCTTCAAGCGCCGGCTACCTCAGCGTCGAAGTTCGGTGTAAGAGCATCCAACACCTGCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCgg                    <  1:1223996/150‑1 (MQ=255)
                                                                                                                                        ggtgctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGc              <  1:28119/150‑1 (MQ=255)
                                                                                                                                        ggtgctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGc              <  1:1218726/150‑1 (MQ=255)
                                                                                                                                          tgctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCaa            <  2:426044/150‑1 (MQ=255)
                                                                                                                                           gctTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAg           >  1:715792/1‑150 (MQ=255)
                                                                                                                                             ttGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAGCt         >  1:2044856/1‑150 (MQ=255)
                                                                                                                                                aGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAGCtcat      <  1:2295008/150‑1 (MQ=255)
                                                                                                                                                 gggCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAGCtcatc     >  1:1078327/1‑150 (MQ=255)
                                                                                                                                                  ggCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTc                                                     >  2:418989/1‑101 (MQ=255)
                                                                                                                                                  ggCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTc                                                     <  1:418989/101‑1 (MQ=255)
                                                                                                                                                    cTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCa                                                            <  1:1511821/92‑1 (MQ=255)
                                                                                                                                                    cTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCa                                                            >  2:1511821/1‑92 (MQ=255)
                                                                                                                                                    cTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAACCGGTCACCGACCATCGGCAGGCGGTCCAGCAAGCTCATCAcc  <  1:1212296/150‑1 (MQ=255)
                                                                                                                                                     ttCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAgg                                                 >  1:198127/1‑102 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCCATTCAGCCTGTTCGGCCTGCTCAGCTACGTCGAACCCGTGCTGCTGGTAGGCGTGGCGTTGCTGCTGGGCGAAACCATCGGCCCGGATCAGTGGCTGACCTACCTGCCAATCTGGGCCGCAGTACTGGTGCTGGTGCTTGAGGGCTTCAAGCACCTGCTACGTCAGCGTCGCCGTTCGGTGTAAGCCGAACCTGGCGTACCCGTCGCTCTTCAACCGCCTCCACGGTCAGCATCCAGCCATTCCAGGCTAGCCGGTCACCGACCATCGGCAGGCGGTCCAGCAAGCTCATCACC  >  NC_002947/3893381‑3893678

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