Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,657,604 T→G 23.5% intergenic (‑86/+66) epd ← / ← tktA D‑erythrose 4‑phosphate dehydrogenase/transketolase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,657,6040TG23.5% 491.1 / 171.0 249intergenic (‑86/+66)epd/tktAD‑erythrose 4‑phosphate dehydrogenase/transketolase A
Reads supporting (aligned to +/- strand):  ref base T (96/93);  new base G (29/29);  total (125/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-01

Only 100 of 253 total aligned reads displayed.

GCAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGC  >  NC_002947/5657455‑5657751
                                                                                                                                                      |                                                                                                                                                    
gcAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCCCGCCCACACGGttt                                                                                                                                                      >  1:3461550/1‑150 (MQ=255)
    gACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTTCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCAt                                                                                                                                                   <  1:3780276/150‑1 (MQ=255)
      cacaGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATga                                                                                                                                                 >  2:5300267/1‑149 (MQ=255)
       acaGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGGCGGTGTCGCCTTCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGaaa                                                                                                                                               <  1:4091723/150‑1 (MQ=21)
       acaGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCaa                                                                                                                                               <  2:1442730/150‑1 (MQ=255)
        caGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAAc                                                                                                                                              <  1:1829701/150‑1 (MQ=255)
         aGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAAcc                                                                                                                                             >  2:8096811/1‑150 (MQ=255)
          gCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACct                                                                                                                                            <  1:999883/150‑1 (MQ=255)
             gCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGGCGGTGTCGCCTTCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAAccctgt                                                                                                                                         <  1:8594334/150‑5 (MQ=18)
              ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                        <  2:2209972/150‑1 (MQ=255)
              ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                        <  1:6067171/150‑1 (MQ=255)
              ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                        <  1:6635605/150‑1 (MQ=255)
              ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                        <  1:6312249/150‑1 (MQ=255)
              ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                        <  1:7218303/150‑1 (MQ=255)
               cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgtg                                                                                                                                       <  2:5686266/150‑1 (MQ=255)
               cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgtg                                                                                                                                       <  1:7985427/150‑1 (MQ=255)
                gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTgg                                                                                                                                      >  2:7202457/1‑150 (MQ=255)
                gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTgg                                                                                                                                      <  1:5575458/150‑1 (MQ=255)
                 aTACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGGCGGTGTCGCCTTCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTggca                                                                                                                                    <  1:7687757/150‑3 (MQ=17)
                 aTACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGGCGGTGTCGCCTTCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGa                                                                                                                                    <  1:5143380/150‑1 (MQ=17)
                  tACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGa                                                                                                                                    <  2:7003436/150‑1 (MQ=38)
                  tACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGa                                                                                                                                    <  2:4225419/150‑1 (MQ=38)
                  tACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGa                                                                                                                                    <  1:6369670/150‑1 (MQ=255)
                   aCGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAg                                                                                                                                   <  2:1630818/150‑1 (MQ=255)
                    cGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGc                                                                                                                                  <  1:3881316/150‑1 (MQ=255)
                    cGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGc                                                                                                                                  <  1:5752043/150‑1 (MQ=255)
                        cATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGc                                                                                                                              <  1:5240065/150‑1 (MQ=255)
                         aTAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGGCGGTGTCGCCTTCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCtt                                                                                                                            <  2:1861227/150‑1 (MQ=17)
                          tAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCtt                                                                                                                            <  2:1778079/150‑1 (MQ=255)
                              cGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGccc                                                                                                                        >  2:8412152/1‑150 (MQ=255)
                                ttGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCcgc                                                                                                                      >  2:9037905/1‑150 (MQ=255)
                                ttGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCcgc                                                                                                                      >  2:1141693/1‑150 (MQ=255)
                                     tGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCCCAGGGGTTGGTGCAACCGCGGTGGGAGCGGGCGTGCCCGCGAATg                                                                                                                 >  1:1659812/1‑150 (MQ=38)
                                     tGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCCCTCCCACAGGGTTTAAAGACACCACTGGTGGAGCGGGGTTTCCCGCGAAta                                                                                                                 >  1:2690309/1‑149 (MQ=255)
                                       cAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGaa                                                                                                              <  2:6473339/150‑1 (MQ=255)
                                          cTTGGTTGGGAGGGTGGGGCATGCGGGGCTCACGAGACCGTTCAATTGTGGTCCGCCTGCGGAATGTGGGGTGTGGCATTCGGGGGCAAGCCCGCTCCCACCGCGGTTGCATGAAACCCTGGGGGAGCGGGCTTGCCCGCGAATgccaca                                                                                                            <  2:1258890/150‑6 (MQ=12)
                                                 gggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACAccg                                                                                                    >  2:3434333/1‑150 (MQ=255)
                                                   gACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGTGGGCTTGCCCGCGAATGAAGGCGACAccgcc                                                                                                  >  2:2127938/1‑150 (MQ=255)
                                                      ggTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACAccgcccgc                                                                                               >  2:6374285/1‑150 (MQ=255)
                                                       gTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCa                                                                                              <  1:6581556/150‑1 (MQ=255)
                                                       gTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCACGGTGGGAGCGGGCTTGCCCGCGCATGAAGGCGACAccgcccgcc                                                                                              >  1:7373427/1‑149 (MQ=255)
                                                        tGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCGGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAg                                                                                             >  2:7945833/1‑150 (MQ=255)
                                                         ggggCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCa                                                                                                                                                    <  1:1024522/95‑1 (MQ=255)
                                                          gggCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGATTTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGTTTTCATGCAACCTCTGTGGGTGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGAt                                                                                           <  1:924318/150‑1 (MQ=255)
                                                            gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTGATGCAACCTCAGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCCGATCt                                                                                         >  1:538442/1‑150 (MQ=255)
                                                            gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCGCTGTGGGAGCGGGCTTGCCCGCGAAAGAAGGCGACACCGCCCGCAGATCt                                                                                         >  1:6044145/1‑150 (MQ=255)
                                                                 cGGGGCTCACGAGACAGTTCAATTGGCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCACGCCCGCTCCCACTGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTaaa                                                                                    <  1:2598787/150‑1 (MQ=255)
                                                                  ggggCTCACGAGACTGTTCAATTGTCGTCAGCCTGCGGATTTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAAt                                                                                   <  1:8861994/150‑1 (MQ=255)
                                                                    ggCTCACGAGACAGTTCATTTGTCGTCAGCCTGCGGAATTTGCGGTGGGGCATTCGCGGGCAAGCCCGCTCCCACTGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCa                                                                                 <  1:1610972/150‑1 (MQ=255)
                                                                    ggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCGGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCGTTAAATCa                                                                                 >  2:2061790/1‑150 (MQ=255)
                                                                       tCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTc                                                                              >  2:6023952/1‑150 (MQ=255)
                                                                         aCGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCCCGCCCCCAGGGTTTAATGCAAACACAGGGGGGGCGGGCTTGCCCCCCCATGAAGGCGACACCCCCCCCAGATCTTTAAAACAGTcat                                                                            >  1:1259057/1‑148 (MQ=255)
                                                                           gagaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCggg                                                                                                                               >  2:2416510/1‑98 (MQ=18)
                                                                           gagaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTcc                                                                          <  1:5699559/150‑1 (MQ=255)
                                                                                 gTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCAACCACAGGGTTTCAAGCAACCTCTGTGGGAGCGGGCATGCCCGCGAATGAAGGCGACACCGCCCGCAGATATTTAAATCAGTAATCCagcagc                                                                    >  2:1193552/1‑150 (MQ=255)
                                                                                  ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCACTGGGGGAGCGGGCGTGGCCCCCCAAGAAGGCGACACCGCCCGCAGATCTTTAAATAAGTCTTCCagcagca                                                                   >  1:2002600/1‑149 (MQ=38)
                                                                                  ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCGCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCCGATCTTTAAATCAGTCTTCCcgcagca                                                                   >  1:5280374/1‑149 (MQ=255)
                                                                                  ttCAATTGTCCTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGct                                                                   >  2:8279180/1‑150 (MQ=255)
                                                                                    cAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGctct                                                                 <  2:5163742/150‑1 (MQ=255)
                                                                                     aaTTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCtt                                                                >  1:4466046/1‑150 (MQ=255)
                                                                                       ttGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCg                                                              <  1:4235727/150‑1 (MQ=255)
                                                                                       ttGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCg                                                              <  1:4365106/150‑1 (MQ=255)
                                                                                          tcgtcAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCa                                                           <  2:1735271/150‑1 (MQ=255)
                                                                                            gtcAGCCTGCGGAATTTGCGGTGTGGCGTTCGCGGGCAAGCCCCCTCCCACAGGGTTTGATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACAACGCCCGCAGATCTTTAAAAAAGTATTCCAGCAGCTCGTAGGCAGt                                                         >  1:7605405/1‑150 (MQ=37)
                                                                                              cAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTgg                                                                                                                                      >  2:3325865/1‑72 (MQ=14)
                                                                                              cAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTgg                                                                                                                                      <  1:3325865/72‑1 (MQ=14)
                                                                                                  cTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAg                                                   <  2:7409369/150‑1 (MQ=255)
                                                                                                   tGCGGAATTTGCGGTGTGGCATTCGCGCGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtg                                                                                                                                         <  2:9053107/64‑1 (MQ=255)
                                                                                                             gCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCa                                        <  2:6832922/150‑1 (MQ=255)
                                                                                                               ggTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg                                      <  2:4735888/150‑1 (MQ=255)
                                                                                                                    ggCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGaa                                 <  2:2690504/150‑1 (MQ=255)
                                                                                                                                   gCCCGCTCCCCCAGGGTTCCTTCCACCCTCTGTGGGCGCGGGCGTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGaa                  <  1:3965610/150‑1 (MQ=255)
                                                                                                                                     ccGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACa                >  2:6238909/1‑150 (MQ=255)
                                                                                                                                     ccGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACa                >  2:4648626/1‑150 (MQ=255)
                                                                                                                                      cGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACaa               >  1:82171/1‑150 (MQ=255)
                                                                                                                                      cGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACaa               >  1:82311/1‑150 (MQ=255)
                                                                                                                                       gCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAAt              >  1:2789578/1‑150 (MQ=38)
                                                                                                                                       gCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAAt              <  2:3461550/150‑1 (MQ=255)
                                                                                                                                        cTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATg             >  2:4470124/1‑150 (MQ=255)
                                                                                                                                         tCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATTAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGc            >  1:4594200/1‑150 (MQ=38)
                                                                                                                                           ccACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGCCACACCGCAAATTCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGcc           >  1:6467821/1‑150 (MQ=17)
                                                                                                                                           ccACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTtc                                            >  1:6983721/1‑116 (MQ=25)
                                                                                                                                           ccACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg          >  1:1515634/1‑150 (MQ=38)
                                                                                                                                           ccACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg          >  1:2761456/1‑150 (MQ=38)
                                                                                                                                           ccACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCTAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg          >  1:8129283/1‑150 (MQ=255)
                                                                                                                                           ccACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg          >  2:5361016/1‑150 (MQ=255)
                                                                                                                                            cacaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg                                      <  2:3868087/121‑1 (MQ=255)
                                                                                                                                            cacaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGa         >  2:7587816/1‑150 (MQ=255)
                                                                                                                                            cacaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGa         >  2:1130994/1‑150 (MQ=255)
                                                                                                                                             acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTc                                                               >  2:4204671/1‑95 (MQ=255)
                                                                                                                                             acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAAc                          <  1:8968092/132‑1 (MQ=255)
                                                                                                                                             acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAAc                          >  2:8967924/1‑132 (MQ=255)
                                                                                                                                             acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa        >  2:3023492/1‑150 (MQ=255)
                                                                                                                                             acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa        <  1:8664831/150‑1 (MQ=255)
                                                                                                                                               agaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGc      >  2:5878044/4‑150 (MQ=255)
                                                                                                                                               aGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGc      >  1:6879406/1‑150 (MQ=255)
                                                                                                                                                 aggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGCCACACCGCAAATTCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGcc     >  1:8379544/2‑150 (MQ=18)
                                                                                                                                                  ggTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTc                                                                           >  1:2245115/1‑79 (MQ=255)
                                                                                                                                                   gTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgc  <  1:8960315/150‑1 (MQ=255)
                                                                                                                                                    tttCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgc  <  1:6675127/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GCAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGC  >  NC_002947/5657455‑5657751

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: