Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913422,4320TG25.0% 767.2 / 61.3 354intergenic (+73/‑86)proY/malZproline‑specific permease/maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base T (74/190);  new base G (88/0);  total (164/190)
Fisher's exact test for biased strand distribution p-value = 4.43e-38
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 381 total aligned reads displayed.

TCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCAT  >  NC_000913/422283‑422582
                                                                                                                                                     |                                                                                                                                                       
tCGCTTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcct                                                                                                                                                         <  7:452059/150‑1 (MQ=255)
tCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCGCCGgcctc                                                                                                                                                        >  8:529231/1‑151 (MQ=255)
tCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcct                                                                                                                                                         <  3:1513832/150‑1 (MQ=255)
tCGCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcct                                                                                                                                                         <  4:558436/150‑1 (MQ=255)
  gcgTGGATTGTTGAGCGGTTGATTGGCTGGATGTTTAAACGACGCCAAGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGTTTCAGGAGGAATTCCTGATCCGGATATTCCGAAAGTTATCCGCCCCCGtcctcct                                                                                                                                                      <  2:1034764/151‑1 (MQ=255)
         ttgttgTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCGCGGtcctgggcccccc                                                                                                                                                >  6:1575806/1‑141 (MQ=255)
         ttgttgTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCGCCCCCac                                                                                                                                               >  5:886445/1‑150 (MQ=255)
            ttgtGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCGCCTCGCCCCAat                                                                                                                                             >  4:709834/1‑150 (MQ=255)
                   gTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTc                                                                                                                                      <  5:1324704/150‑1 (MQ=255)
                   gTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCa                                                                                                                                     <  1:1402658/151‑1 (MQ=255)
                     tGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAg                                                                                                                                    <  1:1672228/150‑1 (MQ=255)
                     tGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAg                                                                                                                                    <  3:1047790/150‑1 (MQ=255)
                      gATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGGTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCGTCCTCCCCCAAATATCCTTCAgtt                                                                                                                                  >  3:77664/1‑149 (MQ=255)
                      gATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCCCATATCCGGCAga                                                                                                                                   >  7:1281494/1‑150 (MQ=255)
                       aTTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCGGTCGGGCTCCCCCAAatatagataatatg                                                                                                                                 >  8:345733/1‑141 (MQ=255)
                         tGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCAAATATCCTTCAgatgat                                                                                                                               >  2:710286/1‑151 (MQ=255)
                          ggCGGGATGTTTAAACGCCGCCACGATCCTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGatgatg                                                                                                                              <  7:671640/151‑1 (MQ=255)
                           gCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGAtgatga                                                                                                                             <  5:669693/151‑1 (MQ=255)
                             tGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTCTCCGCCCCCGTCCGCCGCCCGCAACTATCCTTCAGATGATGAgg                                                                                                                           >  4:1056956/1‑150 (MQ=255)
                               gATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGa                                                                                                                         <  1:783891/151‑1 (MQ=255)
                                aTGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGAt                                                                                                                        <  6:1410076/151‑1 (MQ=255)
                                aTGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGGCGCCGCCCCCAAATATCCTTCAGATGGTGAGTGAt                                                                                                                        >  7:1215150/1‑151 (MQ=255)
                                  gtgtAAACGCCGCGACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATc                                                                                                                       <  7:307569/147‑1 (MQ=255)
                                     tAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTg                                                                                                                    <  1:156185/150‑1 (MQ=255)
                                         cgccgcCTCGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCACGAGGAAATCCTGATCCGGCTATTCCGAAAGATATCCGCCCCCGTCCTCCTCCCCCAAATCTCCTTCAGATGATGAGTGATCCTGCATTa                                                                                                               <  4:971692/151‑1 (MQ=255)
                                         cgccgcCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTa                                                                                                               <  6:774610/151‑1 (MQ=255)
                                           ccgccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                                                                                                             <  5:764321/151‑1 (MQ=255)
                                            cgccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCGCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGc                                                                                                            <  5:822950/151‑1 (MQ=255)
                                              ccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTa                                                                                                          <  3:292696/151‑1 (MQ=255)
                                               cACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTAt                                                                                                         <  6:1067562/151‑1 (MQ=255)
                                               cACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATACTTCAGATGATGAGTGTTCCTGCGTTAgggggt                                                                                                         >  7:1072020/1‑147 (MQ=255)
                                               cACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTGCTGTAAGCCTGTTCAGGAGGAAATCCTCATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTAt                                                                                                         <  6:1326578/151‑1 (MQ=255)
                                                     cGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAg                                                                                                   <  5:1087909/151‑1 (MQ=255)
                                                      gTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCTCCTCGCCCAAAAATCCTTCAGATGATGAGTGATGCTGCATTAGGCTATGGCAAgg                                                                                                  >  1:609809/1‑151 (MQ=255)
                                                      gTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAgg                                                                                                  <  8:1479643/151‑1 (MQ=255)
                                                         aGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCCACTATCCTTCAGATGATGAGTGATCCTGCATTAGGGTATGGCATGGTgt                                                                                               >  2:649164/1‑150 (MQ=255)
                                                          gctggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGAt                                                                                              <  2:427901/151‑1 (MQ=255)
                                                           ctggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAAAATACTTCAGATGATGAGTGATCCTGCATTAGGGGATGGCAAGGTGATc                                                                                             >  8:105017/1‑151 (MQ=255)
                                                            tggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCCAATATCCTTCAGATGATGAGGGCTCCTGCATTAGGCTATGGCAAGGTGATc                                                                                             >  1:98609/1‑150 (MQ=255)
                                                              gctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGa                                                                                          <  7:192022/151‑1 (MQ=255)
                                                                tgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCGCCGGCCGCCTCCCCCCAATATGCTTCAGAAGAGTAGTGATGCTTCATTAGGGTATGGGAAGCTgtagagatt                                                                                        >  3:201533/1‑143 (MQ=255)
                                                                tgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCGTCCGCCCCCAAATATCCTTAAGATGATGAGGGATCCTGCATTAGGCTATGGCAAGGTTAGCAGAtt                                                                                        >  8:642533/1‑151 (MQ=255)
                                                                tgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGGAAGGTGATCAGAtt                                                                                        >  7:987987/1‑151 (MQ=255)
                                                                 gAAAACCAGTAATCCTTGCGTTCTGTAAGGCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAt                                                                                         <  2:47347/149‑1 (MQ=255)
                                                                 gAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAtt                                                                                        <  5:485512/150‑1 (MQ=255)
                                                                 gAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAt                                                                                         <  4:1354556/149‑1 (MQ=255)
                                                                   aaaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCCAATATCCTTCACATGATGCGGGATCCTGCCTTAGGCTGTGGCATGGGGATCAGAtttt                                                                                      >  5:1382713/1‑150 (MQ=255)
                                                                    aaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtca                                                                                    <  5:334953/151‑1 (MQ=255)
                                                                      ccAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAGATCCTGATCCGTCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtcatc                                                                                  <  5:536689/151‑1 (MQ=255)
                                                                      ccAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtcatc                                                                                  <  5:977149/151‑1 (MQ=255)
                                                                        aGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATcac                                                                                <  2:1161481/151‑1 (MQ=255)
                                                                           aaTCCTTGCGGTCTGTAAGCCGGCTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAg                                                                              <  8:222233/150‑1 (MQ=255)
                                                                            aTCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAggg                                                                            <  1:1289625/151‑1 (MQ=255)
                                                                              ccTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGa                                                                          <  5:1629535/151‑1 (MQ=255)
                                                                              ccTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAgggg                                                                           <  3:1329063/150‑1 (MQ=255)
                                                                               cTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAAAATCATACAGATGAGGAGTGATCCTGCATTATGGTATGGCATGGTGATCAGATTTTGATACAAGGGGa                                                                          >  4:442166/1‑150 (MQ=255)
                                                                                 tgagTTATGTAATCCGGTTGAGGAGGAAATCCTGATCGTGCTATTCCGTAAGTTATGCGCCCCAGTCCTCCTCCCCCAAATATCGTTGAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAAt                                                                        <  4:509247/147‑1 (MQ=255)
                                                                                 tGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAAt                                                                        <  3:368201/150‑1 (MQ=255)
                                                                                  gCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCCAATATCCTGCAGATGCTGCGTGATGGGGGATTAGGCGATGGCGAGGGGGTCAGAGTGTGATGACAGGGGaagtt                                                                      >  5:1144767/1‑148 (MQ=255)
                                                                                   cGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATCCTTCAGATGATGAGTGAGCCTGCATTTGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTat                                                                     >  5:273279/1‑151 (MQ=255)
                                                                                     ttCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCGTCAACCAAAAATCGTTCAGATTATGAGTGATGCTGAATTATGCAATGGAAATGTGATCAGATTTTAATAACAGGGGAATAatg                                                                    >  6:95880/1‑150 (MQ=255)
                                                                                      gcTGTAAGCCGGCTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTatgat                                                                  <  6:529784/150‑1 (MQ=255)
                                                                                         gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTatgat                                                                  <  2:5487/148‑1 (MQ=255)
                                                                                         gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGtt                                                               <  8:1102160/151‑1 (MQ=255)
                                                                                         gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGt                                                                <  6:1196589/150‑1 (MQ=255)
                                                                                            aGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGCCCCGGCCGCCTCCCCCAAATATCCTTCAGATGGTGAGTGCTCGTGCGTTAGGGGATGGCGAGGTGGGCAGATTTTGGTCACAGGGGAGGTTTGATGTTaa                                                             >  8:1316087/1‑150 (MQ=255)
                                                                                            aGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTaaa                                                            <  3:895718/151‑1 (MQ=255)
                                                                                            aGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCAAATATCCTTGAGATGATGAGTGATCCTGCATTAGGGTATGGGAAGGGGATCAGATTTTCATCAGAGGGGATTTATTATGttta                                                             >  6:148940/1‑148 (MQ=255)
                                                                                               cgggtgaggaggAAATCCTGATGCGGCTATTCCGAAAGTTATCCGCCACCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGATTTATGATGTTAAATGc                                                         <  3:954441/145‑1 (MQ=255)
                                                                                                ggTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCGCCGGCCTCCTCCCCCAAATATACTTCAGATGCTGAGTGGTGGGGCAGGAGGCGATGGGAAGGGGATCCGATTTTGATCACAGGGGAATTATGATGTTAATTga                                                         >  6:403794/1‑149 (MQ=255)
                                                                                                ggTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCGCCCCCAAAAATACTACAGATGATGAGTGGTGCTGCATTAGGCTTTGGCAAGGTGATCAGATATTCTTAAGAGGGGAAGTATTATGTTAAAtta                                                         >  1:401086/1‑148 (MQ=255)
                                                                                                 gTTCAGGATGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATGCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCAt                                                       <  4:797681/151‑1 (MQ=255)
                                                                                                   tCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCCCCCCCGGCCGCCTCCCCCACAAATCCTAAATATGATGAGTGAGAGTGGATTAGGCTATGGCGAGGTGATCAGAGTTTGTTGACAGGGGAATTATGATTTTTTATGCATgg                                                    >  2:353478/1‑151 (MQ=255)
                                                                                                    cAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGc                                                   <  4:381179/151‑1 (MQ=255)
                                                                                                     aggaggAAATCATGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCa                                                  <  6:333053/151‑1 (MQ=255)
                                                                                                         ggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCGCCCCCAAAAAACCGTCCCAAGAAGAGTGATCCTCCATTAGGCTATGTTAATGTGGTCAGATTATCATAACACGGTAATTATTATGTTAAATTAATGGAAcatc                                              >  5:1477524/1‑148 (MQ=255)
                                                                                                             aTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCCGCCGCCGCCCCCAAATATCCTTCAGATGTTGAGTGCTCCTGGATTAGGCTATGGCAAGGTGAGCAGATTTTCATCACAGGGGTATTATGTTGTTTAATTCATGGCACCTGCCgg                                          >  3:1176844/1‑151 (MQ=255)
                                                                                                               ccTGATCCGGCTATTCCGAAAGTTATCCGCGCGCGGCCGCCGCCCCCCAAGATCCTTCAGATGATGAGTGATGCTGGATTAGGGTATGGCAAGGTGATCAGATTTTGATCACAGGGGAATTATGTTGTTTAATGCATGGCaggtgcagggg                                        >  4:1460985/1‑141 (MQ=255)
                                                                                                               ccTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCggg                                         <  2:1277702/150‑2 (MQ=255)
                                                                                                               ccTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGt                                         <  1:108286/150‑1 (MQ=255)
                                                                                                               ccTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCGCCCCCAAATATCATTGAGATGATGAGTGGTCCTGCGGTTGGGTATGGCAAGGTGGTAAGATTTTCTTCACAGGGGAATTTTGATGTTAAATTCATGGCCCCTGCCGGTg                                        >  8:443336/1‑151 (MQ=255)
                                                                                                                 tGATCCGGCTATTCCGAAAGTTATCCGCGGCCGGGCGCCGCCCCCCACCATCCTGCAGATGATGAGTGATCCGGCATTATGATATTGATAGGTGAGCAGATTTTCTTCCGAGGGGAAGTATGTTTTTAAATTAATGGCACATTCCCGTGcc                                      >  3:1076086/1‑151 (MQ=255)
                                                                                                                        gCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCAttt                               <  7:1405172/151‑1 (MQ=255)
                                                                                                                           aTTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTg                              <  6:270942/149‑1 (MQ=255)
                                                                                                                           aTTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGtt                            <  2:811726/151‑1 (MQ=255)
                                                                                                                           aTTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGtt                            <  1:1009574/151‑1 (MQ=255)
                                                                                                                            ttCCGAAAGTTATCCGCCCCCGGCCGCCTGCCCCAAAAAAAAATCAGATGATGAATTATACTGCATTAGGCTATTGAAAGGTGATCAGATTTTCATAAAAGGGGAAGTATTATGTAAAATGAAAGGAACCTGACGGTGACCACATTAGta                            >  8:1573294/1‑149 (MQ=255)
                                                                                                                               cGAAAGTTATCCGCCCCCGTGCTCCTCCCCCAAATTACATAATTATAATAAGGGATCGTAAATTAGGATATGGGAAGGTGTTCAGATTTTATTCAAAGGCTAATTTTCTTTTTAAATCATAGGAACCGGCCGGGGCCGCCTATAGTTaaa                         >  7:900674/1‑150 (MQ=255)
                                                                                                                               cGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAAAATACTTCAGATGATGAGTGATCCTGCATTAGGGGATGGGAAGGTGATGAGATTTTAATGACAGGGGAATTATGATGTTAAATTAATGGCACCTGCCGGTGCGCCCAATTGTTaaa                         >  6:683740/1‑150 (MQ=255)
                                                                                                                                 aaaGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaa                      <  3:506893/151‑1 (MQ=255)
                                                                                                                                  aaGTTATCCGCCCCCGGCCTCCGCCCCCAAATATACTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTAATCACAGGGGAATTATGATGTTAAATGAATGGGACCTGCCGGGGCCCCCCGTTGTTAAAcaa                      >  8:910276/1‑150 (MQ=255)
                                                                                                                                      tATCCGTCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagca                  <  4:1631399/150‑1 (MQ=255)
                                                                                                                                       aTCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagcaaa                <  3:642279/151‑1 (MQ=255)
                                                                                                                                           gCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATc            <  4:742665/151‑1 (MQ=255)
                                                                                                                                             ccccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCa           >  2:1517694/1‑150 (MQ=255)
                                                                                                                                              cccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGCTCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGGGCCCCCATTTGTTAAACAAAGCAAAGATCaa          >  1:290980/1‑150 (MQ=255)
                                                                                                                                                 gTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGc      >  2:1418781/1‑151 (MQ=255)
                                                                                                                                                  gcctcctccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCt     <  1:114531/150‑1 (MQ=255)
                                                                                                                                                   cctcctccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTc    <  2:949289/151‑1 (MQ=255)
                                                                                                                                                   cctcctccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTc    >  2:1653221/1‑151 (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: