Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,673,054 G→C 8.6% intergenic (+80/‑87) 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,0540GC 8.6% 585.9 / 32.6 183intergenic (+80/‑87)prpF/prpDaconitate isomerase/2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (69/98);  new base C (4/12);  total (73/110)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 192 total aligned reads displayed.

GGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGC  >  NC_002947/2672905‑2673202
                                                                                                                                                      |                                                                                                                                                    
ggCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCTCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAgg                                                                                                                                                      <  1:3720455/150‑1 (MQ=255)
ggCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAgg                                                                                                                                                      <  2:1809282/150‑1 (MQ=255)
ggCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAgg                                                                                                                                                      <  2:1595303/150‑1 (MQ=255)
   aaTCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGtt                                                                                                                                                   >  2:6018814/1‑150 (MQ=255)
   aaTCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCGTTCCCCCACCCCCCCGGAATAAAAGGGtt                                                                                                                                                   >  1:2872129/1‑150 (MQ=255)
   aaTCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCGTGGCCGCGCCCCCACGTGAACACAggggt                                                                                                                                                   >  1:5250274/1‑148 (MQ=255)
    aTCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGttt                                                                                                                                                  <  2:4454627/150‑1 (MQ=255)
       aTGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGaa                                                                                                                                               <  1:6157953/150‑1 (MQ=255)
        tGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAAt                                                                                                                                              >  2:2744368/1‑150 (MQ=255)
            ccGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtg                                                                                                                                          <  2:4049084/150‑1 (MQ=255)
            ccGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtg                                                                                                                                          <  2:1033741/150‑1 (MQ=255)
             cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                         <  1:1304432/150‑1 (MQ=255)
             cGTAGTGCTCGCGTGCTGATGGAGGGCGGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                                         <  2:6429662/150‑1 (MQ=255)
              gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                                        <  2:4570860/150‑1 (MQ=255)
                aGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGat                                                                                                                                      <  2:2627255/150‑1 (MQ=255)
                    cTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATcc                                                                                                                                  <  2:4808031/150‑1 (MQ=255)
                     tCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCACCCACAGGTAATCACAGGGTTTGAAAGCGGTGATATccc                                                                                                                                 >  1:2805447/1‑150 (MQ=255)
                      cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                <  2:3803938/150‑1 (MQ=255)
                      cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                <  1:5102678/150‑1 (MQ=255)
                      cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                <  2:3803891/150‑1 (MQ=255)
                        cgTGCTGATGGAGGGCTGTGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                              <  2:6337139/150‑1 (MQ=255)
                        cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                              <  2:1094584/150‑1 (MQ=255)
                        cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                              <  2:2010768/150‑1 (MQ=255)
                         gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                             <  1:6320933/150‑1 (MQ=255)
                         gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                             <  1:4443408/150‑1 (MQ=255)
                          tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTgg                                                                                                                            <  2:1354512/150‑1 (MQ=255)
                          tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTgg                                                                                                                            <  2:3788971/150‑1 (MQ=18)
                          tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTgg                                                                                                                            <  2:1699905/150‑1 (MQ=18)
                            cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                          <  1:5957563/150‑1 (MQ=18)
                             tGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCATGCCCCCGCCCACAGGTAAACCCACGGGTTGAATGCTTTGCAATCCCGGTGGGAg                                                                                                                         >  2:1846836/1‑150 (MQ=255)
                              gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGGTTGAATGCAGGGATATCCCTGTGGGAGc                                                                                                                        >  2:3400154/1‑150 (MQ=255)
                               ttGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCg                                                                                                                       <  1:3951940/149‑1 (MQ=255)
                               aTGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCg                                                                                                                       <  1:4922909/150‑1 (MQ=255)
                                tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                      <  1:4989777/150‑1 (MQ=255)
                                 ggAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCCTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCggg                                                                                                                     <  1:5042057/150‑1 (MQ=255)
                                   aGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGtt                                                                                                                                                   <  2:5087259/118‑1 (MQ=255)
                                    gggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCaa                                                                                                                  <  1:5048163/150‑1 (MQ=255)
                                    gggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCaa                                                                                                                  <  1:5715232/150‑1 (MQ=255)
                                    gggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCaa                                                                                                                  <  1:1743894/150‑1 (MQ=255)
                                    gggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTGTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCaa                                                                                                                  <  2:1667512/150‑1 (MQ=17)
                                      gCTGGGTTCGCGTGCCTGGCGATAGCTTCTAACGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGc                                                                                                                <  1:187420/150‑1 (MQ=255)
                                      gCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGc                                                                                                                <  1:4133586/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)
                                             tCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGaa                                                                                                         <  2:1110972/150‑1 (MQ=255)
                                                   gCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCATGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGATCCCACAGTTAATCACAGGTTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTg                                                                                                   <  2:5592854/150‑1 (MQ=255)
                                                     cTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGca                                                                                                 >  2:5681013/1‑150 (MQ=255)
                                                          gATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCCCACGGAAAAACCGGGGTTTTATTCGGTGATATACCGGGGGGCGCGGGCCACCCCCCGCAGGGCGGCCcagca                                                                                            >  1:4194442/1‑150 (MQ=255)
                                                                ttCTAACGCACCACGGGCGCTCGGCTTCAAGCGGGGCCGCTTTGCGGCCGTTCGCGGGCTGGCCCGCTCCCACCGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGTGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCg                                                                                      <  2:6626980/150‑1 (MQ=37)
                                                                 tCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTGTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGTAATCACCGGTTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCgg                                                                                     <  2:389690/150‑1 (MQ=255)
                                                                      aGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAAtt                                                                                <  2:4538768/150‑1 (MQ=255)
                                                                        caccacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGTTATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAAATCt                                                                              >  1:4496689/1‑150 (MQ=255)
                                                                           cacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                           <  2:5018642/150‑1 (MQ=255)
                                                                           cacTGGCGCTCGGATTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGATTGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                           <  2:1551699/150‑1 (MQ=255)
                                                                               ggCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGTAAACACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGa                                                                       <  2:2158485/150‑1 (MQ=255)
                                                                                gcgcTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTGGCCCGCTCCCACCGGTAATCACAGGGTTTGATTGCTGTGATATCCCTGTGGGAGGGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAg                                                                      <  2:3568387/150‑1 (MQ=255)
                                                                                   cTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCt                                                                                                                                           <  2:2325446/78‑1 (MQ=255)
                                                                                   cTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCt                                                                                                                                           >  1:2325446/1‑78 (MQ=255)
                                                                                    tCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCAAACGTAAACACCGGGTTTGAATGCTGGGATATACCAGTGGGAGCCGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTgg                                                                  >  2:2705940/1‑150 (MQ=38)
                                                                                         ttCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCAGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATg                                                             >  2:541535/1‑150 (MQ=255)
                                                                                            aaGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTGGCCCGCTCCCACAGGTAATCACTGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAgc                                                          <  2:3950314/150‑1 (MQ=255)
                                                                                                 gggCCGCTTGGCGGCCCTACGCGGGCTGGCCCGCTCCCCCAGGAAAGCCCAGGGTTTGAAGGCGGTGATATCCCTGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAAGTCTGAAGGGAGCTGGATATGAGCGCCaa                                                     <  2:1424457/150‑1 (MQ=37)
                                                                                                    ccGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGc                                                                                                                        >  1:3067913/1‑80 (MQ=255)
                                                                                                    ccGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCCCAGGGAATCACAGGGGTTGGATGCGGTGTTATACCAGTGGGGGCGGGCCAGCCCGCCCAGGGGGGCACAGCAGCCCCGGCCAATGTGGAGGGGGCGGGGTATGGGCGCCCACGt                                                  >  2:614504/1‑150 (MQ=255)
                                                                                                             ggCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCaa                                         <  2:5712957/150‑1 (MQ=255)
                                                                                                              gCCCTGCGCGGGCTTGCCCGCTCCCACAGTAAAGCACAGGTTTTGATTGTTTTTATACCCCTGTGGGCGGGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAAc                                        <  1:1053987/150‑1 (MQ=255)
                                                                                                                  ttCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGTATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCGCAACGATc                                    >  2:986916/1‑150 (MQ=255)
                                                                                                                  ttCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTATGAAGGGCGCTGGATATGAGCGCCAACGTAGACCTCAACGATc                                    >  2:472685/1‑150 (MQ=255)
                                                                                                                            tGCCCGCTCCCCCAGGTAATCACCGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                          <  2:5943118/150‑1 (MQ=255)
                                                                                                                            ggCCCGCTCCCCCCGGTAATCCCAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGTTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                          <  2:6435317/149‑1 (MQ=255)
                                                                                                                            gccccGCTCCCACCGTAAATCACTGGTTTTTATTCGTTTTATCCCCCTGGGGGCGCGGCCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                          <  1:2965624/148‑1 (MQ=255)
                                                                                                                             gCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCACTTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc                         >  2:5261322/1‑150 (MQ=255)
                                                                                                                             gCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc                         >  1:3233614/1‑150 (MQ=255)
                                                                                                                              cccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACg                        >  2:3694176/1‑150 (MQ=255)
                                                                                                                              cccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACg                        >  1:4605980/1‑150 (MQ=255)
                                                                                                                               ccGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa                       >  1:5655560/1‑150 (MQ=18)
                                                                                                                               ccGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa                       >  1:3572460/1‑150 (MQ=17)
                                                                                                                                     ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                 >  2:5691255/1‑150 (MQ=255)
                                                                                                                                     ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                 >  1:3773720/1‑150 (MQ=255)
                                                                                                                                     ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                 >  2:2958574/1‑150 (MQ=255)
                                                                                                                                     ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                 >  1:3868098/1‑150 (MQ=255)
                                                                                                                                     ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                 >  1:6281966/1‑150 (MQ=255)
                                                                                                                                     ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                 >  1:6109973/1‑150 (MQ=255)
                                                                                                                                      cccaGGTAATCACAGGTTTTGATTGCTGTGCTAGCCCTGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTg                <  2:2872129/148‑1 (MQ=255)
                                                                                                                                       acaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAAc                                        <  1:2669900/125‑1 (MQ=255)
                                                                                                                                       acaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGc               >  2:3140567/1‑150 (MQ=255)
                                                                                                                                         agttAACCACTGGGTTTGAATGCTTTCATATCCCTGTGGGAGCGGGCACCCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTa             <  1:3784701/147‑1 (MQ=255)
                                                                                                                                           gTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTAc                         <  2:2483714/136‑1 (MQ=255)
                                                                                                                                           gTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACa           >  1:5398069/1‑150 (MQ=255)
                                                                                                                                             aaTCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACaaa         >  2:6254138/1‑150 (MQ=255)
                                                                                                                                               aTCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAc        >  1:5280554/1‑150 (MQ=255)
                                                                                                                                                 cactGGGTTTGATTGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAccc      <  2:3669953/150‑1 (MQ=255)
                                                                                                                                                 cactGGGTTTGAAGGCTGTGATCTCCGTGTGGGTGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAccc      <  2:5790264/150‑1 (MQ=255)
                                                                                                                                                  acaGGGTTTGAATGCTGTGATATCCCTGGGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt     >  2:2276267/1‑150 (MQ=255)
                                                                                                                                                    aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt     >  1:5771550/1‑148 (MQ=255)
                                                                                                                                                    aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg   >  1:4713426/1‑150 (MQ=255)
                                                                                                                                                    aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg   >  2:4339499/1‑150 (MQ=255)
                                                                                                                                                     gggTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCACTTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGc  >  2:1741073/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGC  >  NC_002947/2672905‑2673202

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: