Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913540,5580TG44.6% 109.1 / 46.0 93intergenic (+50/‑7)allB/ybbYallantoinase/putative uracil/xanthine transporter
Reads supporting (aligned to +/- strand):  ref base T (13/38);  new base G (41/0);  total (55/38)
Fisher's exact test for biased strand distribution p-value = 5.54e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACG  >  NC_000913/540409‑540705
                                                                                                                                                      |                                                                                                                                                    
tATCACGAAAACCATCTTACGTGGTGATGTGATTTACTACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTGATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgtg                                                                                                                                                     <  1:336000/151‑1 (MQ=255)
     cGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAGGTGTGTTTa                                                                                                                                               <  2:1435624/151‑1 (MQ=255)
       aaaaCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTa                                                                                                                                               <  1:1138064/150‑1 (MQ=255)
       aaaaCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTa                                                                                                                                               <  2:1145014/150‑1 (MQ=255)
       aaaaCCATCTTACGTGGTGATGTGATTTACGACATTGAACAAGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTCCCGGTTACGGGGGGGTTAt                                                                                                                                              >  2:1075756/1‑151 (MQ=255)
        aaaCCCTCTTTCGGGGGGATGTGATTTACGGCATTGAACAGGGCTTTCCTGTTGGGGCGAAAGGGCAAATTTTCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATg                                                                                                                                             <  2:224803/151‑1 (MQ=255)
           ccATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGCTTCTCCGGTTAAGGGGTGGTTATGttt                                                                                                                                          >  1:546674/1‑150 (MQ=255)
           ccATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTc                                                                                                                                          <  2:1295439/151‑1 (MQ=255)
             aTCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCCGTTAAGGGGTGTTTTTGTTCaa                                                                                                                                        >  2:472973/1‑151 (MQ=255)
               cTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTTTCCCGCTTAAGGGGGGttttgtttacat                                                                                                                                       >  1:806643/1‑142 (MQ=255)
               cTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGCTCATTCTCCGGTTAAGGTGTGTTTATGTTCAAtt                                                                                                                                      >  1:1179662/1‑151 (MQ=255)
                ttACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAAttt                                                                                                                                     <  2:33874/151‑1 (MQ=255)
                 tACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCCGTTAAGGGGGGGTTTTTTTCAAtttt                                                                                                                                    >  2:961303/1‑151 (MQ=255)
                  aCGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATg                                                                                                                                             <  1:1384090/141‑1 (MQ=255)
                  aCGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGGGTTTATgctggctcttt                                                                                                                                   >  2:1384778/1‑141 (MQ=255)
                   cGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGc                                                                                                                                  <  1:729874/151‑1 (MQ=255)
                   cGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGc                                                                                                                                  <  2:569292/151‑1 (MQ=255)
                     tggtgATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAg                                                                                                                                <  2:243188/151‑1 (MQ=255)
                       gtgatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTAAGGGGGGTTTATTTTCAAtttttccgg                                                                                                                               >  2:565180/1‑145 (MQ=255)
                       gtgatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGGTAAGGGGGGGTTTTTTTCAATTTTTCAGt                                                                                                                               >  2:525887/1‑150 (MQ=255)
                         gatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGATTCTCCCGGTATGGGGGGTTTATGTTTTAttttttaggtag                                                                                                                            >  1:959296/1‑143 (MQ=255)
                              gatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTTAGGGGTGGTTTTGGTCAATTTTGCAGTCAGCcgcg                                                                                                                      >  2:51870/1‑151 (MQ=255)
                               atTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCCGGTTAGGGGGGGTTTTGTTCAATTTTGCTGCCAGCcgcca                                                                                                                     >  1:493722/1‑149 (MQ=255)
                                    cGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGcc                                                                                                                <  2:1486175/151‑1 (MQ=255)
                                      aCATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTCCGGTTAAGGGGTGGTTTTGTTCAATTTTGCAGTCAGCCGCCAAAGc                                                                                                                 >  1:554847/1‑148 (MQ=255)
                                        aTTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTTAGGGGTGTTTTTGTTCAATTTTGCCGTCAGCCGCGAAAGCCTGtt                                                                                                            >  2:220169/1‑151 (MQ=255)
                                        aTTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTTTGTTCAATTTTGCCGGCAGCCCCGAAAGCCTGtt                                                                                                            >  2:464313/1‑151 (MQ=255)
                                           gAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCTGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGTGATTCTTCTGCTTAGGTGTGGTTTTGTTTTTTTTTTCAGCCAGCCGCGcaattccgcttt                                                                                                          >  1:650346/1‑141 (MQ=255)
                                           gAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAAGTGTGTTTATGTTCAATTTTTCAGTCAGCCGCGAAAGCCTGTTAt                                                                                                          >  1:666234/1‑150 (MQ=255)
                                           gAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGGGTTTTTGTTCCATTTTGCAGTCAGCCGCGAAAGCCTGTTAt                                                                                                          >  2:315956/1‑150 (MQ=255)
                                                gggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTCCGGGTAAGGTGGGTTTTTGTTCTCTTTTGCAGTCAGCCGCGAAAGCCTGTTTTCAgggt                                                                                                    >  1:435458/1‑149 (MQ=255)
                                                gggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTTAATTTTGTAGTCAGCCGCGCAAGCCTGTTATCAGGAt                                                                                                    >  1:55637/1‑151 (MQ=255)
                                                gggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTTCATTTTTCAGTCAGCCGCGAAAGCCTGTTTTCAgggt                                                                                                    >  2:67653/1‑149 (MQ=255)
                                                gggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGTGTTTTTGTTTAATTTTTCAGTCAGGCCCGCAAGGCCGTTATCAGGAt                                                                                                    >  1:705916/1‑151 (MQ=255)
                                                gggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGTGGTTTTGTTTAATTTTTCAGTTAGCCGCGGAAGCCCGTTAACAgggt                                                                                                    >  2:1186137/1‑149 (MQ=255)
                                                 ggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAgg                                                                                                      <  2:1080752/148‑1 (MQ=255)
                                                     tCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAg                                                                                               <  2:355654/151‑1 (MQ=255)
                                                      cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGt                                                                                              <  2:1330692/151‑1 (MQ=255)
                                                        cTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTgg                                                                                            <  2:949444/151‑1 (MQ=255)
                                                                ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCt                                                                                     <  2:1476070/150‑1 (MQ=255)
                                                                      ggTCAATTTATCCTTAAACATCCGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCAGACAGTCTCCCTGTAACGTGTGTTTGTTATGCATTTTTCATTCACGCCGAACAGCTGGTTATCAGGGTCTCACTGGCTTTTTTTTCtttt                                                                              >  1:958265/1‑151 (MQ=255)
                                                                      ggTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGATTTTCCCGGTTAAGGGTGGTTATGTTTAAATTTTCAGTAAGCCGCGCAAACCCGGTATCAGGGTTTCAGGGGtttttttttttttt                                                                              >  1:1292223/1‑151 (MQ=255)
                                                                      ggTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGGTAAGGGGGGTTTTTGTTCAATTTTTCAGTCAGCCCCGGAAGACCGTTTTTAGGATTTTCGTGGtttttttttttttt                                                                              >  2:480144/1‑151 (MQ=255)
                                                                      ggTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTTTCAGGGTTTTAGTGGGTtttttttttttt                                                                              >  2:189299/1‑151 (MQ=255)
                                                                      ggTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGGTTTTGTTCAATTTTTCAGTCAGCCCCGAAAGCCTGTTTTTAGGATTTCAGTGGGTtttttttttttt                                                                              >  2:1507152/1‑151 (MQ=255)
                                                                      ggTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGGGTTTATTTTTAATTTTTCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGGGGGTtttttttttttt                                                                              >  2:1090615/1‑151 (MQ=255)
                                                                       gTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTAttt                                                                               <  2:919038/149‑1 (MQ=255)
                                                                         cAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGGGTTTTTGTTTAATTTTTCAGGCAGCCGCGAAAGCCCGTTATTAGGATTTCAGTGGGTTTTTTTTTTTTTTTg                                                                           >  2:436762/1‑151 (MQ=255)
                                                                           aTTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTTAGGGGGGTTTTTGTTCAATTTTTCCGTCAGCCCCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTTTTTATTTTTTGCa                                                                         >  2:1502640/1‑151 (MQ=255)
                                                                           aTTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCAATTTTCGGTTAAGGGGTGTTTTTGTTTTAATTTGTCGTCAGCCGCGGAAGCCTGTTATCAGGATTTTAGGGGTTTTTTTTTTTTTTTTGCa                                                                         >  1:879871/1‑151 (MQ=255)
                                                                                 ccTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGATTCTTCGGTTAAGGGGGGGTTATGGTCAATTTTGCTGTCGCCCCCGAAAGCCTGTTGTTAGGATTTCCGCGGTTTTTCTTTATTTTCTGCaaacccg                                                                   >  1:1067596/1‑146 (MQ=255)
                                                                                 ccTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACgg                                                                   <  1:1122473/151‑1 (MQ=255)
                                                                                 ccTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGGGTTTTTGTTCCATTTTGCAGTCCGCCGCGGAAGCCCGTTTTCAGGGTTTTAGTGGTTTTTCTTTTTTTTTTTCAAAACgg                                                                   >  1:603535/1‑151 (MQ=255)
                                                                                      aaCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTTAGGGGGGTTTATGTTCAATTTTGCAGTCAGCCCCGAAAGCCTGTTTTCAGGATTTTAGGGGTTTTTTTTTTTTTTTTTCAACACGGGtgtg                                                              >  2:1151191/1‑151 (MQ=255)
                                                                                       aCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTgg                                                             <  1:809536/151‑1 (MQ=255)
                                                                                       aCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTgg                                                             <  2:1363844/151‑1 (MQ=255)
                                                                                       aCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGTGTTTTTGTTCAATTTTTCAGGCAGCCGCGAAAGCCTGTTTTTAGGATTTTAGTGGGTTTTTTTTTTTTTTTTCAACACGGGTGTgg                                                             >  2:787709/1‑151 (MQ=255)
                                                                                        cATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTAAGGTGTGTTTATGTTTAATTTTTCAGTCAGCCGCGCAAGCCTGTTATCAGGATTTCAGGGGTTTTTCTTTTTTTTTTTCAACACGGTTGTGGt                                                            >  2:1390695/1‑151 (MQ=255)
                                                                                        cATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGt                                                            <  1:395337/151‑1 (MQ=255)
                                                                                         aTCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGtt                                                           <  2:807102/151‑1 (MQ=255)
                                                                                         aTCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGtt                                                           <  2:138386/151‑1 (MQ=255)
                                                                                         aTCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGtt                                                           <  2:1195631/151‑1 (MQ=255)
                                                                                         aTCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGtt                                                           <  1:182856/151‑1 (MQ=255)
                                                                                          tCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTtc                                                          <  1:953388/151‑1 (MQ=255)
                                                                                           cagcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTTAGGGGTGTTTTTGTTTAATTTTTCAGTCAGCCGCGAAAGCCTGTTATCAGGGTTTCAGTGGTTTTTCTTTATTTTTTTCAACACGGGTGTGGTtc                                                          >  1:308525/1‑150 (MQ=255)
                                                                                             gcagTAATTTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGGGTTTTTGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTTTCAGGGTTTCTGTGGTTTTTTTTTTTTTTTTTCCACACGGTTGTGGGtcct                                                        >  2:1297149/1‑150 (MQ=255)
                                                                                             gcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTTTCCGGTTTAGGGGTGGTTTTGTTCAATTTTTTCGTCAGCCGCGAAAGCCTGTTATCAGGGTTTTAGGGGTTTTTTTTTTTTTTTTTCAACACGGTTGTGGGtcct                                                        >  2:1068201/1‑150 (MQ=255)
                                                                                             gcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCTCCGGGTTAGGTGTGGTTATGTTCAATTTTTCAGGCAGCCGCGAAAGCCCGTTATCAGGGTTTCTGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTttct                                                        >  1:1470496/1‑150 (MQ=255)
                                                                                             gcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCCGTTATGGGGGGGTTTTGTTCAATTTTGCAGGCAGGCCCCGAAGGCTGTTATCCGGATTTTTGTGGGTTTTCTTTTTTTTTTTCCACACGGTTTTGGGtccc                                                        >  1:1297357/1‑149 (MQ=255)
                                                                                             gcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTTCGGGTAAGGGGTGGTTATGTTCAATTTTGCCGTCCGCCGCGAAAGCCTGTTATCAGGGGTTTAGTGGGTTTTCTTTTTTTTTTGTAACACGGGTGTGGTtcct                                                        >  1:1396913/1‑150 (MQ=255)
                                                                                             gcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGGGTTTATGTTCAATTTTTCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCTGTGGGTTTTCTTTTTTTTTTTCAACACGGTTTTGGGtcct                                                        >  2:1418743/1‑150 (MQ=255)
                                                                                                  aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGATTTCCCGGGTAAGGGGGGGTTATTTTCAATTTTGCAGTCAGCCGGGAAAACCTGTTATCAGGATTTTAGTGGGTTTTTTTTTTTTTTTGCAACACGGTTGTGGTTCCTCCTACg                                                  >  2:505850/1‑151 (MQ=255)
                                                                                                  aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTTAGGGGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTTTCAGGATTTCAGTGGGTTTTTTTTTTTTTTTTCAACACGGGTGTGGGTCCTCCTACg                                                  >  2:1451561/1‑151 (MQ=255)
                                                                                                  aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTTCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGGTTTTTTTTTTTTTTTTCAACACGGTTGTGGGTTCTCCTACg                                                  >  2:342433/1‑151 (MQ=255)
                                                                                                  aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGTGTTTATGTTCAATTTTTCAGTCAGCCGCGAAAGCCTGTTATCAGGGTTTTAGTGGGTTTTTTTTTTTTTTTGCAACACGGTTGTGGGTCCTCCTACg                                                  >  2:1487105/1‑151 (MQ=255)
                                                                                                   aTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGc                                                 <  1:1053661/151‑1 (MQ=255)
                                                                                                   aTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGc                                                 <  1:749336/151‑1 (MQ=255)
                                                                                                   aTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTGTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGc                                                 <  1:1390005/151‑1 (MQ=255)
                                                                                                    tCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTTAGGGGTGTTTTTTTTCAATTTTTCAGTCAGCCGCGAAAGCCTGTTTTCAGGGTTTCAGTGGTTTTTTTTTTTTTTTTTCAACACGGTTGTGGTTTCTCCTACg                                                  >  2:631243/1‑149 (MQ=255)
                                                                                                         ccccTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGGGTTTATGTTCCATTTTTCAGGCAGCCGCCAAAGCCTGTTCTTAGGGTTTTCGGGGTTTTTCTTTTTTTTTTGCAACACGGGTGTGGGTCCTCCTTCCCTACttt                                           >  1:12427/1‑151 (MQ=255)
                                                                                                         ccccTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGTGTTTTTGTTCAATTTTTCAGTCAGCCGCGAAAGCCTGTTTTCAGGGTTTCAGTGGTTTTTTTTTTTTTTTTTCAACCCGGGTGTGGTTCCTCCTACGCTACttt                                           >  2:672431/1‑151 (MQ=255)
                                                                                                          cccTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTc                                          <  2:419839/151‑1 (MQ=255)
                                                                                                             tGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGc                                       <  1:85844/151‑1 (MQ=255)
                                                                                                              gCAATGCCCGTCCTTGCGGCGGGCCTTCTCCGGGTAAGGGGGGTTTTTGTTTAATTTTGCAGGCAGCCCCGAAAACCTCTTTTTCGGGTTTTCGTGGTTTTTTTTTCTTTTTTGCCACACGGTTGTTGTTCCTCCTTCGCttctttctt                                        >  1:1437107/1‑148 (MQ=255)
                                                                                                              gCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTCCATTTTGCAGTCAGCCGCGAAAGCCTGTTTTCAGGGTTTTAGTGGTTTTTCTTTTTTTTTTTCAACACGGTTGTGGTTCCTCCTCCGCTCCTTTCTg                                        >  2:748779/1‑149 (MQ=255)
                                                                                                                        tCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGc                            <  1:1389965/151‑1 (MQ=255)
                                                                                                                                 cggGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGt                    <  2:254776/150‑1 (MQ=255)
                                                                                                                                 cgcgCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTACAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGt                    <  1:384738/147‑1 (MQ=255)
                                                                                                                                    gCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                 <  1:824647/150‑1 (MQ=255)
                                                                                                                                               gTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACg        >  1:95924/1‑149 (MQ=255)
                                                                                                                                               gTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGc       <  2:459148/150‑1 (MQ=255)
                                                                                                                                                ttAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCt      <  2:1005807/150‑1 (MQ=255)
                                                                                                                                                ttAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCt      <  1:525683/150‑1 (MQ=255)
                                                                                                                                                   aGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACg  >  2:247733/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACG  >  NC_000913/540409‑540705

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: