Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913422,4290TG29.9% 587.3 / 107.1 354intergenic (+70/‑89)proY/malZproline‑specific permease/maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base T (58/188);  new base G (103/2);  total (164/190)
Fisher's exact test for biased strand distribution p-value = 1.25e-43
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 375 total aligned reads displayed.

TTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCT  >  NC_000913/422281‑422579
                                                                                                                                                    |                                                                                                                                                      
  tCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCGCCGgcctc                                                                                                                                                    >  8:529231/1‑151 (MQ=255)
      gTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGGCGCGtgctccggc                                                                                                                                                >  2:1009988/1‑148 (MQ=255)
       tGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGACGCCACGATCGTCAGCTGGCTCAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATACGGCTATTCCGAAAGTTATCCGCGCGCCGCCGCCTccc                                                                                                                                               >  5:1213607/1‑151 (MQ=255)
          aTTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCACGAGGAAATCCTGATCCGGCTAATCCGAAAGTTATCCGCCCCCGTCCGCCGcaccc                                                                                                                                             >  6:392795/1‑150 (MQ=255)
          aTTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTCTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCAGTCCTCCAccccc                                                                                                                                             >  6:517131/1‑150 (MQ=255)
           ttgttgTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCGCCCCCac                                                                                                                                           >  5:886445/1‑150 (MQ=255)
               tgtgCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAatat                                                                                                                                       >  1:1291994/1‑151 (MQ=255)
                gtgCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATc                                                                                                                                      <  6:1551841/151‑1 (MQ=255)
                gtgCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCCACTATc                                                                                                                                      >  4:1414042/1‑151 (MQ=255)
                gtgCTGTTGATTGGCTGGATGTTTAAACGCCGCAACGATCGTCAGATGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAatat                                                                                                                                       <  2:1561266/150‑1 (MQ=255)
                       tGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGGCCCCCAATATCCTTc                                                                                                                                  >  5:1520027/1‑148 (MQ=255)
                        gATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAgat                                                                                                                              <  3:1083913/151‑1 (MQ=255)
                          ttGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCGCCGGCCGCCTCCCCCAAAGATCCTTCAgatg                                                                                                                             >  7:1367738/1‑150 (MQ=255)
                            ggCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCGCCCCCAAAAATCCTTCAgatgagg                                                                                                                          >  3:878128/1‑149 (MQ=255)
                               tGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCCAATATCCTTAAGATGATGAGt                                                                                                                       >  8:133133/1‑151 (MQ=255)
                               tGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCCAATATCCTTCAGCTGATgtgt                                                                                                                       >  2:755180/1‑151 (MQ=255)
                                ggATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGt                                                                                                                       <  7:676900/150‑1 (MQ=255)
                                 gATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGa                                                                                                                     <  1:783891/151‑1 (MQ=255)
                                 gATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGa                                                                                                                     <  3:1222554/151‑1 (MQ=255)
                                  aTGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                                                                                                    <  6:1054175/151‑1 (MQ=255)
                                  aTGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                                                                                                    <  6:1410076/151‑1 (MQ=255)
                                    gTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATc                                                                                                                   <  7:1097063/150‑1 (MQ=255)
                                    gTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGGCGCCGCCCCCCAATATCCTTCAGAGGATGAGTGCTcg                                                                                                                  >  2:668223/1‑150 (MQ=255)
                                      ttAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCGCCGCCCCCACATAACCTTCAgatgcggcgggatgctg                                                                                                                >  8:1039361/1‑138 (MQ=255)
                                       tAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTg                                                                                                                <  1:463551/150‑1 (MQ=255)
                                       tAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTg                                                                                                                <  1:156185/150‑1 (MQ=255)
                                             ccgccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                                                                                                         <  7:1253206/151‑1 (MQ=255)
                                             ccgccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCACAAATCCTTCAGATGAGGCGTGATCCTGCAGTAgg                                                                                                         >  4:386148/1‑151 (MQ=255)
                                               gccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCAAATATCCTTCAGAGGGTGAGTGCTCCTGCATTAGGcg                                                                                                       >  5:143614/1‑150 (MQ=255)
                                                  aCGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATg                                                                                                    <  1:1687967/151‑1 (MQ=255)
                                                    gATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGTCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCACAACTCTACTTCAGATGAGGTGTGATGCGGCAGTAGGCTAGGGc                                                                                                  >  4:995430/1‑151 (MQ=255)
                                                     aTCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGCCGCGGCCGCCTGCACCAAATATAAGTAAGAGGATGATTGATGATGCAGTAGGCTAttgt                                                                                                  >  3:694558/1‑147 (MQ=255)
                                                        gTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCTCCTCGCCCAAAAATCCTTCAGATGATGAGTGATGCTGCATTAGGCTATGGCAAgg                                                                                              >  1:609809/1‑151 (MQ=255)
                                                          cAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGt                                                                                             <  6:148353/150‑1 (MQ=255)
                                                           aGCTGCCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGa                                                                                           <  8:459541/151‑1 (MQ=255)
                                                            gcttgctgAAAACCTGTATTGCTTGCGTTCTGTAATCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTGCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGAt                                                                                          <  7:323028/151‑1 (MQ=255)
                                                            gctggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGAt                                                                                          <  2:427901/151‑1 (MQ=255)
                                                             ctggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAAAATACTTCAGATGATGAGTGATCCTGCATTAGGGGATGGCAAGGTGATc                                                                                         >  8:105017/1‑151 (MQ=255)
                                                               ggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCCAAGATCCTTCAGATGATGAGTGAGCCTGCATTAGTCTATGGCAAGGTGGTc                                                                                         >  6:641000/1‑149 (MQ=255)
                                                                gctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGa                                                                                      <  8:873570/151‑1 (MQ=255)
                                                                  tgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCGTCCGCCCCCAAATATCCTTAAGATGATGAGGGATCCTGCATTAGGCTATGGCAAGGTTAGCAGAtt                                                                                    >  8:642533/1‑151 (MQ=255)
                                                                    aaaaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGACTCCGGGCGCGAAAGTCCTTGAGATGTTGAGTGATCCGGGCTTAGGGTCTTGAAAGGAGTTCAGAtttt                                                                                  >  2:997747/1‑151 (MQ=255)
                                                                    aaaaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAtttt                                                                                  <  3:93029/151‑1 (MQ=255)
                                                                    aaaaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAtttt                                                                                  <  4:1271177/151‑1 (MQ=255)
                                                                    aaaaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCCCCCCCGTCCTCCGCCCCCAAAACTCCAGAAGATGAGGATTTATCATGCCTTAGGGTATTTCTAGGTGATCAGAtttt                                                                                  >  3:1393699/1‑151 (MQ=255)
                                                                     aaaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATACTTCAGAGGATGAGTGATCCTGCATTAGGGGATGGCAAGGGGATCAGAttttg                                                                                 >  4:590149/1‑150 (MQ=255)
                                                                      aaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCAAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtca                                                                                <  6:301278/151‑1 (MQ=255)
                                                                       aCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtcat                                                                               <  3:1527080/151‑1 (MQ=255)
                                                                       aCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATACTTCAGATGATGAGTGATCCTGCATTAGGGGATGGCGAGGTGATCAGATTTtca                                                                                >  7:940982/1‑150 (MQ=255)
                                                                        ccAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtcatc                                                                              <  5:977149/151‑1 (MQ=255)
                                                                        ccAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtcatc                                                                              <  5:780583/151‑1 (MQ=255)
                                                                         cAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtcatca                                                                             <  2:476478/151‑1 (MQ=255)
                                                                         cAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTAATCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCCAATATCCTGCCGATGATGAGTGCTCCGGCATTAGGGTATGGCGCGGTGCTCAGATTTtcatc                                                                              >  2:188209/1‑150 (MQ=255)
                                                                          aGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATcac                                                                            <  2:1161481/151‑1 (MQ=255)
                                                                          aGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtcatctc                                                                            >  5:284619/1‑149 (MQ=255)
                                                                           gTAATCCTTGCGTTCTGTCAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATcaca                                                                           <  4:446044/151‑1 (MQ=255)
                                                                             aaTCCTTGCGGTCTGTAAGCCGGCTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAg                                                                          <  8:222233/150‑1 (MQ=255)
                                                                                ccTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGa                                                                      <  6:750955/151‑1 (MQ=255)
                                                                                   tgagTTATGTAATCCGGTTGAGGAGGAAATCCTGATCGTGCTATTCCGTAAGTTATGCGCCCCAGTCCTCCTCCCCCAAATATCGTTGAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAAt                                                                    <  4:509247/147‑1 (MQ=255)
                                                                                   tGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCGCCCCCAAAAATCCTGCAGATGATGAGTGATCCTGCATTAGGCTATGGCTAGGTGAGGAGATTTTGATGAGAGGGGGAtt                                                                   >  4:1003350/1‑151 (MQ=255)
                                                                                          tGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCGGGCGCCCCCCAATATCCTTCAGATGATGAGTGGTCCGGGGTTGGGCTATGGCGAGGTGGTGCGAGTTTCTTCGCAGGGGGATTTTGATGt                                                            >  5:177882/1‑151 (MQ=255)
                                                                                           gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTatgatg                                                             <  6:1178831/149‑1 (MQ=255)
                                                                                           gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTatgat                                                              <  6:851022/148‑1 (MQ=255)
                                                                                           gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTatgat                                                              <  2:5487/148‑1 (MQ=255)
                                                                                           gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGtt                                                           <  1:306609/151‑1 (MQ=255)
                                                                                           gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATCCTTGAGGGGCGGAGTGGGGCTGCAGTTGGCTGTGGGCAGGTGCTCAGGGTTTCTTCACAGGGGGTGTATGATGtt                                                           >  3:306667/1‑151 (MQ=255)
                                                                                            tAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCCAATATCCTTGAGATGATGTGTGTTCCTGGAGTTGGCTAGGGCTGGGTGTTGAGCTTTTGTTCACAGGGGACTTATGATGttt                                                          >  5:674047/1‑150 (MQ=255)
                                                                                             aaGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATACTTCAGATGATGAGTGATCCTGCAGTAGGGTATGGCAAGGTGATCAGATTTTGATAAAAGGGGGATTATGATGTTat                                                         >  6:1299009/1‑150 (MQ=255)
                                                                                             aaGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCGCCCCCAAATATAGTTGTGATGGTGAGTGATCCTGGATTTGGCGATGGCATGGTGATGAGATTTTCTTACGAGGGGAATGATgttgtttt                                                         >  3:597778/1‑149 (MQ=255)
                                                                                              aGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTaaa                                                        <  5:64125/151‑1 (MQ=255)
                                                                                               gCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCGCCCCCCAAAATCCTTCTGATGATGAGTGCTGCGGCGTTTGGGGATGGGAAGGTGTTGCGATTTTTATGACAGGGGAATTATGATGTTAAAt                                                       >  3:409974/1‑151 (MQ=255)
                                                                                                 cGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCTCCCCCAAATAACAATAAAATGATGAGTGATCCTGCATTATGCTATGGCAAGGTGATCAGATTTTGATCACAGGGGAATTATTATGTTAAATg                                                      >  7:209566/1‑150 (MQ=255)
                                                                                                  ggTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCa                                                    <  1:1441300/151‑1 (MQ=255)
                                                                                                  ggTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATAACCATCAGATGATGAGTGATACTGCATTAGGCTATGGGAAGGTGATCAGATTTTAATCACAGGGGAATTATGATGTTAAAtta                                                     >  4:715317/1‑148 (MQ=255)
                                                                                                   gTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCAt                                                   <  8:412705/151‑1 (MQ=255)
                                                                                                     tCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATg                                                  <  4:77488/150‑1 (MQ=255)
                                                                                                      cAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATgg                                                 <  3:251102/150‑1 (MQ=255)
                                                                                                         gaggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCGGCGCCCCCCAATATCCTTCAGATGATGCGTGCTCGGGCGGTAGGGGATGGGAAGGGGAGCAGCGTTTGATCGCAGGGGAGTGATGATGTTTAATGCATGGCAcc                                             >  7:1470603/1‑151 (MQ=255)
                                                                                                             aaaTCCTGATCCGGCTATTCCAAAAGTTATCCGCCCCCGGGCGCCTCCCCCAAATATACTTCAGATGATGAGTGATGCTGCATTAGGGTATGGGAAGGGGATCAGATTTTGATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGcc                                         >  8:1492757/1‑151 (MQ=255)
                                                                                                                 ccTGATCCGGCTATTCCGAAAGTTATCCGCGCGCGGCCGCCGCCCCCCAAGATCCTTCAGATGATGAGTGATGCTGGATTAGGGTATGGCAAGGTGATCAGATTTTGATCACAGGGGAATTATGTTGTTTAATGCATGGCaggtgcagggg                                     >  4:1460985/1‑141 (MQ=255)
                                                                                                                 ccTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCGCCTCCCCCAAATATACGTCAGATGATGAGTGATCCTGCATTAGGGTATGGGAAGGTGATAAGATTTTCATAACAGGGGAATTATGATGTTAAATTAATGGCAGCTGCCGGTg                                     >  4:872422/1‑151 (MQ=255)
                                                                                                                 ccTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCggg                                      <  2:1277702/150‑2 (MQ=255)
                                                                                                                 ccTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGt                                      <  1:108286/150‑1 (MQ=255)
                                                                                                                 ccTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGt                                      <  2:87128/150‑1 (MQ=255)
                                                                                                                         ggCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCCCAAATCCTTCAGATGATGAGTGATCCTGCATTAGGGTATGGCAAGGGGATCAGATTTTCATCACAGGGGAATTATTATGTTAAATGAATGGCACCTGCCGGTGCCCCCCtt                             >  3:1049271/1‑151 (MQ=255)
                                                                                                                          gCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCAttt                            <  5:490044/151‑1 (MQ=255)
                                                                                                                             aTTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGtt                         <  1:1009574/151‑1 (MQ=255)
                                                                                                                              ttCCGAAAGTTATCCGCCCCCGGCCGCCTGCCCCAAAAAAAAATCAGATGATGAATTATACTGCATTAGGCTATTGAAAGGTGATCAGATTTTCATAAAAGGGGAAGTATTATGTAAAATGAAAGGAACCTGACGGTGACCACATTAGta                         >  8:1573294/1‑149 (MQ=255)
                                                                                                                                ccGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTaa                       <  7:1012489/150‑1 (MQ=255)
                                                                                                                                       ttATCCGCCCCCGGCCGCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCAATAGGCTATGGCAAGGTGGTCAGATTTTAATCACAGGGGAAATATGATGTTAAATGAATGGCACCTCCCCGAGCCCCCAAATGTTAAAaaaaga                >  6:1670326/1‑149 (MQ=255)
                                                                                                                                        tATCCGTCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagca               <  4:1631399/150‑1 (MQ=255)
                                                                                                                                         aTCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagcaaa             <  7:453381/151‑1 (MQ=255)
                                                                                                                                         aTCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagcaa              <  2:403360/150‑1 (MQ=255)
                                                                                                                                             gCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATAAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCACCATTTGTTAAACAAACCAAAGAt          >  5:137098/1‑150 (MQ=255)
                                                                                                                                              cccccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATc         >  2:1100659/1‑150 (MQ=255)
                                                                                                                                              cccccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCa        <  5:109766/151‑1 (MQ=255)
                                                                                                                                                  cGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTg    >  8:1175954/1‑151 (MQ=255)
                                                                                                                                                   gTCCTCCTCCCCCAGATATCCTTCAGATGATGAGTGATCCTGCATGAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTG‑CCCCATTTTTTTAAcaaacgaaagcagcactggg  >  5:358515/1‑141 (MQ=255)
                                                                                                                                                   gTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGc   >  2:1418781/1‑151 (MQ=255)
                                                                                                                                                   gTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGc   >  1:1151826/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TTTCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCT  >  NC_000913/422281‑422579

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: ‑