Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,655,641 Δ1 bp 16.5% intergenic (‑96/+56) 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,655,6410T.16.5% 184.0 / 41.3 79intergenic (‑96/+56)epd/tktAD‑erythrose 4‑phosphate dehydrogenase/transketolase A
Reads supporting (aligned to +/- strand):  ref base T (24/37);  new base . (7/6);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.74e-01

GGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACC  >  NC_002947/5655494‑5655789
                                                                                                                                                    |                                                                                                                                                    
ggCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACctc                                                                                                                                                     >  2:173252/1‑149 (MQ=255)
 gCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACctct                                                                                                                                                    <  7:73618/149‑1 (MQ=255)
 gCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACctct                                                                                                                                                    <  1:300745/149‑1 (MQ=255)
    gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgtg                                                                                                                                                 <  4:130567/149‑1 (MQ=255)
    gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGAAGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAAATTTCGGAGTGGCATTAGCGGGCAAGCCCGCTCCCACAGAGGTTTCATGAAACCCtgtg                                                                                                                                                 >  5:188093/1‑149 (MQ=37)
     aTACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTgg                                                                                                                                                <  3:369310/149‑1 (MQ=255)
     aTACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTgg                                                                                                                                                <  5:106796/149‑1 (MQ=255)
     aTACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTgg                                                                                                                                                <  5:240868/149‑1 (MQ=255)
      tACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTggg                                                                                                                                               <  1:173252/149‑1 (MQ=255)
        cGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAg                                                                                                                                             <  7:49036/149‑1 (MQ=255)
         ggCCTTAACTGTTTTGTGCAACTTTGTAGGGACTGTGGGGCATTCGGGGCTCACTTGACAGTTCAATTGTAGTCAGCCTGCGGAATTTTCGGTGTTTCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTTCATGAAACCCTGTGGGAGc                                                                                                                                            <  6:188093/149‑1 (MQ=21)
         ggCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGc                                                                                                                                            <  6:54236/149‑1 (MQ=255)
           ccATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCgg                                                                                                                                          <  8:287750/149‑1 (MQ=255)
           ccATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCgg                                                                                                                                          <  8:282594/149‑1 (MQ=255)
           ccATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGAAAGTTCAATTGTCGTCAGCCTTCGGAATTTGAGGTGTTGCATTAGATGGCAATCCCGCTCCAACAGGGTTTCATGCAACCTATGTGGGAGAgg                                                                                                                                          >  5:198463/1‑149 (MQ=255)
                  cGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGGAAGCCCGCTCCCAAAGGGTTTAATGGAAACTGGGTGGGTGAGGGGtttgc                                                                                                                                   >  5:16159/1‑146 (MQ=255)
                      gAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCGACCACAGGGGTTGATGGAACCGGTGTGGGGGCGGGGGTGCCcgc                                                                                                                                >  7:143449/1‑148 (MQ=255)
                       aGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTGCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGa                                                                                                                              <  1:43589/149‑1 (MQ=255)
                            aaCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGaa                                                                                                                         <  4:220368/149‑1 (MQ=255)
                            aaCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCCCGCCCCCAGGGGTTGAAGGAAACACTGTGGGGGCCGGCGTGCCCCCGCAAGaa                                                                                                                         >  3:191538/1‑149 (MQ=255)
                             aCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGa                                                                                                                                              >  1:336344/1‑127 (MQ=37)
                             aCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGa                                                                                                                                              <  2:336344/127‑1 (MQ=37)
                                    aGGGACGGTGGGGCATTCGGGCCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCCTGCAACCGCTGTGGGAGCGGGCTTGGCCGCGAATTAAGGCGccac                                                                                                                 >  2:347397/1‑149 (MQ=38)
                                          ggTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATg                                                                                                                           >  1:110892/1‑133 (MQ=255)
                                          ggTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATg                                                                                                                           <  2:110892/133‑1 (MQ=255)
                                              gggCATTCGGGGCTCACGAGACAGTTCAATTGTCTTCAGCCTGCTGAATTTGCGGTGTGTCATTCGCGTGCAATCCCGCTCCAACATTTTTTCATTCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAg                                                                                                        <  6:198463/148‑1 (MQ=255)
                                                gcctgcGGGCCTCACGAGACAGTTCAATTGGCGTCAGCCTGCGTAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCCCCCACAGGGTTTCATGCAACCTCTGGGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATc                                                                                                     <  1:347397/144‑1 (MQ=38)
                                                gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCACCCACAGGGTTTCATGCAACCGATGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACCCCGCCCGCAGAtg                                                                                                     >  8:223060/1‑148 (MQ=255)
                                                       gggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCa                                                                                                         >  1:243133/1‑138 (MQ=255)
                                                       gggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCa                                                                                                         <  2:243133/138‑1 (MQ=255)
                                                               gagaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTc                                                                                      >  8:351547/1‑149 (MQ=255)
                                                               gagaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTc                                                                                      <  7:223060/149‑1 (MQ=255)
                                                                 gaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCca                                                                                    <  4:155214/149‑1 (MQ=255)
                                                                 gaCAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCca                                                                                    <  7:351547/149‑1 (MQ=255)
                                                                      ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGGAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCagcagc                                                                               >  6:206507/1‑149 (MQ=255)
                                                                      ttCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCagcagc                                                                               >  5:132103/1‑149 (MQ=255)
                                                                        cAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAAAGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTACAGCACct                                                                              >  8:176519/1‑148 (MQ=255)
                                                                             ggcgtcCGCCTGCGGAATTTGCGGGGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCg                                                                         <  5:30781/146‑1 (MQ=255)
                                                                              tcgtcAGCCTTCGTAATTTGCTGTGTGCCATTCGCGGGAAAGCCCGCTCCAAACGGGTTTCATGCAACCTCTTTGGTAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCgg                                                                        <  6:16159/148‑1 (MQ=255)
                                                                              gcggcAGCCTGCGGAATTTGCGGGGTGGCTTTCGGGGGAAACCCCGCTCCCACTGGGTTTCATCCACCCTCTGTGGTTGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCgg                                                                        <  8:143449/144‑1 (MQ=255)
                                                                                gtcAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCGCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAg                                                                     >  3:88568/1‑149 (MQ=255)
                                                                                gtcAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGATCCCACTGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCa                                                                      <  4:191538/148‑1 (MQ=255)
                                                                                gtcAGCCGGCGGAATTTGCGGTGTGGCAGTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATCCAACCTCTGGGGGAGCGGGCTTGCCCGCGAATGAAGGCGaca                                                                                                                  <  1:145389/104‑1 (MQ=37)
                                                                                    gCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAAt                                                                                                                            >  5:354825/1‑90 (MQ=255)
                                                                                    gCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAAt                                                                                                                            <  6:354825/90‑1 (MQ=255)
                                                                                      cTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTAcc                                                                 >  4:302717/1‑147 (MQ=255)
                                                                                       tGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGGTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAg                                                              >  2:111651/1‑149 (MQ=255)
                                                                                        gCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAg                                                              >  8:171253/1‑148 (MQ=255)
                                                                                         cGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGa                                                            <  7:171253/149‑1 (MQ=255)
                                                                                                    gtgtGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAgg                                                 <  7:176519/149‑1 (MQ=255)
                                                                                                     tgtgGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATCCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATGAGt                                                                                          <  3:50648/107‑1 (MQ=255)
                                                                                                       tgGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATATTTAAATAAGTCTTCCAGCAGCTATTAGGCAGTAACGAGGATGTTATCCAGGGt                                               >  2:123518/1‑148 (MQ=37)
                                                                                                         gCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGa                                             <  1:123518/148‑1 (MQ=255)
                                                                                                           aTTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGc                                          >  2:99215/1‑149 (MQ=255)
                                                                                                                         ccGCTCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGaa                             >  1:101793/1‑148 (MQ=37)
                                                                                                                         ccGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGCTCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGaa                             >  5:18335/1‑148 (MQ=255)
                                                                                                                         ccGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAAc                            <  3:302717/149‑1 (MQ=255)
                                                                                                                          cGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGaa                             >  4:104981/1‑147 (MQ=255)
                                                                                                                          cGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGaa                             <  3:104981/147‑1 (MQ=255)
                                                                                                                           gCTCCCACAGGGTTTCGTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACaa                          >  1:54555/1‑149 (MQ=255)
                                                                                                                           gCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACaa                          <  5:214498/149‑1 (MQ=255)
                                                                                                                             tCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaa                                       <  7:269295/134‑1 (MQ=37)
                                                                                                                             tCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaa                                       >  8:269295/1‑134 (MQ=37)
                                                                                                                             tCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCca                                                                                    >  8:246618/1‑89 (MQ=255)
                                                                                                                             tCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCca                                                                                    <  7:246618/89‑1 (MQ=255)
                                                                                                                               ccACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACGGTGGATGTTCTGCAGGGTGAAGCCaaa                                      <  7:247487/133‑1 (MQ=25)
                                                                                                                               ccACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCa                                                   >  5:159553/1‑120 (MQ=37)
                                                                                                                               ccACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCa                                                   <  6:159553/120‑1 (MQ=37)
                                                                                                                               ccACAGAGGTTGCATGAAACCCTGTGGGAGAGGGCTTGCACGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACAGAGGATCTTCTCCAGGGAGAAGCCaaa                                      >  8:247487/1‑133 (MQ=25)
                                                                                                                               ccACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAg                                                                     >  2:142035/1‑102 (MQ=255)
                                                                                                                               ccACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAg                                                                     <  1:142035/102‑1 (MQ=255)
                                                                                                                               ccACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGcc                      <  5:206507/149‑1 (MQ=255)
                                                                                                                                cacagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGcc                      >  4:241529/7‑148 (MQ=38)
                                                                                                                                cacaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATgcc                                                                                                                         <  1:103022/49‑3 (MQ=255)
                                                                                                                                cacaGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGAAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGcc                      >  3:55681/1‑148 (MQ=255)
                                                                                                                                 acagaggTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCg                     >  6:291088/6‑148 (MQ=38)
                                                                                                                                    gggTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGc                 >  3:15445/1‑149 (MQ=255)
                                                                                                                                        tttCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCgg                                                                        <  3:74618/91‑1 (MQ=255)
                                                                                                                                        tttCATGCAACCTATGTGGGAGAGGGCTTGCCCGCCAAAGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTAATCCAGCAGCTCATccg                                                                        >  4:74618/1‑89 (MQ=255)
                                                                                                                                           cATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGgcgc           <  5:291088/148‑1 (MQ=255)
                                                                                                                                            aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGACTGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCtt                                 <  4:256192/126‑1 (MQ=255)
                                                                                                                                            aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCtt                                 >  3:256192/1‑126 (MQ=255)
                                                                                                                                            aTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGcc          >  3:106490/1‑149 (MQ=255)
                                                                                                                                                   ctctGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATg                                                          <  6:295630/94‑1 (MQ=255)
                                                                                                                                                   ctctGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATg                                                          >  5:295630/1‑94 (MQ=255)
                                                                                                                                                    tctGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCAcc  >  3:34428/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGCCGGCGCCGACTCACC  >  NC_002947/5655494‑5655789

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: