Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,673,062 T→C 8.6% intergenic (+88/‑79) 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,0620TC 8.6% 520.9 / 27.0 178intergenic (+88/‑79)prpF/prpDaconitate isomerase/2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (66/93);  new base C (3/12);  total (71/107)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-01

Only 100 of 182 total aligned reads displayed.

TGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACG  >  NC_002947/2672913‑2673210
                                                                                                                                                      |                                                                                                                                                    
tGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAAt                                                                                                                                                      <  2:4583029/150‑1 (MQ=255)
   gCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCt                                                                                                                                                   <  2:2314355/150‑1 (MQ=255)
    ccGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtg                                                                                                                                                  <  2:1033741/150‑1 (MQ=255)
     cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                                 <  1:1304432/150‑1 (MQ=255)
     cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                                 <  2:1785597/150‑1 (MQ=255)
     cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                                 <  2:1668060/150‑1 (MQ=255)
     cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                                 <  1:4105748/150‑1 (MQ=255)
     cGTAGTGCTCGCGTGCTGATGGAGGGCGGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                                 <  2:6429662/150‑1 (MQ=255)
      gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                                                <  1:5388809/150‑1 (MQ=255)
       tAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGa                                                                                                                                               <  1:2925446/150‑1 (MQ=255)
       tAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCGTGGCCGCGCCCCCCCGTGAAAACAGGGTTTGGATGCGGTGa                                                                                                                                               >  1:4143562/1‑150 (MQ=255)
       tAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCGTGCCCCCACCCCCACGTAAAACCCGCGGTTGAATGCtgtgt                                                                                                                                               >  1:5192586/1‑149 (MQ=255)
        aGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGat                                                                                                                                              <  2:2627255/150‑1 (MQ=255)
            cTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATcc                                                                                                                                          <  2:4808031/150‑1 (MQ=255)
             tCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCACCCACAGGTAATCACAGGGTTTGAAAGCGGTGATATccc                                                                                                                                         >  1:2805447/1‑150 (MQ=255)
              cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                        <  1:5752760/150‑1 (MQ=255)
              cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                        <  2:3803891/150‑1 (MQ=255)
              cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                        <  2:3803938/150‑1 (MQ=255)
                cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                                      <  2:1094557/150‑1 (MQ=255)
                cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                                      <  2:365617/150‑1 (MQ=255)
                 gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                                     <  1:4443408/150‑1 (MQ=255)
                  tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGGCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTgg                                                                                                                                    <  2:1509608/150‑1 (MQ=18)
                  tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTgg                                                                                                                                    <  2:3788971/150‑1 (MQ=18)
                  tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTgg                                                                                                                                    <  2:1710024/150‑1 (MQ=18)
                  tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTgg                                                                                                                                    <  2:5085372/150‑1 (MQ=18)
                    cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                                  <  1:901237/150‑1 (MQ=18)
                     tGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCATGCCCCCGCCCACAGGTAAACCCACGGGTTGAATGCTTTGCAATCCCGGTGGGAg                                                                                                                                 >  2:1846836/1‑150 (MQ=255)
                       ttGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCg                                                                                                                               <  1:3951940/149‑1 (MQ=255)
                       aTGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCg                                                                                                                               <  1:4922909/150‑1 (MQ=255)
                       aTGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCg                                                                                                                               <  1:5694320/150‑1 (MQ=18)
                        tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                              <  1:5819465/150‑1 (MQ=255)
                            gggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCaa                                                                                                                          <  1:5048163/150‑1 (MQ=255)
                            gggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCaa                                                                                                                          <  1:1743894/150‑1 (MQ=255)
                            gggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCaa                                                                                                                          <  2:1097132/150‑1 (MQ=255)
                              gCTGGGTTCGCGTGCCTGGCGATAGCTTCTAACGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGc                                                                                                                        <  1:187420/150‑1 (MQ=255)
                               cTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGcc                                                                                                                       <  2:6417252/150‑1 (MQ=255)
                               cTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGGGATATACCGGTGGGAGCGGGCAAGcc                                                                                                                       >  2:1536915/1‑150 (MQ=255)
                                     tCGCGTGCGGGGCTATAGCTCCTAAAGCACCAGTGGCGGTCGGCTTCAAGTGGGGCCGCTTTGCGGCCCTTCGGGGGTTTGCCCGCTCCCACTGGTAATCACTGGTTTTGATTGGTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGaa                                                                                                                 <  2:4596375/150‑1 (MQ=37)
                                     tCGCGTGCCTGGCGATTGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGaa                                                                                                                 <  2:468351/150‑1 (MQ=255)
                                     tCGCGTGCCTGGCGATAGCTTCTACAGCACCACGGGGGCTCGGCTTCAAGCGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGTGCGGGCAAGCCCGCGaa                                                                                                                 <  2:6526399/150‑1 (MQ=38)
                                        cgTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTGGGATTAATGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATAACAGGGTTTGATTGCTGTTATTTCCCTTTGGGAGCGGGAAAGCCCGCGAAggg                                                                                                              <  2:1403456/150‑1 (MQ=37)
                                           gCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCATGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGATCCCACAGTTAATCACAGGTTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTg                                                                                                           <  2:5592854/150‑1 (MQ=255)
                                             cTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGca                                                                                                         >  2:440882/1‑150 (MQ=255)
                                             cTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGca                                                                                                         >  2:5681013/1‑150 (MQ=255)
                                                  gATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCCCACGGAAAAACCGGGGTTTTATTCGGTGATATACCGGGGGGCGCGGGCCACCCCCCGCAGGGCGGCCcagca                                                                                                    >  1:4194442/1‑150 (MQ=255)
                                                         tCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTGTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGTAATCACCGGTTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCgg                                                                                             <  2:389690/150‑1 (MQ=255)
                                                          cTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGATTGCCCGCTCCCACAGTTAATCACAGGTTTTTAATGCTGTTATATCCCTGTGGGTGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGc                                                                                            <  2:4279853/150‑1 (MQ=255)
                                                              aGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAAtt                                                                                        <  2:4538768/150‑1 (MQ=255)
                                                                caccacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGGGTGCCCGCTCCCCCAGGTCATCAAAGGGGTTGAATTGGGGTTTATACCGGGGGGCGGGGGGCAGCCCGCCAAGGGCGGCACAGCCGCCCCCGCCAGTCt                                                                                      >  1:2713847/1‑150 (MQ=255)
                                                                caccacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGTTATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAAATCt                                                                                      >  1:4496689/1‑150 (MQ=255)
                                                                   cacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCCCACGGAATCCCACGGTTTTAAAGCGGTGATAAACCGGGGGGAGCCGGCAAGCCCCCGAAGGGCGGCACCGCCGCCCCCGCAAGTCGGaa                                                                                   >  2:630003/1‑150 (MQ=255)
                                                                   cacTGGCGCTCGGATTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGATTGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                                   <  2:1551699/150‑1 (MQ=255)
                                                                    acTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCACAGGGAATCACAGGGTTTGAATGCTGTGATATCCCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCCGCACTTCTGGAg                                                                                  >  2:2417346/1‑150 (MQ=37)
                                                                            tCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCAAACGTAAACACCGGGTTTGAATGCTGGGATATACCAGTGGGAGCCGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTgg                                                                          >  2:2705940/1‑150 (MQ=38)
                                                                                 ttCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCAGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATg                                                                     >  2:541535/1‑150 (MQ=255)
                                                                                    aaGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTGGCCCGCTCCCACAGGTAATCACTGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAgc                                                                  <  2:3950314/150‑1 (MQ=255)
                                                                                            ccGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGc                                                                                                                                >  1:3067913/1‑80 (MQ=255)
                                                                                                tttGCGGCCCTCCGCGGGCTTGCCCGCTCCCCCCGTTAATCACCGGTTTTTAATGCTTTTATACCCGTGGGGGAGCGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGAc                                                      <  1:6018814/150‑1 (MQ=255)
                                                                                                     ggCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCaa                                                 <  2:5712957/150‑1 (MQ=255)
                                                                                                     ggCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCaa                                                 <  2:2247230/150‑1 (MQ=255)
                                                                                                      gCCCTGCGCGGGCTTGCCCGCTCCCACAGTAAAGCACAGGTTTTGATTGTTTTTATACCCCTGTGGGCGGGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAAc                                                <  1:1053987/150‑1 (MQ=255)
                                                                                                         cTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGCATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGAt                                             >  2:2786025/1‑150 (MQ=255)
                                                                                                          ttCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGTATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCGCAACGATc                                            >  2:986916/1‑150 (MQ=255)
                                                                                                          ttCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTATGAAGGGCGCTGGATATGAGCGCCAACGTAGACCTCAACGATc                                            >  2:472685/1‑150 (MQ=255)
                                                                                                                   ttGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt                                   >  2:401593/1‑150 (MQ=255)
                                                                                                                    tGCCCGCTCCCCCAGGTAATCACCGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                                  <  2:5943118/150‑1 (MQ=255)
                                                                                                                    ggCCCGCTCCCCCCGGTAATCCCAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGTTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                                  <  2:6435317/149‑1 (MQ=255)
                                                                                                                    ggCCCGCTCCCCCAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                                  <  2:1597050/149‑1 (MQ=255)
                                                                                                                     gCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc                                 >  1:3233614/1‑150 (MQ=255)
                                                                                                                      cccGCTCCCACAGGTAATCACTGGTTTTTATTGCTGTGATATCCGTGTGGGAGGGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACg                                <  1:4993082/150‑1 (MQ=255)
                                                                                                                      cccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACg                                >  2:3694176/1‑150 (MQ=255)
                                                                                                                      cccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACg                                >  1:4605980/1‑150 (MQ=255)
                                                                                                                       ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                                   >  1:5588287/1‑98 (MQ=255)
                                                                                                                       ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa                               >  2:2056663/1‑150 (MQ=255)
                                                                                                                       ccGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa                               >  1:3572460/1‑150 (MQ=17)
                                                                                                                         gCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGAcc                             >  2:6145723/1‑150 (MQ=255)
                                                                                                                           tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                           >  2:3685987/1‑150 (MQ=255)
                                                                                                                           tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                           >  2:1182953/1‑150 (MQ=255)
                                                                                                                           gcccACCGGTAATCACTGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                           <  2:5250274/149‑1 (MQ=255)
                                                                                                                             ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                         >  2:5691255/1‑150 (MQ=255)
                                                                                                                             ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                         >  1:3868098/1‑150 (MQ=255)
                                                                                                                             ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                         >  1:6281966/1‑150 (MQ=255)
                                                                                                                             ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                         >  1:3773906/1‑150 (MQ=255)
                                                                                                                             ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                         >  1:6109973/1‑150 (MQ=255)
                                                                                                                             ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                         >  2:2958574/1‑150 (MQ=255)
                                                                                                                              cccaGGTAATCACAGGTTTTGATTGCTGTGCTAGCCCTGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTg                        <  2:2872129/148‑1 (MQ=255)
                                                                                                                               acaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGc                       >  2:3140567/1‑150 (MQ=255)
                                                                                                                                 agttAACCACTGGGTTTGAATGCTTTCATATCCCTGTGGGAGCGGGCACCCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTa                     <  1:3784701/147‑1 (MQ=255)
                                                                                                                                  ggTAATCACTGGTTTTGAATGCTGTTATATCCCTGTGGGAGCGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTAc                    <  2:3418373/150‑1 (MQ=255)
                                                                                                                                   gTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                                  <  1:4979809/135‑1 (MQ=255)
                                                                                                                                   gTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                                  >  2:4979809/1‑135 (MQ=255)
                                                                                                                                   gTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc                                 >  1:2483714/1‑136 (MQ=255)
                                                                                                                                   gTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc                                 <  2:2483714/136‑1 (MQ=255)
                                                                                                                                       aTCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAc                >  1:5280554/1‑150 (MQ=255)
                                                                                                                                         cactGGGTTTGAAGGCTGTGATCTCCGTGTGGGTGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAccc              <  2:5790264/150‑1 (MQ=255)
                                                                                                                                           ccggttttGATTGCTGTTATTGCCCGGTGGGTGGGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg            <  1:5984366/145‑1 (MQ=255)
                                                                                                                                            aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg           >  1:4713426/1‑150 (MQ=255)
                                                                                                                                                tttGATTGCTGTTATATCCCTGTGGTAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGa       <  2:2980082/150‑1 (MQ=255)
                                                                                                                                                 ttGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGAt      >  1:6571937/1‑150 (MQ=255)
                                                                                                                                                     aTGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACg  <  2:6789944/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACG  >  NC_002947/2672913‑2673210

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: