Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913422,4320TG26.4% 745.0 / 74.9 372intergenic (+73/‑86)proY/malZproline‑specific permease/maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base T (77/194);  new base G (97/0);  total (178/194)
Fisher's exact test for biased strand distribution p-value = 6.44e-41
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 389 total aligned reads displayed.

TCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCAT  >  NC_000913/422283‑422582
                                                                                                                                                     |                                                                                                                                                      
tCGCGTGGATTGTTGTGATGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcct                                                                                                                                                        <  5:1804193/150‑1 (MQ=255)
   cgTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctc                                                                                                                                                    <  5:1644486/151‑1 (MQ=255)
    gTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGCGCCGtccgccgcc                                                                                                                                                   >  8:1882025/1‑151 (MQ=255)
      ggATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCGcccc                                                                                                                                                 >  8:457185/1‑151 (MQ=255)
       gATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCCTTCTGTAAGCCGGTTAAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCTccccc                                                                                                                                                >  1:360769/1‑151 (MQ=255)
        aTTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCTcccccc                                                                                                                                               >  7:149964/1‑150 (MQ=255)
           gttgtGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCaaa                                                                                                                                             <  3:1376125/150‑1 (MQ=255)
              gtgCTGTTGATTGGCTGGATCTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATc                                                                                                                                         <  2:647648/151‑1 (MQ=255)
                gCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCt                                                                                                                                       <  7:110338/151‑1 (MQ=255)
                gCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCGGCCCCCCAATATccg                                                                                                                                       >  8:1576042/1‑150 (MQ=255)
                   gtcgATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCa                                                                                                                                    <  7:1729064/148‑1 (MQ=255)
                   gTTGATTTGCTTGATGTTTAAACGCCGCCACGATCTTCAGCTGGCTTAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCtt                                                                                                                                      <  1:152311/149‑1 (MQ=255)
                   gTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCGTGCGCCCCCAAATATCCTTc                                                                                                                                     >  7:1824110/1‑150 (MQ=255)
                     tGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAg                                                                                                                                   <  1:929155/150‑1 (MQ=255)
                      gATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAgat                                                                                                                                 <  3:233681/151‑1 (MQ=255)
                           gCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCTCCTCCCCCAACGATCCGGCAgctgaggg                                                                                                                            >  6:1513854/1‑148 (MQ=255)
                           gCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCCAATATACGTCAGatgatgt                                                                                                                            >  3:1166622/1‑150 (MQ=255)
                                aTGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGTAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCAAATATGCTTCAGATGATGAGTggt                                                                                                                       >  3:585588/1‑149 (MQ=255)
                                aTGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                                                                                                       <  1:1684858/151‑1 (MQ=255)
                                 tGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCCCCGCCGGGCTCGTCGCCCCAATATCAGTCAGAAGCTGAGTGCTc                                                                                                                      >  7:356044/1‑151 (MQ=255)
                                 tGTTTAAACGCCGCCACGATCGTCAGCTGACTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCGCGTCCGGCTCCCCCCAATATGCTTCAGAGGAGGAGTGTTc                                                                                                                      >  1:445919/1‑151 (MQ=255)
                                  gTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGGCTCCCCCAAATATCCTTCAGATGATGAGTGAtac                                                                                                                     >  7:849584/1‑149 (MQ=255)
                                    ttAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCt                                                                                                                    <  7:635449/150‑1 (MQ=255)
                                     tAAAAGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCCACCACCGCCCCCAAATAACCTTCAGATCATGAGGGATCCTgg                                                                                                                  >  8:1401149/1‑150 (MQ=255)
                                      aaaCGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCa                                                                                                                 <  4:152303/151‑1 (MQ=255)
                                      aaaCGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCAAATATCCGTCAGCTGATGAGTGATGCTgg                                                                                                                  >  5:1338255/1‑149 (MQ=255)
                                          gccgccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCGGCCTCCCCCAAATAACCTTCAGATGATGAGTGATCCTggagtat                                                                                                             >  6:428168/1‑145 (MQ=255)
                                          gccgccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCCGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGGCGCCCCCAAATATCCTTCAGATGATGAGTTATCCTGCAtt                                                                                                               >  2:1457423/1‑149 (MQ=255)
                                           ccgccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                                                                                                            <  6:1829572/151‑1 (MQ=255)
                                            cgccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGc                                                                                                           <  8:674073/151‑1 (MQ=255)
                                              ccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTa                                                                                                         <  5:1730216/151‑1 (MQ=255)
                                                aCGATCGTCAGCTGGCTGAAAACCATTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCGCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATg                                                                                                       <  7:557090/151‑1 (MQ=255)
                                                aCGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATg                                                                                                       >  7:1000030/1‑151 (MQ=255)
                                                aCGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCAGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATg                                                                                                       <  7:1559229/151‑1 (MQ=255)
                                                  gATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGc                                                                                                     <  4:1711539/151‑1 (MQ=255)
                                                  gATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAAAATCCTTCAAATGAAGAGTGATACAGCATTAGGCTATGGc                                                                                                     >  4:142277/1‑151 (MQ=255)
                                                       tCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCGGCCGCCACCAAATATCCGGCAGATGCTGAGTGATCCTGCAGTAGGGTATGGGAAggg                                                                                                >  6:1402821/1‑150 (MQ=255)
                                                          gctggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGCCCCGGCCTCCGCCCCCACAGATCCTGCAGAGGATGAGTGGTGCTGCGTTAGTCGAGGGCAAgggggt                                                                                             >  1:534558/1‑147 (MQ=255)
                                                             ggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATc                                                                                            <  3:384944/149‑1 (MQ=255)
                                                                tgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCGCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAtt                                                                                       >  1:1341918/1‑151 (MQ=255)
                                                                 gAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAt                                                                                        <  3:337724/149‑1 (MQ=255)
                                                                 gAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAt                                                                                        <  6:817429/149‑1 (MQ=255)
                                                                  aaaaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAtttt                                                                                     <  8:1226010/151‑1 (MQ=255)
                                                                   aaaCCAGTAATCCTTGCGTGCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtc                                                                                    <  1:659829/151‑1 (MQ=255)
                                                                     aCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtcat                                                                                  <  5:1368156/151‑1 (MQ=255)
                                                                     aCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtcat                                                                                  <  2:1721640/151‑1 (MQ=255)
                                                                         gTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATcaca                                                                              <  1:1746406/151‑1 (MQ=255)
                                                                           aaTCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCTCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAg                                                                             <  4:965138/150‑1 (MQ=255)
                                                                           aaTCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAgg                                                                            <  1:1839831/151‑1 (MQ=255)
                                                                                  gCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGCCCCGTCCGCCTCCCCCCAATATCCGTCAGATGATGAGTGCTGCTGCATTAGGCTATGGCAAGGGGATCAGATTTTCATAACAGAGGAAtt                                                                      >  8:739895/1‑150 (MQ=255)
                                                                                  gCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAAtt                                                                      <  5:1116260/150‑1 (MQ=255)
                                                                                  gCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTa                                                                     <  1:1014971/151‑1 (MQ=255)
                                                                                   cGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGCCCCGGCCTCCTCCCCCAAATATACTTCAGATGATGAGTGATGGTGGATTAGGCTATGGCTAGGTGGTCAGATTTTAATGAGAGGGGAATTat                                                                    >  2:1155690/1‑151 (MQ=255)
                                                                                   cGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTat                                                                    <  3:704735/151‑1 (MQ=255)
                                                                                    gTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCGGCGGCCCCCAATATCCTTCAGATGGTGAGTGGTCCTGGAGTAGGCTATGGGCAGGTGCTCAGATTTTCTTGACAGGGGAGttttt                                                                   >  3:1081717/1‑148 (MQ=255)
                                                                                        ngTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGGCCCACCGCCTCCGCCAAATATCAGTGAGAAGAAGAGTGTTAGTGGAATTGTCGATGGCAAGGGGATGAGATTTTCAAAACAGGGGTATTttgattt                                                               >  2:418758/2‑149 (MQ=255)
                                                                                         gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGtt                                                              <  8:923823/151‑1 (MQ=255)
                                                                                         gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGt                                                               <  2:795886/150‑1 (MQ=255)
                                                                                         gTAAGCCGGATCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGtt                                                              <  2:236623/151‑1 (MQ=255)
                                                                                           aaGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTaa                                                            <  5:226593/151‑1 (MQ=255)
                                                                                           aaGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTaa                                                            <  4:122513/151‑1 (MQ=255)
                                                                                             gCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAAt                                                          <  2:593258/151‑1 (MQ=255)
                                                                                                ggTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAAGATCCTACAGATGATGAGTGCTCGGGCATTAGGCGATGGCAAGGGGATCAGAGTGTAGTAACAGGTGAATTATGTTGTTtaaatc                                                        >  3:1809170/1‑147 (MQ=255)
                                                                                                  ttCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCTCCCCCCAATAAACTTCAGATGATGAGTGATCCTGCATTATTCTATGGCTAGGTGAGGAGATTTTAATCACAGGGGATATATTATGTTAAATGCATg                                                     >  6:1215336/1‑151 (MQ=255)
                                                                                                    cAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCGCCCCCACATATCCTGAAGATGATGAGGGATCGTGGAGTAGGCTATGGCAAGGTGATGAGCTTTTCATGAGAGGGGAATTATGATGTTTAATGCATGGc                                                   >  6:1613179/1‑151 (MQ=255)
                                                                                                     aggaggAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTCAAATGCATGGCa                                                  <  1:844369/151‑1 (MQ=255)
                                                                                                     aggaggAAATCCTGATCCGGCTATTCCGAAAGATATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCa                                                  <  2:999770/151‑1 (MQ=255)
                                                                                                        aggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCGCCCCCAAATATCCTTCAGAAGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCACATTTTAATCACAGGGGAATTATGAGTTTTAATTGATGGCAc                                                 >  4:247880/1‑149 (MQ=255)
                                                                                                           aaaTCCTGATCCGGCTATTCCGAAAGTTATCCGCGCGCGGCGGGCGGGCGCACAGATCGTTCAGCAGGAGAGTCGTCCTGCATTAGGCGAGGGCAAGGGGGGCGGATTTTGCTCACGGGGGGATTTTGTTGTTAGTTGCATGGcgggtggg                                            >  2:1806320/1‑144 (MQ=255)
                                                                                                               ccTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTAGGGCAAGGGGATCAGATTTTCATCACACGGGAATTATGATGTTAAATTCATGGCCCCTGCCgggg                                        >  7:1365400/1‑149 (MQ=255)
                                                                                                                  gATCCGTCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGc                                       <  4:1845278/149‑1 (MQ=255)
                                                                                                                    tCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCCAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCGATGGGAAGGGGATCAGCTTTTGATAACAGGGGAATTATGATGTTAAATTCAATGCCCCGGCCGGTGccccc                                   >  1:1391215/1‑151 (MQ=255)
                                                                                                                     ccGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCTCCCCCAAATATACTTCAGATGATGAGTGGTCCTGCATTAGGCTATGGGAAGGTGCTCAGATTTTCATCACAGGGGAATTATGATGTTAAATGAATGGAACCTGCCGGTGccccc                                   >  3:1557193/1‑150 (MQ=255)
                                                                                                                     ccGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCAAATATACTTCAGATGATGAGTGATCCTGCATTAGGCTATGGAAAGGTGATCAGATTTTCATAACAGGGGAATTATGATGTTAAATTCATGGCACCTGCCGGTGccccc                                   >  5:151817/1‑150 (MQ=255)
                                                                                                                     ccGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAAAATCCTTAAGATGATGAGTGATCCTGCAGTAGGGTATGGCAAGGTGATCAGATTTTAAACACAGGGGAATTATGATGTTAAATGCAAGGCACCTGCCGGGGccccc                                   >  4:107524/1‑150 (MQ=255)
                                                                                                                     ccGGCTATTCCGAAAGTTATCCGCCCCCCGCCTCCTCCCCCAAATATCCTTAAGATGATGAGTGATCCTGATTTAGGCGATGGCAAGGTGATCAGATTTTAATCAAAGGGGAATTTTGATTTTATATTCATGGCACCTCCCGGGGcgccc                                   >  7:458455/1‑150 (MQ=255)
                                                                                                                      cGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCAAAAATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGGAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCCCCTGCCGGGGgccccct                                 >  8:39558/1‑149 (MQ=255)
                                                                                                                        gCTATTCCGAAAGGTATCCGCCCCCGGCCGCGGCCCCCCCCTATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCAAAGGGGAATTATGATGTTTAATGAATGGCACCTGCCGCTGCCCCCAttt                               >  3:1334300/1‑151 (MQ=255)
                                                                                                                          tATTCCGAAAGTTATCCGGCGCCGGCCGCCGCCCCCCAAGATCCTTCAGATGATGAGTGATGGTGCATTAGGGGATGGGAAGGTGAGGAGATTTTGATGACAGGGGAATTTTGATGTTAAATTGATGGCAGCTGCCGGTGGCGCCAGTTGt                             >  6:841849/1‑151 (MQ=255)
                                                                                                                              ccGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATACTTCAGATGATGAGTGATCATGCATTATGCTATTGCAAGGTGATCAGATTTTATTCAAAGCGGAATTATTATGTTTAATTAATGGCACCTTCCGCGGCCCCCAATAGTTaa                          >  5:439576/1‑150 (MQ=255)
                                                                                                                              ccGAAAGTTATCCGCCCCCGGCCTCCGCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGTCTATGGCAAGGTGATCAGATTTTAATCAAAGGGGAATTATGATGTTAAATTAATGTCAACTGACGGTGCCCCCATTTGTTaa                          >  1:269805/1‑150 (MQ=255)
                                                                                                                                 aaaGTTATCCGCCCCCGGGCGCCTCCCCCAAAAATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGGAAGGTGATCAGATTTTCATGACACGGGCATTTTTATGTTTTATGAATGGCGCCTGCCGGGGCCGCCCTTTGTTAAAca                       >  8:112037/1‑150 (MQ=255)
                                                                                                                                  aaGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaa                     <  5:879482/151‑1 (MQ=255)
                                                                                                                                     ttATCCGCCCCCGGCCGCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTAATCACAGGGGAATTATGATGTTAAATGAATGGCAACTGCCGGTGCCCCCATGTGTTAAAcaaaga                   >  6:1025835/1‑149 (MQ=255)
                                                                                                                                      tATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCTGGTGCCCCCATTTGTTAAAcaaagca                  <  7:901874/150‑1 (MQ=255)
                                                                                                                                      tATCCGCCCCCGGCCTCCTCCCCCAAATATCCGTCAGATGATGAGTGATCCGGCATTAGGCTATGGCCAGGTGATGAGATTTTCATCACAGGGGAATTATGATGTTAAATGCAAGGCACCTGCCGGGGCCCCCAGTTGTTAAAcaaagaa                  >  8:27052/1‑148 (MQ=255)
                                                                                                                                       aTCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagcaa                 <  3:820566/150‑1 (MQ=255)
                                                                                                                                         ccGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGa              <  5:1445207/151‑1 (MQ=255)
                                                                                                                                         ccGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGa              <  6:291031/151‑1 (MQ=255)
                                                                                                                                         ccGCCCCCGGCCTCCGCCCCCAAATATCCTTCAGATGATGAGTGATCCTGGATTAGGCTATGGCAAGGTGATCAGATTTTCATAACAGGGGAATTATTATGTTAAATGCATGGCACCTGCCGGTGCGCCCATTTGTTAAAcaaagcaa                 >  6:550383/1‑148 (MQ=255)
                                                                                                                                          cGCCCCCGTCCTCCTCCCCCAAATATACTTCAGATGATGAGTGATCCTGCATTAGGCTATGTCAAGGTGATCAGATTTTAATAACAGGGGAATTATGATGTTAAATGCATGGAACCTGCCGGTGCCCCCAATTGTTAAAaaaagaaaaca              >  7:1192871/1‑148 (MQ=255)
                                                                                                                                           gCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCAATTGTTAAAcaaagcaa                 >  8:1683057/1‑146 (MQ=255)
                                                                                                                                             ccccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCa           >  7:998173/1‑150 (MQ=255)
                                                                                                                                               ccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACt        <  5:608492/151‑1 (MQ=255)
                                                                                                                                               ccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACt        <  7:510781/151‑1 (MQ=255)
                                                                                                                                                ggTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTATGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATTTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTg       <  1:1193208/150‑1 (MQ=255)
                                                                                                                                                 gTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGc      >  4:365616/1‑151 (MQ=255)
                                                                                                                                                  tcctcctccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCt     >  6:74839/1‑151 (MQ=255)
                                                                                                                                                   cctcctccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTc    <  4:705328/151‑1 (MQ=255)
                                                                                                                                                   cctcctccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTc    <  3:946722/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCAT  >  NC_000913/422283‑422582

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