Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,673,071 Δ1 bp 9.2% intergenic (+97/‑70) prpF → / → prpD aconitate isomerase/2‑methylcitrate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,673,0710A. 9.2% 533.9 / 50.4 173intergenic (+97/‑70)prpF/prpDaconitate isomerase/2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (92/63);  new base . (6/10);  total (98/75)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-01

Only 100 of 179 total aligned reads displayed.

CTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGC  >  NC_002947/2672925‑2673218
                                                                                                                                                    |                                                                                                                                                   
cTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATcc                                                                                                                                                  <  1:2486685/150‑1 (MQ=255)
 tCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATccc                                                                                                                                                 <  1:6339694/150‑1 (MQ=255)
   gcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtg                                                                                                                                               <  2:8668073/150‑1 (MQ=255)
    cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                                              <  2:6337549/150‑1 (MQ=255)
     gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGCATTCAAACCCTGTGATTACCTGTgg                                                                                                                                            <  1:4916224/150‑1 (MQ=25)
      tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTgg                                                                                                                                            <  2:7268204/150‑1 (MQ=255)
      tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTgg                                                                                                                                            <  1:6918780/150‑1 (MQ=255)
      tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTgg                                                                                                                                            <  1:6551344/150‑1 (MQ=18)
        cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                                          <  2:8031394/150‑1 (MQ=18)
         tGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGGTTGAATGCTGTGATATCCCGGTGGGAg                                                                                                                                         >  2:1478553/1‑150 (MQ=255)
         tGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGAATCACAGGGGTTGAATGCGGTGATATCCCTGTGGGAg                                                                                                                                         >  1:7240420/1‑150 (MQ=255)
          gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                                                      <  1:1777545/150‑1 (MQ=21)
          gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGc                                                                                                                                        <  2:4811268/150‑1 (MQ=18)
          gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCGTGCCCCCGCCCCCCCGGAAAAAAAGGGGTTGTATTGTGTGATATCCCGGTGGGGGc                                                                                                                                        >  1:5942817/1‑150 (MQ=255)
            tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                                                      <  2:1350570/150‑1 (MQ=18)
            tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                                                      <  2:7818233/150‑1 (MQ=18)
             ggAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTCAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCggg                                                                                                                                     <  1:6706792/150‑1 (MQ=255)
             ggAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCggg                                                                                                                                     <  2:7686608/150‑1 (MQ=255)
              gAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGc                                                                                                                                    <  2:1980104/150‑1 (MQ=255)
                   cTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGcc                                                                                                                               >  2:2437845/1‑150 (MQ=255)
                       gTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATTCTGTGATATACCGGTGGGAGCGGGCAAGCCcgcg                                                                                                                           >  2:1061464/1‑150 (MQ=255)
                       gTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGTGCGGGCAAGCCcgcg                                                                                                                           >  2:8386696/1‑150 (MQ=255)
                           gcgTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCCCTACCACAGGTAAAAACAGGGTTTGAAAGCGGTGATAAACCTGTGGGAGCGGGCAACCCCGCGAAgg                                                                                                                       >  2:3945009/1‑150 (MQ=255)
                             gTGCCTGGCGATCGCTACTAAAGCACCACGGGCGCTCGGCTCCAAGCGGGGCCGCTTTGCGGCCCTTCGCGGGCTGGCCCGCTCCCCCAGGTAATCCCTGGGTTTGACTGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGc                                                                                                                     <  2:3560846/150‑1 (MQ=37)
                              gcccgggCGATCGCTAATAAAGCACCAGTGGCGGTCGGCTCCACGTGGGGCGGCTTTGCGGCCCTTCGCGGGCTTGCCCGCCCCCACCGGTAACCACCGGGTTTGAATGCGGTGATAGCCGTGTGGGCGCGGGCACGCCCGCGAAGGGCt                                                                                                                    <  2:4573530/145‑1 (MQ=255)
                               gCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGAAAGCCCGCGAAGGGCTg                                                                                                                   <  2:2040176/150‑1 (MQ=255)
                                ccTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGc                                                                                                                  >  2:8250300/1‑150 (MQ=255)
                                     cGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCACCCACAGGTAATCCCAGGGTTTGAATGCTGTGATATCCCGGTGGGAGCGGGCAAGCCCGCGAAGGGCGGCAcagc                                                                                                             >  2:2198684/1‑150 (MQ=255)
                                         aGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCACAGGTAATCACAGGGGTTGTATGCTGTGATATCCCTGGGGGAGCGGGCAAGCCCGCGAAGGGCTGCACCGCAGcc                                                                                                         >  2:3370791/1‑150 (MQ=38)
                                           cTTCTAACGCACCACTGGCGCTCGGGTTCACGCGGGGCCGCTGTGCGGCCCTGCGCGGGCTTGCCCGCTCCCACCGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGcccc                                                                                                       <  2:6953653/150‑1 (MQ=25)
                                             tCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTGGCCCGCTCCCACAGGTAATCACAGGGTTGGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCgg                                                                                                     <  2:539944/150‑1 (MQ=255)
                                               tAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCCCAGGTAATCACAGGGGTTGAATGCTGTGATATCCCGGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCa                                                                                                   >  2:1582688/1‑150 (MQ=255)
                                                   gCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTc                                                                                               >  2:2700757/1‑150 (MQ=255)
                                                       cccTGGCGCTGGGATGCAAGGGGGGCGGATTTGCGGCCATTCGCGGGCTTGCCCGCTCCCACAGTTAATCACTGGTTTTGATTGGTGTGATATCCGTTTGGGAGGGTCCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                                           <  1:8076591/148‑1 (MQ=255)
                                                       cacTGGCGCTCGGCTTCATGTGGGGCGGCTTGGCGGCCCTTCGGGGGCTTGCCGGCTCCCACGGGTAAACACCGGTTTTGAATGCTGTGATATCCCTGTGGGAGCGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                                           <  2:2348543/150‑1 (MQ=255)
                                                       cacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGTATTGTGTGATATACCGGTGGGTGCGGGCAAGCCCGCGAATGGCTGCACAGCAGCCCCGGCAAttattaa                                                                                           >  2:1907279/1‑145 (MQ=37)
                                                       cacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                                      <  1:5813098/107‑1 (MQ=255)
                                                       cacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                                      >  2:5813098/1‑107 (MQ=255)
                                                                 cGGCTGCAAGAGGGGCCGGTTTGCGGCCCTTGGCGGGCTTGCCCGCTCCCACAGTTAATCACAGGGTTTGAATGCTGTTATATCCCTGTGGGAGCGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGa                                                                                 <  2:6614064/150‑1 (MQ=38)
                                                                      tcatgcggggCCGCTTTGCGGCCCTGCGGGGGCTTGCCCGCTCCCCCTGTTAATCCCTGGGTTTGAATGGTGTGATTTCCCTGTGGGCGCGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGa                                                                            <  1:3617110/144‑1 (MQ=255)
                                                                        atgtggggCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGGTATCACAGGGTTTTAAGGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAgc                                                                          <  2:7761105/146‑1 (MQ=25)
                                                                               gCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACg                                                                   >  2:5119573/1‑150 (MQ=255)
                                                                                  gCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGGGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGCGCTGGATATGAGCGCCAACGTAg                                                                >  2:7880655/1‑150 (MQ=255)
                                                                                  gCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCGGGGAGATCCCGGGGGGAGCGGGCCAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCCCTTCTGACGGGAGCTGGATATGAGCGCCAACGtcg                                                                >  2:7880397/1‑148 (MQ=37)
                                                                                    tttGCGGCCCTTCGGGGGCTGGCCCGCTCCCACCGCTACACACCGGTTTTGAATGGGGTTATACCCCGGTGGGCGGGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGAc                                                              <  1:3370791/150‑1 (MQ=255)
                                                                                    tttGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGTTAATCCCTGGTTTTTATTGCTGTGATATCCCTGTGGGAGGGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGAc                                                              <  1:3503866/150‑1 (MQ=38)
                                                                                          gCCCTTCGGGGGCTTGCCCGCTCCCACAGGTAAACACTGGTTTTGAATGCTGTGATAGCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAAc                                                        <  2:113603/150‑1 (MQ=38)
                                                                                             cTTCGCGGGCTTGCCCGCGCCCACAGGGCTATCACAGCATTCAAACCCGGTGATTACCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCCACGTAGACCTCAACGAt                                                     >  2:506483/1‑150 (MQ=14)
                                                                                                     gcgggCCCGCGCCCCCAGTTATTCACCGGGTTTGAAGGCTGTGCTAGCCCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGa                                             <  2:3654772/146‑1 (MQ=255)
                                                                                                       ttGCCCGCTCCCACAGGTAATCACCGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt                                           >  1:2125943/1‑150 (MQ=255)
                                                                                                       ttGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt                                           >  2:2882792/1‑150 (MQ=255)
                                                                                                       ggccccGCCCCCCCCGGAAACCACTGGTTTTGTTTGGGGTGATAGCCGTGTGGGCGGGGGCAACCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt                                           <  1:5221786/147‑1 (MQ=255)
                                                                                                        tGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                                          >  2:12231/1‑150 (MQ=255)
                                                                                                        ggCCCGCGCCCCCAGGTAAACACCGGTTTTGAATGCTGGGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                                          <  2:6045785/149‑1 (MQ=39)
                                                                                                         gCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc                                         <  2:4286743/150‑1 (MQ=255)
                                                                                                         gCCCGCTCCCACAGGTAAACACTGGGTTTGATTGCTTTGATAGCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc                                         <  2:449336/150‑1 (MQ=39)
                                                                                                           ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa                                       >  1:332676/1‑150 (MQ=255)
                                                                                                           ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa                                       >  2:8211635/1‑150 (MQ=255)
                                                                                                           ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGACCGGATTACGa                                       >  2:2238621/1‑150 (MQ=255)
                                                                                                           ccGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACgc                                       >  2:6107030/1‑149 (MQ=18)
                                                                                                             gCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGAcc                                     >  1:5132962/1‑150 (MQ=255)
                                                                                                             gCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGAcc                                     >  1:2621857/1‑150 (MQ=255)
                                                                                                               tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                                   <  2:2105376/150‑1 (MQ=255)
                                                                                                                 ccACAGGTACTCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCGGAAGGGAGCTGGATATGAGCGCCAACGTCGACCTCAACGATCGCCCGGATTACGACCGGGt                                 >  2:304781/1‑150 (MQ=255)
                                                                                                                 ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATCTGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                                 >  2:4639919/1‑150 (MQ=255)
                                                                                                                 ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                                 >  1:4497434/1‑150 (MQ=255)
                                                                                                                 ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                                 >  1:2054438/1‑150 (MQ=255)
                                                                                                                 ccACAGGTAATCACAGGGTTTGAATGCTATGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                                 >  1:7692562/1‑150 (MQ=255)
                                                                                                                  cacaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTg                                >  1:3192688/1‑150 (MQ=255)
                                                                                                                  cacaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTg                                >  2:6976087/1‑150 (MQ=255)
                                                                                                                  cacaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTg                                >  1:8742187/1‑150 (MQ=255)
                                                                                                                        gTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACa                           >  1:5586221/1‑150 (MQ=255)
                                                                                                                              cacaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAccc                      >  2:5411490/1‑150 (MQ=255)
                                                                                                                              cacaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAccc                      >  1:6417400/1‑150 (MQ=255)
                                                                                                                               acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt                     >  2:2400045/1‑150 (MQ=255)
                                                                                                                               acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt                     >  2:476191/1‑150 (MQ=255)
                                                                                                                                ctggttttGAATGCTGTGTTATCCCTGTGGGTGGGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg                    <  2:5942817/145‑1 (MQ=255)
                                                                                                                                ctgggTTTGAATGCTGTTATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg                    <  2:1814862/148‑1 (MQ=255)
                                                                                                                                ctgggTTGGAAGGCGGGGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTCCAAACCCTg                    <  2:3718132/148‑1 (MQ=39)
                                                                                                                                 aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg                   >  2:1561116/1‑150 (MQ=255)
                                                                                                                                 aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg                   >  2:4120508/1‑150 (MQ=255)
                                                                                                                                 aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGACTACGACCGGGTGCTACAAACCCTgg                   >  1:711827/1‑150 (MQ=255)
                                                                                                                                   ggTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGAAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGcc                 >  2:5199971/1‑150 (MQ=255)
                                                                                                                                     tttGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGa               >  1:8283454/1‑150 (MQ=255)
                                                                                                                                     tttGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGa               <  2:3852009/150‑1 (MQ=255)
                                                                                                                                      ttGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGAt              >  2:8207271/1‑150 (MQ=255)
                                                                                                                                        aaacccTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTAc           >  1:1274694/6‑150 (MQ=255)
                                                                                                                                          aggCTGTGATATCCCTGTGGGAGCGGGCACGCCCGCGACGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACg          <  2:5593810/148‑1 (MQ=255)
                                                                                                                                           tGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGc         >  1:2920164/1‑150 (MQ=255)
                                                                                                                                            ccTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGccc       >  2:320806/2‑150 (MQ=255)
                                                                                                                                             cTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGccc       >  1:1326323/1‑150 (MQ=255)
                                                                                                                                              tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt      >  1:4062531/1‑150 (MQ=255)
                                                                                                                                              tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt      >  1:6658456/1‑150 (MQ=255)
                                                                                                                                              tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt      >  2:6369635/1‑150 (MQ=255)
                                                                                                                                                ttatatCCGTGTGGGAGCGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCg    <  1:512316/148‑1 (MQ=255)
                                                                                                                                                tgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCg    >  1:147168/1‑150 (MQ=255)
                                                                                                                                                 gATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCgg   >  1:4618721/1‑150 (MQ=255)
                                                                                                                                                 gATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCgg   >  1:3853957/1‑150 (MQ=255)
                                                                                                                                                  atatCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGc  >  1:7623327/1‑150 (MQ=255)
                                                                                                                                                  atatCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGc  >  1:7623849/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGC  >  NC_002947/2672925‑2673218

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