Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*U00096422,4290TG20.7% 440.4 / 25.0 191intergenic (+70/‑89)proY/malZputative transporter ProY/maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base T (23/126);  new base G (39/0);  total (65/126)
Fisher's exact test for biased strand distribution p-value = 1.72e-24
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 200 total aligned reads displayed.

TCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTG  >  U00096/422283‑422577
                                                                                                                                                  |                                                                                                                                                    
tCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcc                                                                                                                                                    <  1:408472/149‑1 (MQ=255)
  gcgTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGgccgg                                                                                                                                                  >  2:408615/1‑147 (MQ=255)
  gcgTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctc                                                                                                                                                  >  2:408548/1‑149 (MQ=255)
  gcgTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctc                                                                                                                                                  >  1:408554/1‑149 (MQ=255)
             tgggCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGGCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAat                                                                                                                                       <  2:408475/146‑1 (MQ=255)
              gtgCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAata                                                                                                                                      >  2:408579/1‑149 (MQ=255)
              gtgCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCGCCCCCAAata                                                                                                                                      >  2:4180875/1‑149 (MQ=255)
                gCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATc                                                                                                                                    <  1:408476/149‑1 (MQ=255)
                     tGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCa                                                                                                                               <  1:408478/149‑1 (MQ=255)
                     tGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATCCTGCa                                                                                                                               >  2:408569/1‑149 (MQ=255)
                        ttGGCTGGATGTTTAAACGCCGCCACGAGCGTCAGCTGGCGGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAga                                                                                                                             <  1:408481/148‑1 (MQ=255)
                           gCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAgatgat                                                                                                                         <  1:408483/149‑1 (MQ=255)
                               gATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGt                                                                                                                     <  2:408485/149‑1 (MQ=255)
                               gATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGt                                                                                                                     <  2:408487/149‑1 (MQ=255)
                                aTGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGCGCCGGCCGCCGCCCCCCAATATCCTTCAggtgtggggtg                                                                                                                    >  1:408567/1‑142 (MQ=255)
                                  gTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGAtgatga                                                                                                                       <  1:408488/144‑1 (MQ=255)
                                   tttAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                                                                                                  <  1:408490/148‑1 (MQ=255)
                                     tAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCt                                                                                                               <  1:408493/149‑1 (MQ=255)
                                     tAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCt                                                                                                               <  2:408492/149‑1 (MQ=255)
                                     tAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCt                                                                                                               <  2:408496/149‑1 (MQ=255)
                                      aaaCGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCGCCGGCCGCCGCCCCCCAAAAGACAGCAGAGGATGAGAGAGccgg                                                                                                              >  2:408636/1‑147 (MQ=255)
                                       aaCGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGc                                                                                                             <  1:408497/149‑1 (MQ=255)
                                          gccgccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCAtt                                                                                                          <  1:408500/149‑1 (MQ=255)
                                            cgccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAg                                                                                                        <  1:408502/149‑1 (MQ=255)
                                              ccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGc                                                                                                      <  1:408505/149‑1 (MQ=255)
                                              ccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCAAATATCCTTCAGATGATGAGTGATGCTGCATTAGGc                                                                                                      >  2:408621/1‑149 (MQ=255)
                                              ccACGATCGTCAGCCGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTGCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGc                                                                                                      <  2:408503/149‑1 (MQ=255)
                                                 ctaTCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTAt                                                                                                   <  1:408507/147‑1 (MQ=255)
                                                  gATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCGCCCCCAAATAGCCTTCAGATGATGAGTGATCCTGCATTAGGGTATg                                                                                                  >  2:408595/1‑149 (MQ=255)
                                                     cGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCa                                                                                               <  1:408510/149‑1 (MQ=255)
                                                         aGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGt                                                                                           <  1:408512/149‑1 (MQ=255)
                                                          gctggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTg                                                                                          <  2:408513/149‑1 (MQ=255)
                                                           ctggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGa                                                                                         <  1:408516/149‑1 (MQ=255)
                                                            tggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGAt                                                                                        <  2:408518/149‑1 (MQ=255)
                                                              gctgAAAACCAGTATTCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCa                                                                                      <  1:408519/149‑1 (MQ=255)
                                                              gctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGAt                                                                                        <  1:408521/147‑1 (MQ=255)
                                                               ctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCGCCCCCAAATATCCTTCAGATGATGAGTGATCCGGGAGTAGGGTATGGCAAGGTGATGAg                                                                                     >  2:408677/1‑149 (MQ=255)
                                                                tgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGa                                                                                    <  1:408523/149‑1 (MQ=255)
                                                                  aaaaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAtt                                                                                  <  1:408526/149‑1 (MQ=255)
                                                                  aaaaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAtt                                                                                  >  2:408589/1‑149 (MQ=255)
                                                                  aaaaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATCCTTCAGATGATGAGTGATGCTGCAGTAGGCTATGGGAAGGGGAGCAGAtt                                                                                  >  2:408659/1‑149 (MQ=255)
                                                                        aGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCt                                                                                                                                  <  1:408528/95‑1 (MQ=255)
                                                                        aGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtcatc                                                                            <  1:408530/149‑1 (MQ=255)
                                                                          tAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATcac                                                                          <  2:408532/149‑1 (MQ=255)
                                                                           aaTCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGTCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATcaca                                                                         <  2:408534/149‑1 (MQ=255)
                                                                           aaTCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCa                                                                                                            <  1:408535/114‑1 (MQ=255)
                                                                             tCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCCAATATCGTTCAGATGATGAGTGATCATGCATTAGGGTATGGCAAGGTGATCAGATGTTCATCAAAgg                                                                       >  2:408603/1‑149 (MQ=255)
                                                                              ccTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGGCCGGCCGCCGCCCCCACAGATCCGGCCGAGGATGAGGGATCCTGCGTGAGGGGAGGGGGAGGTGGTCAGATGTTAATCACAggg                                                                      >  2:408690/1‑149 (MQ=255)
                                                                              ccTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAggg                                                                      <  2:408538/149‑1 (MQ=255)
                                                                               cTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGGGCGGCGGCGGCCCCCAAAAAGCCGGCAGATGATGAGTGATGCTGCATGAGGCGGTGGGAAGGTGATGAGATTTTCATCACAgggg                                                                     >  2:408644/1‑149 (MQ=255)
                                                                                ttGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAggg                                                                      <  1:408539/147‑1 (MQ=255)
                                                                                  gCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAAt                                                                  <  2:408541/149‑1 (MQ=255)
                                                                                    gTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCGCGGGCGCCTCGCCCCAAAAGGCAGGAGCTGAGGAGGGCTCCTGCATTAGGCTATGGCTAGGGGTGCTGAGTGTCATCAAAGGGGGATTa                                                                >  2:408613/1‑149 (MQ=255)
                                                                                    gTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCCAATATACTGCAGAGGATGAGGGAGCCGGGATTAGGGGATGGGAAGGTGATCAGATTTTGATCACAgggggttta                                                                >  2:408640/1‑144 (MQ=255)
                                                                                       cTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGCGCCGGCGGCCGCCCCCAAATATCGGTCAGATGATGAGTGATCCTGCAGTAGGGGATGGCGAGGTGAGGAGATTTTGATGACAGGGGAATTatga                                                             >  2:408642/1‑149 (MQ=255)
                                                                                       cTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTatga                                                             <  2:408543/149‑1 (MQ=255)
                                                                                           aaGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGGGCGGCCTCCTCCCCCAAATATCCTTCAGAGGATGAGGGATCCTGAATTAGGCTATGGCAAGGTGATCAGATTTTAATAACAGGGGAATTATGATGtt                                                         >  2:408686/1‑149 (MQ=255)
                                                                                           aaGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGtt                                                         <  1:408546/149‑1 (MQ=255)
                                                                                             gCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCCAAAATGCTTGAGGTGATGAGTGGTGCTGCAGTAGGGGATGGCGAGGTGGTCAGGTTTTGATCACAGGGGGGTTTTGATGtttt                                                       >  2:408627/1‑147 (MQ=255)
                                                                                              ccGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTaaa                                                      <  1:408548/149‑1 (MQ=255)
                                                                                               cGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAAt                                                     <  2:408550/149‑1 (MQ=255)
                                                                                                 gTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGc                                                   <  2:408552/149‑1 (MQ=255)
                                                                                                   tCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGGGCGGGCGCCGCCCCCAACGATCCGTGAGAGGAGGAGTGAGGCGGCATTAGGCTATGGCAAGGTGGTCAGATGTTCATCACAGGGGAATTAGGAAGTTTAATGAAt                                                 >  2:408668/1‑149 (MQ=255)
                                                                                                   tCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCGGGGGCCGCCGCCCCCAAATATCCGAAAGATGATGAGTGAAGATGCATTAGGCTATGGCAAGGTGATCAGATATTAATAAAAGGGGAATTATGATGTTAATTGAAt                                                 >  2:408652/1‑149 (MQ=255)
                                                                                                   tCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCAt                                                 <  2:408555/149‑1 (MQ=255)
                                                                                                   tCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCTCCCCCAAAGATCCTTCAGAGGATGAGTGATCCGGCAGTAGGCTAGGGCAAGGTGAGAAGATTTTCATCACAGGGGAATTAGGAGGTTAAATGcct                                                 >  2:408665/1‑147 (MQ=255)
                                                                                                    cAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATg                                                <  2:408558/149‑1 (MQ=255)
                                                                                                      ggaggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGc                                              <  1:408562/149‑1 (MQ=255)
                                                                                                      ggaggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGc                                              <  2:408560/149‑1 (MQ=255)
                                                                                                         ggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCGCCGCCCCCAAATATACTTCAGATGATGAGTGGTGGTGGGGTAGGCTATGGCAAGGTGAGCAGATTTTGATCACAGGGGAATTATGATGTTTAATGAATGGCAcc                                           >  2:408759/1‑149 (MQ=255)
                                                                                                            aaTCCTGATCCGGCTATTCCGAAAGTTATCCGGCCCCGGCCGCCGCCCCCACAAATACAACAGATGATGCGTGATCCGGCATTAGGGTATGGGAAGGGGGAGAGAATATGATCACAGGGGAATTATGATGTTAAATGCGTGGCTCCTGc                                        >  2:408662/1‑149 (MQ=255)
                                                                                                            aaTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGc                                        <  2:408563/149‑1 (MQ=255)
                                                                                                             aTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGcc                                       <  1:408564/149‑1 (MQ=255)
                                                                                                               ccTGATCCGGCTATTCCGAAAGTTATCCGGGCCCGGCCTCCGCCCCCAAATATCCGTAACATGAAGGGGGAGCCGGCGATTGGCGATGGGGAGGGGGTCAGAGTTTAATAACAGGGGAATTATGAGGTTTATTGAATGGCACCTGCAgg                                     >  2:408697/1‑149 (MQ=255)
                                                                                                               ccTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCTCCGCCCCCCAATATCCGTCAGATGATGAGTGGTGGTGCCTTTGGGGATGGGAAGGGGGGGAGATTTTCATGACAGGGGAAGTATGATGTTTAATGCATGGCACCTGCCgg                                     >  2:408780/1‑149 (MQ=255)
                                                                                                               ccTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCgg                                     <  2:408567/149‑1 (MQ=255)
                                                                                                                cTGATCCGGCTATTCCGAAAGTTATCCGGGGGCGGGCGCGGGGCCCCGGGAGCCGGCAGATGATGAGTGATCCGGCAGGAGGCGAGGGCAAGGTGAGGAGATTTTGATAAGAGGGGAAGTATGATGTTAGATGAATTTCACGTGCCGGt                                    >  2:408655/1‑149 (MQ=255)
                                                                                                                aggcaccGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGt                                    <  1:2092933/144‑1 (MQ=255)
                                                                                                                         cTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGcc                                       <  1:408570/137‑1 (MQ=255)
                                                                                                                          tATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCAttt                          <  1:408572/149‑1 (MQ=255)
                                                                                                                            ttCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGt                        <  2:408573/149‑1 (MQ=255)
                                                                                                                                gAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTaaa                    <  1:408575/149‑1 (MQ=255)
                                                                                                                                  aaGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAca                  <  1:408578/149‑1 (MQ=255)
                                                                                                                                    gTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaa                <  1:408580/149‑1 (MQ=255)
                                                                                                                                      tATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagc              <  1:408582/149‑1 (MQ=255)
                                                                                                                                      tATCCGCCCCCGTCCGCCGCCCCCCAAGATCCTTCAGATGATGAGTGATGCTGCGTTAGGGTATGGCAAGGGGATCAGATTTTGATCAGAGGGGAAGTATGATGTTTAATGGATGGCCCCTGCCGGGGGGCCCAGGTGGTAAAcaaagc              >  2:408710/1‑149 (MQ=255)
                                                                                                                                       aTCCGCGCCCGGCCGCCGCCCCCAAATATCCGTCAGATGATGAGTGATCCTGCATTAGGGTATGGCAAGGGGATCAGATTTTAATCACAGGGGAATTATGATGTTTAATGCATGGCACCTGCCGGTGCCCCCAGGTGTTAAAcaaagca             >  2:408732/1‑149 (MQ=255)
                                                                                                                                        tCCGTCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagca             <  1:408584/148‑1 (MQ=255)
                                                                                                                                         ccGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagcaaa           <  2:408587/149‑1 (MQ=255)
                                                                                                                                         ccGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagcaa            <  1:408586/148‑1 (MQ=255)
                                                                                                                                          cGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagcaaag          <  1:408589/149‑1 (MQ=255)
                                                                                                                                              cccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCa      <  1:408591/149‑1 (MQ=255)
                                                                                                                                               ccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCaa     >  1:408633/1‑149 (MQ=255)
                                                                                                                                               ccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCaa     <  2:408593/149‑1 (MQ=255)
                                                                                                                                                cGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTTTGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAAc    <  1:408594/149‑1 (MQ=255)
                                                                                                                                                cGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAAc    >  2:2794129/1‑149 (MQ=255)
                                                                                                                                                cGTCCTCCTCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGAt        >  1:408731/1‑144 (MQ=255)
                                                                                                                                                 gTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagc              <  1:408597/138‑1 (MQ=255)
                                                                                                                                                 gTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACt   >  2:408706/1‑149 (MQ=255)
                                                                                                                                                  tcctcctccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTg  >  2:408802/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTG  >  U00096/422283‑422577

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: