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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,655,6310TG12.5% 226.9 / 23.3 88intergenic (‑86/+66)epd/tktAD‑erythrose 4‑phosphate dehydrogenase/transketolase A
Reads supporting (aligned to +/- strand):  ref base T (32/45);  new base G (7/4);  total (39/49)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

CAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGC  >  NC_002947/5655483‑5655774
                                                                                                                                                     |                                                                                                                                               
cAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGttt                                                                                                                                                 <  7:216791/149‑1 (MQ=255)
  ggACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCa                                                                                                                                               <  6:301387/149‑1 (MQ=255)
     cacaGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATg                                                                                                                                             <  2:340867/148‑1 (MQ=255)
        aGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCAt                                                                                                                                              <  2:208493/144‑1 (MQ=255)
        aGCGGCCGATACGGCCATAACCGTTGAATGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCAt                                                                                                                                              >  1:208493/1‑144 (MQ=255)
          cGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCa                                                                                                                                           <  1:151865/145‑1 (MQ=255)
           ggCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACctc                                                                                                                                      >  2:173252/1‑149 (MQ=255)
            gCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACctct                                                                                                                                     <  1:300745/149‑1 (MQ=255)
            gCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACctct                                                                                                                                     <  7:73618/149‑1 (MQ=255)
               gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCtgtg                                                                                                                                  <  4:130567/149‑1 (MQ=255)
               gATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGAAGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAAATTTCGGAGTGGCATTAGCGGGCAAGCCCGCTCCCACAGAGGTTTCATGAAACCCtgtg                                                                                                                                  >  5:188093/1‑149 (MQ=37)
                aTACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTgg                                                                                                                                 <  5:240868/149‑1 (MQ=255)
                aTACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTgg                                                                                                                                 <  5:106796/149‑1 (MQ=255)
                aTACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTgg                                                                                                                                 <  3:369310/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)
                     gCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAAc                                                                                                                                         <  2:286928/136‑1 (MQ=255)
                     gCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTAACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAAc                                                                                                                                         >  1:286928/1‑136 (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)
                                   gTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTc                                                                                                                                                >  6:387312/1‑115 (MQ=255)
                                   gTGCAACTTTGTAGGGACGCTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTc                                                                                                                                                <  5:387312/115‑1 (MQ=255)
                                       aaCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGaa                                                                                                          <  4:220368/149‑1 (MQ=255)
                                       aaCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCCCGCCCCCAGGGGTTGAAGGAAACACTGTGGGGGCCGGCGTGCCCCCGCAAGaa                                                                                                          >  3:191538/1‑149 (MQ=255)
                                        aCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGa                                                                                                                               <  2:336344/127‑1 (MQ=37)
                                        aCTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGAGGTTGCATGAAACCCTGTGGGa                                                                                                                               >  1:336344/1‑127 (MQ=37)
                                               aGGGACGGTGGGGCATTCGGGCCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCCTGCAACCGCTGTGGGAGCGGGCTTGGCCGCGAATTAAGGCGccac                                                                                                  >  2:347397/1‑149 (MQ=38)
                                                 ggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTc                                                                                                                                                <  2:301839/101‑1 (MQ=255)
                                                 ggACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTc                                                                                                                                                >  1:301839/1‑101 (MQ=255)
                                                     ggTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATg                                                                                                            <  2:110892/133‑1 (MQ=255)
                                                     ggTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATg                                                                                                            >  1:110892/1‑133 (MQ=255)
                                                         gggCATTCGGGGCTCACGAGACAGTTCAATTGTCTTCAGCCTGCTGAATTTGCGGTGTGTCATTCGCGTGCAATCCCGCTCCAACATTTTTTCATTCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAg                                                                                         <  6:198463/148‑1 (MQ=255)
                                                           gcctgcGGGCCTCACGAGACAGTTCAATTGGCGTCAGCCTGCGTAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCCCCCACAGGGTTTCATGCAACCTCTGGGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATc                                                                                      <  1:347397/144‑1 (MQ=38)
                                                           gCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCACCCACAGGGTTTCATGCAACCGATGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACCCCGCCCGCAGAtg                                                                                      >  8:223060/1‑148 (MQ=255)
                                                                  gggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCa                                                                                          <  2:243133/138‑1 (MQ=255)
                                                                  gggCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCa                                                                                          >  1:243133/1‑138 (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                                                                                                             <  6:354825/90‑1 (MQ=255)
                                                                                               gCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAAt                                                                                                             >  5:354825/1‑90 (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              <  3:104981/147‑1 (MQ=255)
                                                                                                                                     cGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGaa              >  4:104981/1‑147 (MQ=255)
                                                                                                                                      gCTCCCACAGGGTTTCGTGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACaa           >  1:54555/1‑149 (MQ=255)
                                                                                                                                      gCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACaa           <  5:214498/149‑1 (MQ=255)
                                                                                                                                        tCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaa                        >  8:269295/1‑134 (MQ=37)
                                                                                                                                        tCCCACAGAGGTTGCATGAAACCCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCaa                        <  7:269295/134‑1 (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                                                      <  1:142035/102‑1 (MQ=255)
                                                                                                                                          ccACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAg                                                      >  2:142035/1‑102 (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)
                                                                                                                                                     |                                                                                                                                               
CAGGACACAGCGGCCGATACGGCCATAACCGTTGAGTGCAACTTTGTAGGGACGGTGGGGCATTCGGGGCTCACGAGACAGTTCAATTGTCGTCAGCCTGCGGAATTTGCGGTGTGGCATTCGCGGGCAAGCCCGCTCCCACAGGGTTTCATGCAACCTCTGTGGGAGCGGGCTTGCCCGCGAATGAAGGCGACACCGCCCGCAGATCTTTAAATCAGTCTTCCAGCAGCTCTTCGGCAGTACCGAGGATGTTCTCCAGGGTGAAGCCAAACTCTTCGAACAATGCCGAAGC  >  NC_002947/5655483‑5655774

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: