Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913422,4350TG25.3% 814.3 / 67.9 385intergenic (+76/‑83)proY/malZproline‑specific permease/maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base T (75/209);  new base G (96/0);  total (176/209)
Fisher's exact test for biased strand distribution p-value = 5.66e-43
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 407 total aligned reads displayed.

GCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTAC  >  NC_000913/422285‑422585
                                                                                                                                                      |                                                                                                                                                      
 cgTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGtcctcctc                                                                                                                                                       <  4:929604/151‑1 (MQ=255)
    ggATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCGcccc                                                                                                                                                    >  8:448339/1‑151 (MQ=255)
          gtgtgCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCCGTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAata                                                                                                                                              <  3:1353535/150‑1 (MQ=255)
           tgtgCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCCCCCCCGTCCTCCTCCCCCCAatat                                                                                                                                             >  4:454595/1‑151 (MQ=255)
               ntGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTACTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCAGCCCCCAAATATCCtt                                                                                                                                         >  8:1359051/2‑151 (MQ=255)
                tGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCGGCCTCCCCCAAATATCCTTc                                                                                                                                        >  4:1221485/1‑151 (MQ=255)
                    gATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCTCCCCCAAATATGCTTc                                                                                                                                        >  6:918953/1‑147 (MQ=255)
                     aTTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAgatg                                                                                                                                   <  7:1065647/151‑1 (MQ=255)
                       tGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAgatga                                                                                                                                  <  1:1463899/150‑1 (MQ=255)
                        ggCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTCTCCGCCCCCGGCCGCCGCCCCCCCATATccggccgatgcgg                                                                                                                                >  4:582334/1‑140 (MQ=255)
                             gATGTTTAAACGCCGCCACTTTCGTCAGCTGGCTGAAAAGCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCTAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGa                                                                                                                           <  1:579346/151‑1 (MQ=255)
                             gATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGa                                                                                                                           <  7:971870/151‑1 (MQ=255)
                                gTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCATTTATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATc                                                                                                                         <  6:1460953/150‑1 (MQ=255)
                                      aCGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCGCCCCCAAATATCCTTCATATGATGAGTGATCCTGCAtt                                                                                                                  >  1:346228/1‑151 (MQ=255)
                                       cgccggCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTATGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTa                                                                                                                 <  2:1161025/151‑1 (MQ=255)
                                         ccgccACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAgg                                                                                                               <  1:1502208/151‑1 (MQ=255)
                                             cACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTAt                                                                                                           <  1:1586647/151‑1 (MQ=255)
                                              aCGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATg                                                                                                          <  5:1100374/151‑1 (MQ=255)
                                              aCGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATg                                                                                                          <  7:34405/151‑1 (MQ=255)
                                              aCGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATg                                                                                                          <  2:643856/151‑1 (MQ=255)
                                               cGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATgg                                                                                                         <  8:1585534/151‑1 (MQ=255)
                                                gATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGc                                                                                                        <  3:634689/151‑1 (MQ=255)
                                                gATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCGCCCCCAACTATCCTTCAGATGGAGCGGGATCCTGGAGTAGGGTATGGc                                                                                                        >  3:1019231/1‑151 (MQ=255)
                                                gATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCGCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTTGGGTATggg                                                                                                        >  3:1652844/1‑150 (MQ=255)
                                                gATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGGCGCCCCCCAATATCCTTCAGATGATGAGGGGTCCGGGATTAGGGTATggg                                                                                                        >  8:522434/1‑150 (MQ=255)
                                                        gctggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGAt                                                                                                <  5:337571/151‑1 (MQ=255)
                                                           ggctgAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAg                                                                                             <  5:591475/151‑1 (MQ=255)
                                                               gAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGGGCCGGCGTCCGCGGCCCAAGATCCGGCAGAGGATGAGTGATCCTGCATTTGGGTATGGGAAGTTGTTCGGCttt                                                                                         >  1:1354568/1‑151 (MQ=255)
                                                               gAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAt                                                                                           <  4:47029/149‑1 (MQ=255)
                                                                aaaaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAtttt                                                                                        <  5:289121/151‑1 (MQ=255)
                                                                aaaaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCGGTGGGAAGGTGAGGAGCtttt                                                                                        >  6:1563435/1‑151 (MQ=255)
                                                                  aaCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtca                                                                                      <  6:785717/151‑1 (MQ=255)
                                                                   aCCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtcat                                                                                     <  1:983131/151‑1 (MQ=255)
                                                                    ccAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCAAATATCCTTCAGATGATGAGGGCTCCTGCATTAGGCTATGGCAAGGTGATGAGATTTtcatc                                                                                    >  3:575375/1‑151 (MQ=255)
                                                                     cAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTtcatca                                                                                   <  8:698700/151‑1 (MQ=255)
                                                                      aGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATcac                                                                                  <  6:890679/151‑1 (MQ=255)
                                                                      aGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCGGCCGCCCCCAAATATCGTTCAGATGTGGAGTGATCCGGCATTAGGCTTGGGCAAGGTGAGCAGATTTTGATaac                                                                                  >  5:185306/1‑151 (MQ=255)
                                                                      aGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCGCCCCCCAATATCCTTCAGATGATGAGTGATCCTGCATTAGGGGATGGCAAGGGGGGCAGAGTTTCATcac                                                                                  >  2:1669802/1‑151 (MQ=255)
                                                                      aGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCAAAAAAGCTTAAGATGATGAGTGATCCTGCATTAGGCGATGGCACGGGGCTCAGATTGtgatccc                                                                                  >  3:235294/1‑149 (MQ=255)
                                                                      aGTAATCCTTGCGTTCTGTAAGCCGCTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATcac                                                                                  <  4:994404/151‑1 (MQ=255)
                                                                        tAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGCCTATTCCGAAAGTTATCCGCCCCCGGGCGCCGCCCCCAAATATCAGTCAGATGATGAGTGCTGGTGGATTAGGCTTTGGGAAGCTGTTCCGATTTtcatcagac                                                                                >  7:1548551/1‑148 (MQ=255)
                                                                         aaTCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAAAATACTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAACGTGATCAGATTTTCATAACAgg                                                                               >  8:1252093/1‑151 (MQ=255)
                                                                          aTCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGGCCCCGGCCGCCGCCCCCCAATCTCCGTCAGATGATGAGTGCTCCTGCAGTAGGGTATGGCAAGGTGGGGCGAGTTtcatcagcggg                                                                              >  3:1468573/1‑146 (MQ=255)
                                                                          aTCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGGCCGGGGGCCGGCGCCCACTATCCTGCCGATGTTGAGGGCTGCTGGGTTAGGCTAGGGCATGGTGATCTGATTTTCATGACAggg                                                                              >  3:668532/1‑151 (MQ=255)
                                                                           tCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCGCCGCCCCCACATAACAGTCACAAGATGCGTGATCATGCAGTAGGGGATGGCAAGGTGAACAGAgttagatgaccgggt                                                                             >  2:366333/1‑139 (MQ=255)
                                                                            ccTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCGTCCCCCCAATATCCGTCAGATGAAGATTGATGCTGCATTAAGCGATGTCGAGGTGAGCAGATTTTGATCACAgggg                                                                             >  2:770409/1‑150 (MQ=255)
                                                                                gCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGATAGTTATCCGCCCCCGTCCTGCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTTGGCTATGGCAAGGTGATCAGATTTTCATCACAGTGGAATTa                                                                        <  6:595353/151‑1 (MQ=255)
                                                                                gCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGGTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATTGCAAGGTGATCAGATTTTCATCACAGGGGAATTa                                                                        <  8:747150/151‑1 (MQ=255)
                                                                                gCGTGCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAAtt                                                                         <  1:1138722/150‑1 (MQ=255)
                                                                                 cGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTat                                                                       <  3:389026/151‑1 (MQ=255)
                                                                                    tCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTatga                                                                     <  3:1582828/150‑1 (MQ=255)
                                                                                     cTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCTCCCCCAAAAATCCTTCAGATGATGAGTGCTCCTGCATTAGGGGATGGGAAGGTGATCCGAGTTTCATCACAGGGGAGTtattatt                                                                   >  5:1359605/1‑150 (MQ=255)
                                                                                      tGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGGGCGGCGGGCTGCCCCCCAGATCCTTCAGATGAGGCGTGCTGGTGGAGGAGGGTATGGCTGGGGGATGGGAGTTTGATGAGAGGGGTTTTTTGATGt                                                                  >  3:919840/1‑151 (MQ=255)
                                                                                       gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTatgat                                                                    <  7:144362/148‑1 (MQ=255)
                                                                                       gTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGt                                                                  <  4:716058/150‑1 (MQ=255)
                                                                                         aaGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTGCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTa                                                                <  7:825386/150‑1 (MQ=255)
                                                                                               gTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCGCCGCCCCCACATATCCTACGGGGCATGAGTGCGCGGGCAGTAGGGTATGGCCACGGGCTCGGATTTTCTTCACCGGGGAATTATGTTGTTTAATGCAt                                                         >  7:1185055/1‑151 (MQ=255)
                                                                                               gTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTGCTGCCCCCAATATGCTTCAGATGATGAGTGGTGGTGCATTAGGCTATGGCAGGGTGCGGGGATTTTGATGACAGGGGATTTATGGTGTTTAATGCAt                                                         >  7:1288123/1‑151 (MQ=255)
                                                                                               gTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATATCCTTCAGATGCTGGGTGCTCCTCCATTAGGCTATGGGAAGGTGATCAGATTTTCATGACAGGGGAAGTATGATGTTAAATGAAt                                                         >  4:235242/1‑151 (MQ=255)
                                                                                                 tCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCGCCCCTCCTCCGCCCCCAAATATACTTCAGATGATGAGTGATGCTGCATTAGGGTATGGCAAGGTGAGCAGATTTGAATCGCAGGCGAATTATTATTTTTAATGCATgt                                                       >  2:1420343/1‑150 (MQ=255)
                                                                                                 tCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCGGCCCCAAAAATCATTCACATGATGAGTGATGCTGCAGTAGGATATGGCTAGGTGATCAGATTTTAATCACAGGCGATTTATGATGTTAAATTCATgg                                                       >  3:1246756/1‑151 (MQ=255)
                                                                                                  cAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATgg                                                       <  7:699449/150‑1 (MQ=255)
                                                                                                   aggaggAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCCGCCTCCACCCCCAAATATCCTTCAGATGATGAGTGATCCTGCAGTAGGGGATGGCAAGGTGGTCAGAGTTTGATCACAGGGGAAGTATTATGTTAAATGCATGGcc                                                     >  4:345011/1‑150 (MQ=255)
                                                                                                    ggaggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCa                                                     <  6:1347849/150‑1 (MQ=255)
                                                                                                     gaggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCAcc                                                   <  7:91578/151‑1 (MQ=255)
                                                                                                      aggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCGCAGGGGAATTATGATGTTAAATGCATGGCACCt                                                  <  4:521470/151‑1 (MQ=255)
                                                                                                       ggaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTg                                                 <  4:1605276/151‑1 (MQ=255)
                                                                                                         aaaTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGcc                                               <  5:646991/151‑1 (MQ=255)
                                                                                                         aaaTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGcc                                               <  5:448006/151‑1 (MQ=255)
                                                                                                         aaaTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAAAATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGAGTTTCATCAAAGGGGAATTATGATGTTTAATGAATGGCACCTGcc                                               >  7:183743/1‑151 (MQ=255)
                                                                                                           aTCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCgg                                             <  1:1665882/151‑1 (MQ=255)
                                                                                                            tCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCGCCCCCCAATATCCTTCAGCTGATGAGGGATCCTGCATTAGGCGATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGACTGGCAccggccggg                                            >  3:746789/1‑150 (MQ=255)
                                                                                                             ccTTATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTTATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGt                                            <  8:257901/150‑1 (MQ=255)
                                                                                                             ccTGATCCGGCTATTCCGAAAGTTATCCGCGCCCGGCCGCCTCCCCCAAATATACTTCAGATGATGAGTGATGGTGCATTTGGCTATGGGATGGTGATGAGATTTTGATAAGAGGGGATTTATGATGTTTAATGGATGGCCGCTGGCGGTg                                           >  6:1659066/1‑151 (MQ=255)
                                                                                                             ccTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGT‑ATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGtt                                           <  1:675064/150‑2 (MQ=255)
                                                                                                             ccTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCCGCCTCCCCCAAATATCCTTGAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTGATCACAGGGGTATTATGATTTTTAATTCATGGCCCCTGCCgggg                                           >  4:1010955/1‑149 (MQ=255)
                                                                                                               tGATCCGGCTATTCCGAAAGTTATCCGCCCCCGGCGGCAGCCCCCAAAAATAAATCAGATGATGAGTGATCATGCAATAGGGTATTGGAAGGTGATAAGATTTTAATAACAGGCGAATTATAATGTTAAATTAATGGCACCAGCCGGTGcc                                         >  2:196433/1‑151 (MQ=255)
                                                                                                                 aTCCGGCTATTCCGAAAGTTATCCGGGGGCGGGCGGCGGCCCCAACTAGCCTTCAGATGATGAGTGATGCGGCATTAGGGTATGGGAAGGTGATCTGATTTTGATGAGAGGGGAAGTATGATGTTAATTGAATGGCAGATGCGGGTgagcc                                       >  1:831671/1‑147 (MQ=255)
                                                                                                                 aTCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGcccc                                       <  1:1217501/151‑1 (MQ=255)
                                                                                                                     ggCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCAtt                                   <  4:871929/151‑1 (MQ=255)
                                                                                                                     ggCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGACTTTCATCACAGGGGAATTATGACGTTAAATGCATGGCACCTGCCGGTGCCCCCAtt                                   <  5:794887/151‑1 (MQ=255)
                                                                                                                      gCTATTCCGAAAGTTATCCGCCCCCGGCCTCCTCCCCCAAATATACTTCAGATGATGAGTGATGCTGCATTATGCTATGGCAAGGGGCTCAGATTTTGATCACACGGGACTTATGATGTTAAATGAATGGCACCTGCCGGTGGCCCCAttt                                  >  4:286884/1‑151 (MQ=255)
                                                                                                                         aTTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTg                                 <  4:1460500/149‑1 (MQ=255)
                                                                                                                         aTTCCGAAAGTTATCCGCCCCCGGCCGCCGCCCCCAAATAACCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGTTCAGATTTTCATAAAAGGGGAATTATGATGTTAAATGAATGGAACATCCCGGTGCCCCCAgttgtt                               >  7:1626841/1‑151 (MQ=255)
                                                                                                                             cGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATAATTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTAATCACAGGGGAATTATTATTTTAAATTCATGGCACCTGCCGGTGCCCCCATTTGTTaaa                            >  5:1239302/1‑150 (MQ=255)
                                                                                                                              gAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAca                          <  4:1336309/151‑1 (MQ=255)
                                                                                                                              gAAAGTTATCCGCCCCCGGCCGCCTCCCCCAAAAATACTTCAGATGATGAGTGAGCCGGCAGTAGGCTATGGCAAGGTGATCAGAGTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGGGGCCCCCGTTGTGAAAc                           >  1:1269338/1‑150 (MQ=255)
                                                                                                                               aaaGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaa                         <  2:1476332/151‑1 (MQ=255)
                                                                                                                                aaGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaa                        <  7:537327/151‑1 (MQ=255)
                                                                                                                                       ccGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGa                 <  8:1146585/151‑1 (MQ=255)
                                                                                                                                       ccGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGa                 <  6:105364/151‑1 (MQ=255)
                                                                                                                                         gCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCCTGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagcaa                    >  1:1021354/1‑146 (MQ=255)
                                                                                                                                          cccccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCa              <  6:294133/151‑1 (MQ=255)
                                                                                                                                           ccccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCaa             <  3:1133345/151‑1 (MQ=255)
                                                                                                                                           ccccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCa              >  2:546821/1‑150 (MQ=255)
                                                                                                                                                tcctcctccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCt        >  6:127307/1‑151 (MQ=255)
                                                                                                                                                gcctcctccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCt        <  3:848329/150‑1 (MQ=255)
                                                                                                                                                  ctcctccCCCAAATATCCTTCAGATGAAGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTc       >  1:91352/1‑150 (MQ=255)
                                                                                                                                                    cctccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCAtt    >  5:1609403/1‑151 (MQ=255)
                                                                                                                                                    cctccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGACGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCa      <  6:474006/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GCGTGGATTGTTGTGCTGTTGATTGGCTGGATGTTTAAACGCCGCCACGATCGTCAGCTGGCTGAAAACCAGTAATCCTTGCGTTCTGTAAGCCGGTTCAGGAGGAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTAC  >  NC_000913/422285‑422585

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