Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,657,609 C→A 17.8% 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,6090CA17.8% 412.2 / 82.5 169intergenic (‑91/+61)epd/tktAD‑erythrose 4‑phosphate dehydrogenase/transketolase A
Reads supporting (aligned to +/- strand):  ref base C (70/69);  new base A (17/13);  total (87/82)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.03e-01

Only 100 of 170 total aligned reads displayed.

CAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTC  >  NC_002947/5657463‑5657759
                                                                                                                                                   |                                                                                                                                                      
caGCTGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAAc                                                                                                                                                     <  2:2691711/150‑1 (MQ=255)
caGCTGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAAc                                                                                                                                                     <  2:2691947/150‑1 (MQ=255)
caGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAAc                                                                                                                                                     <  2:2691944/150‑1 (MQ=255)
  gCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAAcc                                                                                                                                                    <  2:689075/150‑1 (MQ=255)
     gCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGACTTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCAGGCAACCTCtg                                                                                                                                                >  2:566014/1‑150 (MQ=255)
      ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAAccctgt                                                                                                                                                <  1:1855217/150‑5 (MQ=38)
      ccGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgt                                                                                                                                               <  2:1868493/150‑1 (MQ=255)
       cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCtgtg                                                                                                                                              <  1:1678223/150‑1 (MQ=38)
       cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgtg                                                                                                                                              <  1:947646/150‑1 (MQ=255)
       cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgtg                                                                                                                                              <  1:1657179/150‑1 (MQ=255)
       cGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgtg                                                                                                                                              <  1:754814/150‑1 (MQ=255)
        gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTgg                                                                                                                                             <  1:89262/150‑1 (MQ=38)
        gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTgg                                                                                                                                             <  2:1455907/150‑1 (MQ=255)
        gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTgg                                                                                                                                             <  1:3744455/150‑1 (MQ=255)
        gATACGGCCATAACCGTTAAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTgg                                                                                                                                             <  1:4190581/150‑1 (MQ=255)
         aTACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCCTGCAACCTCgggggg                                                                                                                                            >  2:769273/1‑150 (MQ=255)
         aTACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTggg                                                                                                                                            <  2:3231175/150‑1 (MQ=255)
         aTACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTggg                                                                                                                                            <  2:3231142/150‑1 (MQ=255)
         aTACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCACTGTggg                                                                                                                                            >  1:205854/1‑150 (MQ=255)
          tACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGa                                                                                                                                           <  1:574928/150‑1 (MQ=255)
          tACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGa                                                                                                                                           <  1:569403/150‑1 (MQ=255)
            cGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGc                                                                                                                                         <  2:3915299/150‑1 (MQ=255)
             ggCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCg                                                                                                                                        <  2:2425137/150‑1 (MQ=38)
             ggCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCg                                                                                                                                        <  2:2425151/150‑1 (MQ=38)
              gCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCgg                                                                                                                                       <  2:867712/150‑1 (MQ=37)
               ccATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCggg                                                                                                                                      <  2:2488067/150‑1 (MQ=255)
                cATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGc                                                                                                                                     <  1:2714137/150‑1 (MQ=255)
                   aaCCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTg                                                                                                                                  <  1:133557/150‑1 (MQ=255)
                    aCCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGc                                                                                                                                 <  2:264261/150‑1 (MQ=255)
                             tGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATg                                                                                                                        >  2:278995/1‑150 (MQ=255)
                             tGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATg                                                                                                                        >  2:2547113/1‑150 (MQ=255)
                              gCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGa                                                                                                                       >  2:4061558/1‑150 (MQ=255)
                                  cTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGTTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGc                                                                                                                   <  2:12383/150‑1 (MQ=255)
                                    ttGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGa                                                                                                                 >  2:3094160/1‑150 (MQ=255)
                                       tAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGacac                                                                                                              >  2:1988277/1‑150 (MQ=255)
                                       tAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCGCCCACAGGGGTTGATGCAACCACTGGGGGAGCGGGCTTGGCCCCGCATGAAGGCGaccc                                                                                                              >  1:392628/1‑148 (MQ=255)
                                         gggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACAccg                                                                                                            >  1:4420219/1‑150 (MQ=255)
                                          ggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCACTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACAccgc                                                                                                           <  2:4303784/150‑1 (MQ=255)
                                           gACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCACTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACAccgcc                                                                                                          >  2:1960254/1‑150 (MQ=255)
                                              ggTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCCCAGGGGTTCATGCAACCTCTGGGGGAGCGGGCTTGCCCGCGACTGGAGGCGACCccgcccgc                                                                                                       >  2:3732977/1‑150 (MQ=38)
                                              ggTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCCCTCCCACACGGTTTCATGCAACCACTGGGGGAGCGGGCTTGCCCGCGAATGAAGGCGACAccgcccgc                                                                                                       >  1:1781448/1‑150 (MQ=255)
                                              ggTGGGGCATTCGGGGCTCACGAGACAGGTCAATTGTCGTCAGCCTGCGGACGTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGGTTCATGCAACCTCGGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACAccgcccgc                                                                                                       >  2:3733081/1‑150 (MQ=38)
                                                 ggggCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGa                                                                                                    >  2:1690777/1‑150 (MQ=255)
                                                    gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCt                                                                                                 >  2:2181538/1‑150 (MQ=255)
                                                    gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCCCCCGCAGATCt                                                                                                 >  1:3679303/1‑150 (MQ=255)
                                                     cATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCTTTCGCGGGCAAGCCCGCTCCCACTGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCtt                                                                                                <  1:3465511/150‑1 (MQ=255)
                                                      aTTCGGGGCTCCCGAGACATTTCACTTGTCGTCAGCGTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCttt                                                                                               <  1:2328334/150‑1 (MQ=255)
                                                          ggggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAAt                                                                                           <  2:1596734/150‑1 (MQ=255)
                                                            ggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCa                                                                                         <  2:4038221/150‑1 (MQ=255)
                                                                 aCGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCtt                                                                                    >  2:4003659/1‑150 (MQ=255)
                                                                 aCGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCCCGCCCCCAGCGGTTCATGCAACCACAGGGGGAGCGGGCGTGCCCCCGCCTGAAGGCGCCCCCCCCCCCCGCGCGTTAAATCAGTCtt                                                                                    >  1:484523/1‑150 (MQ=255)
                                                                  cGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCACTGTGGGAGCGGGCTTGCCCGCGAAAGAAGGCGACACCGCCCGCAGAACTTTAAATCAGTCTTc                                                                                   >  1:339903/1‑150 (MQ=255)
                                                                        aGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCcagcag                                                                             <  1:101369/150‑1 (MQ=255)
                                                                          ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCCCAGGGTTTCATGCAAACACGGTGGGAGCGGGCGTTCCCGCGCATGAAGGCGACACCCCCCCCCCATATTTAAATAAGTCTTACAGCAGct                                                                           >  1:2701903/1‑150 (MQ=37)
                                                                          ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGct                                                                           >  2:293536/1‑150 (MQ=255)
                                                                          ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGct                                                                           >  2:4313836/1‑150 (MQ=255)
                                                                          ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTACAGCAGct                                                                           >  1:1039075/1‑150 (MQ=255)
                                                                          ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCGCGGGGGGAGCGGGCGTGCCCGCGAATGAAGGCGACACCGCCCCCAGAACTTTAAATCAGTCTTCCAGCAGct                                                                           >  2:3028626/1‑150 (MQ=38)
                                                                                tGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCgg                                                                     <  2:1175812/150‑1 (MQ=255)
                                                                                  tcgtcAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCa                                                                   <  1:3300210/150‑1 (MQ=255)
                                                                                  tcgtcAGCCTGCGGAATGTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCa                                                                   <  2:4226962/150‑1 (MQ=255)
                                                                                     tcAGCCTGCGGAATTTGCGGTGTGGCATGCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTa                                                                <  1:3526909/150‑1 (MQ=255)
                                                                                          cTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAg                                                           <  1:289055/150‑1 (MQ=255)
                                                                                          cTGCGGAATTTGCGGGGGGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAg                                                           <  1:1189573/150‑1 (MQ=255)
                                                                                            gCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGa                                                         <  2:4271150/150‑1 (MQ=255)
                                                                                              ggAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAAACAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATg                                                       >  1:3372744/1‑150 (MQ=255)
                                                                                                  tttGCGGTGTGGCATTAGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAAAGAATGCGAAACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTtct                                                   >  2:2010670/1‑150 (MQ=255)
                                                                                                         tgtgGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGt                                            <  2:4485851/150‑1 (MQ=255)
                                                                                                         tgtgGCATTCGCGGGCAAGCCAGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGt                                            <  2:2004305/150‑1 (MQ=255)
                                                                                                            ggCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGaa                                         <  1:801314/150‑1 (MQ=255)
                                                                                                                   gcgGGCAAGCCCGCTCCCACAGGGTTTCAAGCAACCACGGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCCCAGATCTTTAAATCAGTCTTCCAGCAGCTCATCGGCAGGACCGAGGATGTTCTCCAGGGTGAAGCCAAAc                                  >  1:2142977/1‑150 (MQ=38)
                                                                                                                     gggCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCCAActc                                >  1:3149077/1‑150 (MQ=255)
                                                                                                                        cAAGCCCGCTCCCACTGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc                             <  1:3034575/150‑1 (MQ=255)
                                                                                                                        cAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc                             <  1:1690777/150‑1 (MQ=255)
                                                                                                                        cAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTc                             <  1:1920063/150‑1 (MQ=255)
                                                                                                                           gCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGaa                          <  2:3072371/150‑1 (MQ=255)
                                                                                                                            cccGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAAc                         >  1:323938/1‑150 (MQ=255)
                                                                                                                            cccGCTCCCACAGGGTTTCAGGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTACGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTATTCGAAc                         >  1:3800974/1‑150 (MQ=255)
                                                                                                                              cGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACaa                       <  1:4294907/150‑1 (MQ=255)
                                                                                                                               gCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAAt                      >  2:3555955/1‑150 (MQ=255)
                                                                                                                                cTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATg                     <  1:3064886/150‑1 (MQ=255)
                                                                                                                                   ccACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg                  >  1:213948/1‑150 (MQ=38)
                                                                                                                                    cacagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGa                 >  1:3780935/7‑150 (MQ=38)
                                                                                                                                    cacaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGa                 >  2:1069173/1‑150 (MQ=255)
                                                                                                                                    cacaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGa                 <  1:3555955/150‑1 (MQ=255)
                                                                                                                                     acagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa                >  2:511778/6‑150 (MQ=38)
                                                                                                                                     acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa                >  2:1261360/1‑150 (MQ=255)
                                                                                                                                     acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa                >  1:1263157/1‑150 (MQ=255)
                                                                                                                                     acaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGaa                <  2:2590427/150‑1 (MQ=255)
                                                                                                                                       aGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGc              >  1:3960050/1‑150 (MQ=255)
                                                                                                                                            tttCATGCACCCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGGACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgc          <  2:951359/150‑1 (MQ=255)
                                                                                                                                            tttCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgc          >  2:1069817/1‑150 (MQ=255)
                                                                                                                                            tttCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTGTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCCGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgc          >  2:1070002/1‑150 (MQ=255)
                                                                                                                                            tttCAGGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGAACTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgc          <  2:1614972/150‑1 (MQ=255)
                                                                                                                                               cATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGCAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCg      >  1:530984/1‑150 (MQ=255)
                                                                                                                                               cATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAAT‑‑‑‑GCCACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTc  >  1:3855292/1‑150 (MQ=38)
                                                                                                                                                 tGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGa     >  1:1947387/1‑150 (MQ=255)
                                                                                                                                                  gCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGAc    >  2:1666578/1‑150 (MQ=255)
                                                                                                                                                  gCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGAc    >  1:3715106/1‑150 (MQ=255)
                                                                                                                                                  gCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGAc    >  1:3715063/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
CAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTC  >  NC_002947/5657463‑5657759

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: