Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,673,073 C→A 12.6% intergenic (+99/‑68) 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,0730CA12.6% 492.0 / 44.7 177intergenic (+99/‑68)prpF/prpDaconitate isomerase/2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (87/66);  new base A (14/10);  total (101/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 182 total aligned reads displayed.

CTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGCTACC  >  NC_002947/2672925‑2673222
                                                                                                                                                      |                                                                                                                                                     
cTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATcc                                                                                                                                                      <  1:2486685/150‑1 (MQ=255)
 tCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCACCCACAGGTAATCACAGGGTTTGAATGCGGTGATATacc                                                                                                                                                     >  2:2206299/1‑150 (MQ=255)
  cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCCCACCCACACGGAAACAAACCGGTTTAAAGCTGTGCTATcacg                                                                                                                                                    >  1:6614064/1‑147 (MQ=255)
   gcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtg                                                                                                                                                   <  2:8668073/150‑1 (MQ=255)
     gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGCATTCAAACCCTGTGATTACCTGTgg                                                                                                                                                <  1:4916224/150‑1 (MQ=25)
      tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCTGCTGTGCAGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTgg                                                                                                                                                <  2:6346370/150‑1 (MQ=17)
      tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTgg                                                                                                                                                <  1:5187294/150‑1 (MQ=255)
        cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                                              <  2:8031394/150‑1 (MQ=18)
         tGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGGTTGAATGCTGTGATATCCCGGTGGGAg                                                                                                                                             >  2:1478553/1‑150 (MQ=255)
          gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                                                          <  1:1777545/150‑1 (MQ=21)
            tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                                          <  1:3084381/150‑1 (MQ=255)
            tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                                          <  2:7377570/150‑1 (MQ=255)
            tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                                                          <  2:7818233/150‑1 (MQ=18)
            tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCCCGCCCCCACGGAAACACAGGGTTTGAATGCGGGGCTATCCCGGGGGGAGCgg                                                                                                                                          >  1:7107972/1‑150 (MQ=255)
              gAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGc                                                                                                                                        <  2:1980104/150‑1 (MQ=255)
                gggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAACGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCaa                                                                                                                                      <  1:8449119/150‑1 (MQ=17)
                  gCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGc                                                                                                                                    <  1:6458850/150‑1 (MQ=255)
                       gTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATTCTGTGATATACCGGTGGGAGCGGGCAAGCCcgcg                                                                                                                               >  2:1061464/1‑150 (MQ=255)
                       gTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGTGCGGGCAAGCCcgcg                                                                                                                               >  2:8386696/1‑150 (MQ=255)
                           gcgTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCACCCACAGGTAATCACAGGGGTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGaacg                                                                                                                           >  2:1213880/1‑148 (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)
                                    gCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                                                 >  1:2432908/1‑119 (MQ=255)
                                    gCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                                                 <  2:2432908/119‑1 (MQ=255)
                                     cGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCACCCACAGGTAATCCCAGGGTTTGAATGCTGTGATATCCCGGTGGGAGCGGGCAAGCCCGCGAAGGGCGGCAcagc                                                                                                                 >  2:2198684/1‑150 (MQ=255)
                                         aGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCACAGGTAATCACAGGGGTTGTATGCTGTGATATCCCTGGGGGAGCGGGCAAGCCCGCGAAGGGCTGCACCGCAGcc                                                                                                             >  2:3370791/1‑150 (MQ=38)
                                           cTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCACAGGTAATAACAGGGTTTTTATGCTGTGATATACCGGTGGGAGCGGGCCAGCCCCCGAAGGGCTGCACAGCAGcccc                                                                                                           >  2:2397193/1‑150 (MQ=37)
                                             tCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTGGCCCGCTCCCACAGGTAATCACAGGGTTGGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCgg                                                                                                         <  2:539944/150‑1 (MQ=255)
                                               tAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCCCAGGTAATCACAGGGGTTGAATGCTGTGATATCCCGGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCa                                                                                                       >  2:1582688/1‑150 (MQ=255)
                                                aaaGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGGTTTAAAGCTGTGATAAACCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCCGCaa                                                                                                      >  1:7932309/1‑150 (MQ=37)
                                                   gCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTc                                                                                                   >  2:2700757/1‑150 (MQ=255)
                                                       caggggCGCTCGGCTACACGGGGGGCGGGTTGGCGGCCCTGGGGGGGCTGGCCGGCGCCCACCGTTAATCACCGGGTTTTAAGGGTGTGATATCCGTGTGGTAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                                               <  1:7119329/146‑1 (MQ=255)
                                                       cacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                                          <  1:5813098/107‑1 (MQ=255)
                                                       cacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                                          <  1:5812365/107‑1 (MQ=255)
                                                       cacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                                          >  2:5813098/1‑107 (MQ=255)
                                                                 cGGCTGCAAGAGGGGCCGGTTTGCGGCCCTTGGCGGGCTTGCCCGCTCCCACAGTTAATCACAGGGTTTGAATGCTGTTATATCCCTGTGGGAGCGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGa                                                                                     <  2:6614064/150‑1 (MQ=38)
                                                                     ttCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCACCCACAGGTAATCACAGGGTTTGAATGCTGTGATATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGGGCTGGATATg                                                                                 >  2:6215556/1‑150 (MQ=255)
                                                                        acgagGGGCCGCTTGGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGTAATCACAGGGTTTGAATGCTGTGTTATCCCTGGGGTAGGGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAgc                                                                              <  2:1886513/148‑1 (MQ=38)
                                                                               gCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACg                                                                       >  2:5119573/1‑150 (MQ=255)
                                                                                 cGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACACGTAAACACAGGGTTTGAATGCTGTGATATCCCGGGGGGGGCGGGCACCCCCCCCAACGGCTGCACAGCAGCCCCGGCAAGTCTGAAGGGCGCTGGCTATGAGCCCCAACGtc                                                                     >  2:5468455/1‑149 (MQ=255)
                                                                                  gCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGGGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGCGCTGGATATGAGCGCCAACGTAg                                                                    >  2:7880655/1‑150 (MQ=255)
                                                                                  gCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCGGGGAGATCCCGGGGGGAGCGGGCCAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCCCTTCTGACGGGAGCTGGATATGAGCGCCAACGtcg                                                                    >  2:7880397/1‑148 (MQ=37)
                                                                                    tttGCGGCCCTTCGGGGGCTGGCCCGCTCCCACCGCTACACACCGGTTTTGAATGGGGTTATACCCCGGTGGGCGGGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGAc                                                                  <  1:3370791/150‑1 (MQ=255)
                                                                                    ttggCGGCCCTTGGGGGGCTTGCCCGCACCCCCATGTAATCACCGGGTTTGAATGCTTTGATATCCCTGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGAc                                                                  <  1:3504004/147‑1 (MQ=255)
                                                                                    tgggCGGCCCTTCGCGGGCTTGCCCGCTCCCCCCGGTAATCCCAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGAc                                                                  <  2:7023825/147‑1 (MQ=255)
                                                                                      tGCGCCCCTTCGGGGGCTGCCCCGCCCCCCCCTTAAATCCCGGGTTTTGATTGCTGTCATATCCGTGGGGGCGCGGCCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCt                                                                <  1:8239793/150‑1 (MQ=255)
                                                                                             cTTCGCGGGCTTGCCCGCGCCCACAGGGCTATCACAGCATTCAAACCCGGTGATTACCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCCACGTAGACCTCAACGAt                                                         >  2:506483/1‑150 (MQ=14)
                                                                                                     ggttGCCGGCTCCCACTGTTAATCACTGTTTTTTATTGTTGTTATTGCCGTGTGGGTGGGGGAACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGa                                                 <  1:5325411/148‑1 (MQ=255)
                                                                                                     gcgggCCCGCGCCCCCAGTTATTCACCGGGTTTGAAGGCTGTGCTAGCCCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGa                                                 <  2:3654772/146‑1 (MQ=255)
                                                                                                       ttGCCCGCTCCCCCAGGTAATCACAGGGTTTGAATGCTGTGATATCCCGGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCACTTCGGAAGGGAGCTGGATATGAGCGCCACCGTAGACCTCACCGATCGCCCGGCtt                                               >  2:3049875/1‑150 (MQ=38)
                                                                                                       ttGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATCTGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt                                               >  2:4841171/1‑150 (MQ=255)
                                                                                                       ttGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt                                               >  2:2882792/1‑150 (MQ=255)
                                                                                                        tGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                                              <  2:6046271/150‑1 (MQ=255)
                                                                                                         gCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc                                             <  2:4286743/150‑1 (MQ=255)
                                                                                                         gCCCGCTCCCACAGGTAAACACTGGGTTTGATTGCTTTGATAGCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc                                             <  2:449336/150‑1 (MQ=39)
                                                                                                           ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa                                           >  1:332676/1‑150 (MQ=255)
                                                                                                           ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa                                           >  2:1070832/1‑150 (MQ=255)
                                                                                                           ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa                                           >  1:6057027/1‑150 (MQ=255)
                                                                                                           ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa                                           >  2:8211635/1‑150 (MQ=255)
                                                                                                             gCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGAcc                                         >  1:2621857/1‑150 (MQ=255)
                                                                                                               tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                                       <  2:2105376/150‑1 (MQ=255)
                                                                                                                 ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATCTGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                                     >  2:4639919/1‑150 (MQ=255)
                                                                                                                 ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                                     >  1:4497434/1‑150 (MQ=255)
                                                                                                                 ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                                     >  2:8181808/1‑150 (MQ=255)
                                                                                                                 ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                                     >  1:2054438/1‑150 (MQ=255)
                                                                                                                 ccACAGGTAATCACAGGGTTTGAATGCTATGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                                     >  1:7692554/1‑150 (MQ=255)
                                                                                                                  cacaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTg                                    >  2:6976076/1‑150 (MQ=255)
                                                                                                                               acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt                         >  1:7589261/1‑150 (MQ=255)
                                                                                                                               acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt                         >  2:476191/1‑150 (MQ=255)
                                                                                                                               acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt                         >  2:2528504/1‑150 (MQ=255)
                                                                                                                                ctggttttGAATGCTGTGTTATCCCTGTGGGTGGGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg                        <  2:5942817/145‑1 (MQ=255)
                                                                                                                                ctgggTTGGAAGGCGGGGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTCCAAACCCTg                        <  2:3718132/148‑1 (MQ=39)
                                                                                                                                caGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg                        >  2:7737105/1‑150 (MQ=255)
                                                                                                                                 aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg                       >  2:1561116/1‑150 (MQ=255)
                                                                                                                                 aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg                       >  2:4120508/1‑150 (MQ=255)
                                                                                                                                   ggTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGAAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGcc                     >  2:5199971/1‑150 (MQ=255)
                                                                                                                                     tttGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGa                   <  2:3852009/150‑1 (MQ=255)
                                                                                                                                     tttGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGa                   >  1:3852009/1‑150 (MQ=255)
                                                                                                                                        gAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTa                >  2:1267923/1‑150 (MQ=255)
                                                                                                                                        aaacccTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTAc               >  1:1274694/6‑150 (MQ=255)
                                                                                                                                          aggCTGTGATATCCCTGTGGGAGCGGGCACGCCCGCGACGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACg              <  2:5593810/148‑1 (MQ=255)
                                                                                                                                          aTGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACg              >  1:6401428/1‑150 (MQ=255)
                                                                                                                                           tGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGc             >  1:2920164/1‑150 (MQ=255)
                                                                                                                                            gCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGcc            >  2:1778892/1‑150 (MQ=255)
                                                                                                                                            ccTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGccc           >  2:320806/2‑150 (MQ=255)
                                                                                                                                             cTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGccc           >  1:1326323/1‑150 (MQ=255)
                                                                                                                                              tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt          >  2:6369635/1‑150 (MQ=255)
                                                                                                                                              tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt          >  1:4062531/1‑150 (MQ=255)
                                                                                                                                              tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt          >  1:6820767/1‑150 (MQ=255)
                                                                                                                                              tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt          >  1:6658456/1‑150 (MQ=255)
                                                                                                                                               tttatatCCTTGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTc         <  1:3438999/147‑1 (MQ=255)
                                                                                                                                                ttatatCCGTGTGGGAGCGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCg        <  1:512316/148‑1 (MQ=255)
                                                                                                                                                tgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCg        >  1:147168/1‑150 (MQ=255)
                                                                                                                                                 gATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCgg       >  1:4618721/1‑150 (MQ=255)
                                                                                                                                                  atatCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGc      >  1:7623849/1‑150 (MQ=255)
                                                                                                                                                     tCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGCTAc   >  2:2271757/1‑150 (MQ=255)
                                                                                                                                                      cccTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGCTAcc  >  2:7398808/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCGGCTACC  >  NC_002947/2672925‑2673222

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: