Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,657,609 C→A 25.1% intergenic (‑91/+61) 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,6090CA25.1% 296.3 / 108.8 159intergenic (‑91/+61)epd/tktAD‑erythrose 4‑phosphate dehydrogenase/transketolase A
Reads supporting (aligned to +/- strand):  ref base C (63/56);  new base A (19/21);  total (82/77)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

Only 100 of 161 total aligned reads displayed.

ACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGA  >  NC_002947/5657460‑5657756
                                                                                                                                                      |                                                                                                                                                   
  acaGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCaa                                                                                                                                                   >  2:2425304/1‑150 (MQ=255)
  acaGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCaa                                                                                                                                                   >  2:592380/1‑150 (MQ=255)
   caGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAAc                                                                                                                                                  <  1:1367977/150‑1 (MQ=255)
    aGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAAcc                                                                                                                                                 >  2:4876244/1‑150 (MQ=255)
     gCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACct                                                                                                                                                <  2:4136915/150‑1 (MQ=255)
       ggCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAAccct                                                                                                                                               <  1:4423924/150‑3 (MQ=255)
       ggCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAAccct                                                                                                                                               <  1:3605569/150‑3 (MQ=255)
         ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                            <  2:4757009/150‑1 (MQ=255)
         ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                            <  1:1941474/150‑1 (MQ=255)
         ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                            <  1:4644787/150‑1 (MQ=255)
         ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                            <  2:2594357/150‑1 (MQ=255)
         ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAACCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAAccctgt                                                                                                                                             <  2:4653024/150‑5 (MQ=38)
          cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCtgtg                                                                                                                                           <  2:3428728/150‑1 (MQ=38)
          cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCtgtg                                                                                                                                           <  2:3344155/150‑1 (MQ=38)
          cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCtgtg                                                                                                                                           <  1:4025181/150‑1 (MQ=38)
          cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgtg                                                                                                                                           >  2:3247150/1‑150 (MQ=255)
           gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTgg                                                                                                                                          >  1:1256680/1‑150 (MQ=255)
           gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTgg                                                                                                                                          <  1:5298589/150‑1 (MQ=255)
             tACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGa                                                                                                                                        <  2:3724537/150‑1 (MQ=38)
             tACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGa                                                                                                                                        <  2:2715680/150‑1 (MQ=38)
             tACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGa                                                                                                                                        <  2:3164433/150‑1 (MQ=255)
              aCGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCTGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAg                                                                                                                                       <  2:3684658/150‑1 (MQ=255)
               cGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGc                                                                                                                                      <  1:1949648/150‑1 (MQ=255)
                 gCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCgg                                                                                                                                    <  1:1078591/150‑1 (MQ=37)
                   cATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAATCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGc                                                                                                                                  >  2:195902/1‑150 (MQ=255)
                   cATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGc                                                                                                                                  <  2:4216552/150‑1 (MQ=255)
                   cATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAACCCCGCGCCCACAGGGGTTAATGCAACCACTGGGGGGGCCGGc                                                                                                                                  >  1:4096265/1‑150 (MQ=255)
                    aTAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                            <  1:1782786/139‑1 (MQ=255)
                      aaCCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTg                                                                                                                               >  2:2214519/1‑150 (MQ=255)
                      aaCCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCGTg                                                                                                                               >  2:638216/1‑150 (MQ=255)
                       aCCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGc                                                                                                                              <  1:2947107/150‑1 (MQ=255)
                                tGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGa                                                                                                                                        <  1:1280917/131‑1 (MQ=37)
                                tGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGa                                                                                                                                        >  2:1280917/1‑131 (MQ=37)
                                 gCAACTTTGTAGGGAGGGTGGGGCATTCGGGGCTCACGAGACAGTACATTTGTCGTCAGCCTGCGGAATTTGCGGGGGGGCATTCGCGGGCAAGCCCGCTCCCACCGGGTTTAATGCAACCTCTGGGGGAGCGGGCTTGCCCGCGAATGa                                                                                                                    <  1:1520276/150‑1 (MQ=38)
                                    aCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAgg                                                                                                                 <  2:2001541/150‑1 (MQ=255)
                                      tttGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCg                                                                                                               >  1:824311/1‑150 (MQ=255)
                                          tAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGGGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGGGGGAGCGGGCTTGCCCGCGAATGAAGGCGacac                                                                                                           >  1:3541912/1‑150 (MQ=255)
                                            gggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACAccg                                                                                                         >  1:4088337/1‑150 (MQ=255)
                                                       gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCt                                                                                              >  2:399073/1‑150 (MQ=255)
                                                         aTTCGGGGCTCACGAGACTGTTCATTTGTGGTCAGCCTGCGGATTTTGCGGTGGGGCATTCGCGGGCACGCCCGCTCCCACTGGGTTTCATGCAACCTCGGTGGGCGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCttt                                                                                            <  1:638216/150‑1 (MQ=255)
                                                           tCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGAAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTaa                                                                                          <  1:4822922/150‑1 (MQ=255)
                                                             ggggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGc                                                                                                                                      <  2:2065131/104‑1 (MQ=21)
                                                             ggggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGc                                                                                                                                      >  1:2065131/1‑104 (MQ=21)
                                                             ggggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGc                                                                                                                                      >  1:2064955/1‑104 (MQ=21)
                                                             ggggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAAt                                                                                        <  1:399073/150‑1 (MQ=255)
                                                                  tCACGAGACAGTCCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTc                                                                                   >  2:2006834/1‑150 (MQ=255)
                                                                      gagaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCg                                                                                                                                     <  2:531952/96‑1 (MQ=255)
                                                                             ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCACAGTGGGAGCGGGCATGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCagcagcg                                                                        >  2:120330/1‑149 (MQ=255)
                                                                                    gtcgtcAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGc                                                                 <  2:321343/150‑1 (MQ=255)
                                                                                     tcgtcAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCa                                                                >  2:2561709/1‑150 (MQ=255)
                                                                                        tcAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCCCTCCCACAGGGTTTCAAGCAACCGCTGTGGGAGCGGGCGTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCATTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTa                                                             >  2:2337106/1‑150 (MQ=255)
                                                                                           gCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCg                                                          <  2:3739013/150‑1 (MQ=255)
                                                                                            ccTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGa                                                         <  1:3247150/150‑1 (MQ=255)
                                                                                               gCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCGTCCAGCAGCTCTTCGGCAGTACCGAGGa                                                      >  2:2743883/1‑150 (MQ=255)
                                                                                               gCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCCCCCGCAGATCTTTAAATCAGTCTTACAGCAGCTCTTCGGCAGTACCGAGGa                                                      >  2:1921176/1‑150 (MQ=255)
                                                                                                 ggAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATg                                                    <  1:592380/150‑1 (MQ=255)
                                                                                                           gtgtGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAggg                                          <  2:4341124/150‑1 (MQ=255)
                                                                                                            tgtgGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGt                                         <  1:1468096/150‑1 (MQ=255)
                                                                                                            tgtgGCAGTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGt                                         <  2:1176508/150‑1 (MQ=255)
                                                                                                               ggCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGaa                                      <  1:3746346/150‑1 (MQ=255)
                                                                                                               ggCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGaa                                      <  2:2880400/150‑1 (MQ=255)
                                                                                                                  aTTCGCGGGCAAGCCCGCTCCCACCGGGTTTCATGCCACCGCGGGGGGCGCGGGCGTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTGTAAATCAGTCTTCCAGCAGCTCGTCGGCAGTACCGAGGATGGTCTCCATGGTGAAGcc                                   >  2:169939/1‑150 (MQ=255)
                                                                                                                     cgcgGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAggg                                          >  1:5364285/1‑140 (MQ=255)
                                                                                                                        gggCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAActc                             >  1:2886945/1‑150 (MQ=255)
                                                                                                                          gCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCGTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCtt                           <  2:4358480/150‑1 (MQ=255)
                                                                                                                           ccaGCCCGCCCCCCCCGGGGGTCACGCCCCCCCCGGGGGAGCGGTCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGACCTTTAAATCAGTCTGCCAGCAGCTCGGCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc                          <  2:5425882/148‑1 (MQ=255)
                                                                                                                           cAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc                          <  2:3112746/150‑1 (MQ=255)
                                                                                                                            aaGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCg                         >  2:4143485/1‑150 (MQ=255)
                                                                                                                            aaGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCg                         >  1:4878484/1‑150 (MQ=255)
                                                                                                                              gCCCGCTCCCACAGGGTTTCATGCAACCTCGGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCCGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGaa                       >  2:4894358/1‑150 (MQ=255)
                                                                                                                               cccGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGa                                                         >  1:1292031/1‑115 (MQ=255)
                                                                                                                               cccGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGa                                                         <  2:1292031/115‑1 (MQ=255)
                                                                                                                               cccGCTCCCACAGGGTTTCATGCAACCACTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAAc                      >  2:1730576/1‑150 (MQ=255)
                                                                                                                                ccGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACa                     >  1:4182915/1‑150 (MQ=37)
                                                                                                                                 cGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACaa                    >  1:1515578/1‑150 (MQ=255)
                                                                                                                                   cTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATg                  >  1:2729594/1‑150 (MQ=255)
                                                                                                                                    tCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCca                                                                              >  2:4652552/1‑89 (MQ=18)
                                                                                                                                    tCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCca                                                                              <  1:4652552/89‑1 (MQ=18)
                                                                                                                                    tCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGc                 >  1:5076410/1‑150 (MQ=38)
                                                                                                                                      ccACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg               >  2:239220/1‑150 (MQ=38)
                                                                                                                                      ccACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg               >  1:4577635/1‑150 (MQ=38)
                                                                                                                                      ccACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg               >  2:1793257/1‑150 (MQ=255)
                                                                                                                                      ccACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg               >  2:2464521/1‑150 (MQ=255)
                                                                                                                                       cacagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGa              >  1:2083316/7‑150 (MQ=38)
                                                                                                                                       cacaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGa              >  1:1005020/1‑150 (MQ=255)
                                                                                                                                        acagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa             >  2:512318/6‑150 (MQ=38)
                                                                                                                                        acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa             >  2:2245514/1‑150 (MQ=255)
                                                                                                                                        acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa             >  1:2647851/1‑150 (MQ=255)
                                                                                                                                        acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa             <  2:2886945/150‑1 (MQ=255)
                                                                                                                                        acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa             <  2:2470550/150‑1 (MQ=255)
                                                                                                                                           gggTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGcc          <  2:1256680/150‑1 (MQ=255)
                                                                                                                                           gggTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGcc          >  2:3296952/1‑150 (MQ=255)
                                                                                                                                              gTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCgg        <  2:5133761/150‑1 (MQ=255)
                                                                                                                                                 tCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgcgc     <  2:4096265/150‑1 (MQ=255)
                                                                                                                                                 tCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgcgc     >  2:3775689/1‑150 (MQ=255)
                                                                                                                                                 tCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgcgc     >  1:494137/1‑150 (MQ=255)
                                                                                                                                                   aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCg   >  2:3582534/1‑150 (MQ=255)
                                                                                                                                                   aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCg   >  2:3506168/1‑150 (MQ=255)
                                                                                                                                                    tGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGa  <  2:4088337/150‑1 (MQ=255)
                                                                                                                                                    tGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGa  >  2:236413/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
ACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGA  >  NC_002947/5657460‑5657756

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: