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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,670,0440A.12.9% 415.4 / 60.9 139intergenic (+97/‑70)prpF/prpDaconitate isomerase/2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (65/52);  new base . (7/11);  total (72/67)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

Only 100 of 139 total aligned reads displayed.

CTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCG  >  NC_002947/2669898‑2670189
                                                                                                                                                    |                                                                                                                                                 
cTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATcc                                                                                                                                                <  2:2311843/150‑1 (MQ=255)
cTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATcc                                                                                                                                                <  2:107173/150‑1 (MQ=255)
  cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                              <  2:1684286/150‑1 (MQ=255)
  cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                              <  2:874999/150‑1 (MQ=255)
  cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                                              <  2:5767420/150‑1 (MQ=255)
   gcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtg                                                                                                                                             <  2:3782190/150‑1 (MQ=255)
    cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                                            <  1:3445542/150‑1 (MQ=255)
    cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                                            <  2:4901230/150‑1 (MQ=255)
    cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                                            <  1:3842530/150‑1 (MQ=255)
    cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                                            <  1:3320515/150‑1 (MQ=255)
    cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                                            <  2:5755069/150‑1 (MQ=255)
    cgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgt                                                                                                                                            <  2:8224094/150‑1 (MQ=255)
     gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                                           <  1:7331163/150‑1 (MQ=255)
        cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                                        <  2:1207112/150‑1 (MQ=18)
        cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                                                        <  2:1841683/150‑1 (MQ=18)
          gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCACGAGGGGCCGCTTTGCGGCCCATCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGc                                                                                                                                      <  2:3839522/150‑1 (MQ=255)
          gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCCCCGGGAAACACACGGGTTGAAAGCAGAGATATACCAGTGGGAGc                                                                                                                                      >  1:5138150/1‑150 (MQ=255)
          gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGc                                                                                                                                      <  2:8342836/150‑1 (MQ=255)
          gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCAAAAGCACCACTCGCGCTCGGCTTCAAGAGGGGCCGCTTTGCAGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGc                                                                                                                                      <  2:4619858/150‑1 (MQ=17)
            tGGAGGGCTGGGTTCGCGTGCCTGGCGCTAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                                                    <  2:1310739/150‑1 (MQ=18)
            tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                                    <  1:8142914/150‑1 (MQ=255)
            tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                                                    <  1:8142900/150‑1 (MQ=255)
            tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                                                    <  2:5298316/150‑1 (MQ=18)
             ggAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCggg                                                                                                                                   <  2:907383/150‑1 (MQ=255)
             ggAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCggg                                                                                                                                   <  2:907436/150‑1 (MQ=255)
              gAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCCCAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGc                                                                                                                                  <  2:6814019/150‑1 (MQ=18)
                gggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCaa                                                                                                                                <  2:2913745/150‑1 (MQ=255)
                 ggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAg                                                                                                                               <  1:4693896/150‑1 (MQ=255)
                 ggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCCCACGTAAACACACGGTTTGGAATGTGGGTTATCACGGGGGGAGCGGGCacg                                                                                                                               >  1:7034878/1‑148 (MQ=255)
                    tGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGGTATCCCTGTGGGAGCGGGCAAGccc                                                                                                                            >  2:1273294/1‑150 (MQ=255)
                    tGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGccc                                                                                                                            <  2:4287624/150‑1 (MQ=255)
                                  tGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGcac                                                                                                              >  2:5919184/1‑150 (MQ=255)
                                   ggCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTTTGAATGCTGTTATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGcaca                                                                                                             <  2:3749592/150‑1 (MQ=255)
                                      gctAGCTTCTACAGCACAAGTGGCGCTCGGCTTCCAGAGGGGCCGCTTTGCTCCCATTGGCGGGATTCACCGCTCCGACAGTTAACCCATGTGTTTGATTGCTGTTATATACCTGGGGTAGCGGTCATGCCCGCGAAGGGCTGCAcagca                                                                                                          <  2:7406123/148‑1 (MQ=255)
                                       aTAGCTTCTAACGCACCACTGGCGCTCGGCTTCAAGAGGGGCGGCTTGGCGGCCCTTCGCGGGCTGGCCCGCTCCCCCAGGTAATCACTGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCAcagcag                                                                                                         <  2:6557346/150‑1 (MQ=38)
                                           cTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCCCCGGGAATCACCGGGGTTGTATGCGGGGGGATCCCGGGGGGGGCGGGCAAGCCCGCGCAGGGCGGCGCAGCCGcccc                                                                                                     >  1:8903177/1‑150 (MQ=255)
                                             tCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCgg                                                                                                   >  2:4814877/1‑150 (MQ=255)
                                              gccgaGCGCCAGTGGTGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCACCCACAGGTAATAACAGGGTTTGAATTCTGTGATATCCCAGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGc                                                                                                  >  2:190250/5‑150 (MQ=37)
                                              caaaaGCACCACTGGCGCTCGGCTACAAGAGGGGCCGCTTGGCGGCCCTTCGGGGGCTTGCCCGCTCCCACAGTTAACCACAGGTTTTGAATGCTGGGATAACCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGc                                                                                                  <  2:5500684/148‑1 (MQ=37)
                                                aaaGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCCTCGCGGGCTTGCCCGCTCCCACCGGTAATCACAGGGTTTGAATGCTGTTATACCCCGGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGcccttcgcg                                                                                                <  2:2747607/150‑7 (MQ=38)
                                                 aaGCACCAGTGGCGCTCGGCTACAGGAGGGGCCGCTTGGCGGCCCTGCGGGGGGTTGCCCGCTCCCACTGGTAATCACTGGGTTGGAAGGCTGTTATATCCGTGTGGGTGCGGGCAAGCCCGCGAAGGGCTGCTGTGCAGcccttcgcgg                                                                                               <  2:7858575/150‑8 (MQ=255)
                                                 aaGCACCACGGGCGCTCGGCTTCAAGTGGGGCCGCTTGGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCTAGCCCGCGAAGGGCTGCTGTGCAGcccttcgcgg                                                                                               <  2:7858636/150‑8 (MQ=37)
                                                  aGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTTTGAATGCGGTGATATCCCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCaaat                                                                                              >  2:1540986/1‑148 (MQ=255)
                                                      ccacTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTGGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGa                                                                                          <  2:8921123/150‑1 (MQ=255)
                                                       cccgggggcTGGGCTTCACGCGGGGCCGGTTTGGGGCCCTGGGGGGGCTTGCCCGCGCCCACGGTTAAGCCCTGGGTTTGTTTGCTGTTTTTTCCCTGTGGGCGGGGGAACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                                         <  2:2566853/143‑1 (MQ=255)
                                                       cacTGGCGCTCGGATTCATGCGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCACCCACCGTTAATCACAGGTTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                                                         <  2:2567575/150‑1 (MQ=37)
                                                          tGGCTCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAggg                                                                                      <  2:5143410/150‑1 (MQ=255)
                                                                tgggCTTCATGTGGGGCGGCTTTGCGGCCCTTGGCGGGGTTGCCGGCTCCCACTGGTAATCACTGGTTTTGAATGCTGTGATTGCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTgg                                                                                <  1:9024983/148‑1 (MQ=25)
                                                                tCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTgg                                                                                <  2:7340876/150‑1 (MQ=255)
                                                                tCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTgg                                                                                >  2:7459876/1‑150 (MQ=255)
                                                                        atgcggggCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCCCCCACAGTAAATAACCGGTTTTGATTGCTGTCATATCCCGGGGGGCGGGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAgc                                                                        <  1:4466228/146‑1 (MQ=255)
                                                                        aatagGGGCAGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAAAAAAAGGGTTAGAATGCGGTGATATCCCTGTGGGACCCGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAgc                                                                        >  2:5175912/4‑150 (MQ=37)
                                                                                  gCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGTAATCACAGGGTTTGAATGCGGGGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCGGCCCAGCCGCCCCGGCACTTCTGAAGGGAGCTGGATATGTGCGCCAACGTAg                                                              >  2:6716790/1‑150 (MQ=37)
                                                                                  gCTTTGCGGCCCTTCGCGGGCGTGCCCGCTCCCACAGGTAAACACAGGGTTTGAATGCTGTGATATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAAATCTGAAGGGAGCTGGATATGAGCGCCAACGTAg                                                              >  2:6406397/1‑150 (MQ=255)
                                                                                             cTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGAt                                                   >  2:9026722/1‑150 (MQ=255)
                                                                                                  cgGGCTGGCCCGCTCCCCCCGGTACGCCCAGGGTTTGCATGCTGTGATATCCCTGGGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGccc                                              <  1:4475976/150‑1 (MQ=255)
                                                                                                   gggCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCg                                             >  1:2800268/1‑150 (MQ=17)
                                                                                                   gggCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCACAACGATCGCCCg                                             >  1:7333969/1‑150 (MQ=17)
                                                                                                     gCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGa                                           >  2:705014/1‑150 (MQ=255)
                                                                                                     gCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGa                                           >  1:817257/1‑150 (MQ=255)
                                                                                                       ttGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt                                         >  1:1326344/1‑150 (MQ=255)
                                                                                                        tGCCCGCTCCCACAGGTAACCACCGGGTTTGAATGCTGTGATATCCCTGGGGGCGGGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTa                                        <  1:2210728/150‑1 (MQ=255)
                                                                                                          cccGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACg                                      >  1:2837963/1‑150 (MQ=18)
                                                                                                           ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGcccc                                                                                                     <  2:4525817/86‑1 (MQ=255)
                                                                                                           ccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa                                     >  1:4841155/1‑150 (MQ=255)
                                                                                                           ccGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa                                     >  2:6774989/1‑150 (MQ=18)
                                                                                                             gCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGAcc                                   >  1:1294990/1‑150 (MQ=18)
                                                                                                               tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                                 >  2:1616541/1‑150 (MQ=255)
                                                                                                               tCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTATGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                                 >  1:3083138/1‑150 (MQ=255)
                                                                                                               acccACCGTTAATCACTGGTTTTGATTGGTTTTTTACCCCTGTGGTTGGGGGAAACCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCgg                                 <  1:8773172/149‑1 (MQ=255)
                                                                                                                 ccACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt                               >  2:3015507/1‑150 (MQ=255)
                                                                                                                   acaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGccc                                              <  2:2697831/133‑1 (MQ=255)
                                                                                                                   acaGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGc                             >  2:6427687/1‑150 (MQ=255)
                                                                                                                     cggTACTCACAGGGTTTGAAGGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTa                           <  2:8903177/149‑1 (MQ=255)
                                                                                                                       gTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACa                         >  1:5521444/1‑150 (MQ=255)
                                                                                                                             cactGGGTTTGAATGCTGTGATTTCCCTGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAccc                    <  2:831733/150‑1 (MQ=255)
                                                                                                                             cacaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAAccc                    >  1:2364869/1‑150 (MQ=255)
                                                                                                                              acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt                   >  2:4551572/1‑150 (MQ=255)
                                                                                                                              acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt                   >  2:2678249/1‑150 (MQ=255)
                                                                                                                              acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCt                   >  1:7060944/1‑150 (MQ=255)
                                                                                                                              acaGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCAGGTGCTACAAACCCt                   >  1:1980654/1‑150 (MQ=255)
                                                                                                                               ccgggTTTTAATGCTGTGATATCCCTGGGGGTGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg                  <  2:1206670/148‑1 (MQ=255)
                                                                                                                               caGGGTTTGAATGCTGTGATATCCCTGTGTGATCGGGCAAGCCCGCGAAGGGCTGCACATCAGACACGGCAAATATGAAGTGAGCTGGATATGAGCGCCAACGTAGACATCAACGAACGCCCGGATTACTAACGGGTGCTACAAACCCTg                  >  1:1286945/1‑150 (MQ=255)
                                                                                                                               caGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGTTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg                  >  2:4344846/1‑150 (MQ=255)
                                                                                                                               caGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTg                  <  2:5795093/150‑1 (MQ=255)
                                                                                                                                cgggTTCGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg                 <  2:3538806/149‑1 (MQ=255)
                                                                                                                                aGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGGGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTgg                 >  1:7107797/1‑150 (MQ=255)
                                                                                                                                   gTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCg              >  2:4225152/1‑150 (MQ=255)
                                                                                                                                    tttGATTGCTGTTATATCCCTGTGGTAGGGTGCATTCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATTTGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTTGCCGa             <  2:6224383/150‑1 (MQ=255)
                                                                                                                                    tttGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGa             <  2:3191391/150‑1 (MQ=255)
                                                                                                                                    tttGAATGCTGGGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCgt             >  1:9135770/1‑149 (MQ=255)
                                                                                                                                          acccTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGc       >  2:850142/4‑150 (MQ=255)
                                                                                                                                             cTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGccc     >  2:4422082/1‑150 (MQ=255)
                                                                                                                                              tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt    >  2:7086417/1‑150 (MQ=255)
                                                                                                                                              tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt    >  1:6073919/1‑150 (MQ=255)
                                                                                                                                              tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt    >  1:3814950/1‑150 (MQ=255)
                                                                                                                                              tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt    >  2:4984221/1‑150 (MQ=255)
                                                                                                                                              tgtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCt    >  1:4487259/1‑150 (MQ=255)
                                                                                                                                               gtgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTTCGCCCTc   >  2:4644898/1‑150 (MQ=255)
                                                                                                                                                tgATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCg  >  1:8826126/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCTACAAACCCTGGCCGATTACGCCCTCG  >  NC_002947/2669898‑2670189

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: