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

TATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACG  >  NC_000913/540409‑540705
                                                                                                                                                      |                                                                                                                                                    
tATCACGAAAACCATCTTACGTGGTGATGTGATTTACTACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTGATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGgtg                                                                                                                                                     <  1:336000/151‑1 (MQ=255)
     cGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTT‑AGGTGTGTTTa                                                                                                                                               <  2:1435624/151‑1 (MQ=255)
       aaaaCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTa                                                                                                                                               <  1:1138064/150‑1 (MQ=255)
       aaaaCCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTa                                                                                                                                               <  2:1145014/150‑1 (MQ=255)
       aaaaCCATCTTACGTGGTGATGTGATTTACGACATTGAACAAGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTC‑CCGGTTACGGGGGGGTTAt                                                                                                                                              >  2:1075756/1‑151 (MQ=255)
        aaaCCCTCTTTCGGGGGGATGTGATTTACGGCATTGAACAGGGCTTTCCTGTTGGGGCGAAAGGGCAAATTTTCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATg                                                                                                                                             <  2:224803/151‑1 (MQ=255)
           ccATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGCTTCT‑CCGGTTAAGGGGTGGTTATGttt                                                                                                                                          >  1:546674/1‑150 (MQ=255)
           ccATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTc                                                                                                                                          <  2:1295439/151‑1 (MQ=255)
             aTCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCCGTTAAGGGGTGTTTTTGTTCaa                                                                                                                                        >  2:472973/1‑151 (MQ=255)
               cTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTTTCCCGCTTAAGGGGGGttttgtttacat                                                                                                                                       >  1:806643/1‑142 (MQ=255)
               cTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGCTCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAAtt                                                                                                                                      >  1:1179662/1‑151 (MQ=255)
                ttACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAAttt                                                                                                                                     <  2:33874/151‑1 (MQ=255)
                 tACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCCGTTAAGGGGGGGTTTTTTTCAAtttt                                                                                                                                    >  2:961303/1‑151 (MQ=255)
                  aCGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATg                                                                                                                                             <  1:1384090/141‑1 (MQ=255)
                  aCGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGGTAAGGGGGGTTTATgctggctcttt                                                                                                                                   >  2:1384778/1‑141 (MQ=255)
                   cGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGc                                                                                                                                  <  1:729874/151‑1 (MQ=255)
                   cGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGc                                                                                                                                  <  2:569292/151‑1 (MQ=255)
                     tggtgATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAg                                                                                                                                <  2:243188/151‑1 (MQ=255)
                       gtgatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTT‑CCGGTTAAGGGGGGTTTATTTTCAAtttttccgg                                                                                                                               >  2:565180/1‑145 (MQ=255)
                       gtgatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTT‑CCGGGTAAGGGGGGGTTTTTTTCAATTTTTCAGt                                                                                                                               >  2:525887/1‑150 (MQ=255)
                         gatgtgatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGATTCT‑CCCGGTATGGGGGGTTTATGTTTTAttttttaggtag                                                                                                                            >  1:959296/1‑143 (MQ=255)
                              gatTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGGTTAGGGGTGGTTTTGGTCAATTTTGCAGTCA‑GCcgcg                                                                                                                      >  2:51870/1‑151 (MQ=255)
                               atTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCCGGTTAGGGGGGGTTT‑TGTTCAATTTTGCTGCCA‑GCcgcca                                                                                                                     >  1:493722/1‑149 (MQ=255)
                                    cGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGcc                                                                                                                <  2:1486175/151‑1 (MQ=255)
                                      aCATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCT‑CCGGTTAAGGGGTGGTTTTGTTCAATTTTGCAGTCA‑GCCGCCAAAGc                                                                                                                 >  1:554847/1‑148 (MQ=255)
                                        aTTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTTAGGGGTGTTTTTGTTCAATTTTGCCGTCA‑GCCGCGAAAGCCTGtt                                                                                                            >  2:220169/1‑151 (MQ=255)
                                        aTTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTTTGTTCAATTTTGCCGGCA‑GCCCCGAAAGCCTGtt                                                                                                            >  2:464313/1‑151 (MQ=255)
                                           gAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCTGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGTGATTCT‑TCTGCTTAGGTGTGGTTTTGTTTTTTTTTTCAGCCA‑GCCGCGcaattccgcttt                                                                                                          >  1:650346/1‑141 (MQ=255)
                                           gAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAAGTGTGTTTATGTTCAATTTTTCAGTCA‑GCCGCGAAAGCCTGTTAt                                                                                                          >  1:666234/1‑150 (MQ=255)
                                           gAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGGTAAGGGGGGTTTTTGTTCCATTTTGCAGTCA‑GCCGCGAAAGCCTGTTAt                                                                                                          >  2:315956/1‑150 (MQ=255)
                                                gggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCT‑CCGGGTAAGGTGGGTTTTTGTTCTCTTTTGCAGTCA‑GCCGCGAAAGCCTGTTTTCAgggt                                                                                                    >  1:435458/1‑149 (MQ=255)
                                                gggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTTAATTTTGTAGTCA‑GCCGCGCAAGCCTGTTATCAGGAt                                                                                                    >  1:55637/1‑151 (MQ=255)
                                                gggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGGGTGTTTATGTTTCATTTTTCAGTCA‑GCCGCGAAAGCCTGTTTTCAgggt                                                                                                    >  2:67653/1‑149 (MQ=255)
                                                gggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGGTAAGGGGTGTTTTTGTTTAATTTTTCAGTCA‑GGCCCGCAAGGCCGTTATCAGGAt                                                                                                    >  1:705916/1‑151 (MQ=255)
                                                gggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGGTAAGGGGTGGTTTTGTTTAATTTTTCAGTTA‑GCCGCGGAAGCCCGTTAACAgggt                                                                                                    >  2:1186137/1‑149 (MQ=255)
                                                 ggCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAgg                                                                                                      <  2:1080752/148‑1 (MQ=255)
                                                     tCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAg                                                                                               <  2:355654/151‑1 (MQ=255)
                                                      cccTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGt                                                                                              <  2:1330692/151‑1 (MQ=255)
                                                        cTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTgg                                                                                            <  2:949444/151‑1 (MQ=255)
                                                                ccGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCt                                                                                     <  2:1476070/150‑1 (MQ=255)
                                                                      ggTCAATTTATCCTTAAACATCCGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCAGACAGTCT‑CCCTGTAACGTGTGTTTGTTATGCATTTTTCATTCACGCCG‑AACAGCTGGTTATCAGGGTCTCACTGGCTTTTTTTTCtttt                                                                              >  1:958265/1‑151 (MQ=255)
                                                                      ggTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGATTTTCCCGGTTAAGG‑GTGGTTATGTTTAAATTTTCAGTAA‑GCCGCGCAAACCCGGTATCAGGGTTTCAGGGGtttttttttttttt                                                                              >  1:1292223/1‑151 (MQ=255)
                                                                      ggTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTT‑CCGGGTAAGGGGGGTTTTTGTTCAATTTTTCAGTCA‑GCCCCGGAAGACCGTTTTTAGGATTTTCGTGGtttttttttttttt                                                                              >  2:480144/1‑151 (MQ=255)
                                                                      ggTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGGGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTTTCAGGGTTTTAGTGGGTtttttttttttt                                                                              >  2:189299/1‑151 (MQ=255)
                                                                      ggTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGGGTGGTTTTGTTCAATTTTTCAGTCA‑GCCCCGAAAGCCTGTTTTTAGGATTTCAGTGGGTtttttttttttt                                                                              >  2:1507152/1‑151 (MQ=255)
                                                                      ggTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGGGGGTTTATTTTTAATTTTTCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGGGGGTtttttttttttt                                                                              >  2:1090615/1‑151 (MQ=255)
                                                                       gTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTAttt                                                                               <  2:919038/149‑1 (MQ=255)
                                                                         cAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGGTAAGGGGGGTTTTTGTTTAATTTTTCAGGCA‑GCCGCGAAAGCCCGTTATTAGGATTTCAGTGGGTTTTTTTTTTTTTTTg                                                                           >  2:436762/1‑151 (MQ=255)
                                                                           aTTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTTAGGGGGGTTTTTGTTCAATTTTTCCGTCA‑GCCCCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTTTTTATTTTTTGCa                                                                         >  2:1502640/1‑151 (MQ=255)
                                                                           aTTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCAATTT‑TCGGTTAAGGGGTGTTTTTGTTTTAATTTGTCGTCA‑GCCGCGGAAGCCTGTTATCAGGATTTTAGGGGTTTTTTTTTTTTTTTTGCa                                                                         >  1:879871/1‑151 (MQ=255)
                                                                                 ccTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGATTCT‑TCGGTTAAGGGGGGGTTATGGTCAATTTTGCTGTCG‑CCCCCGAAAGCCTGTTGTTAGGATTTCCGCGGTTTTTCTTTATTTTCTGCaaacccg                                                                   >  1:1067596/1‑146 (MQ=255)
                                                                                 ccTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACgg                                                                   <  1:1122473/151‑1 (MQ=255)
                                                                                 ccTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGGGTTTTTGTTCCATTTTGCAGTCC‑GCCGCGGAAGCCCGTTTTCAGGGTTTTAGTGGTTTTTCTTTTTTTTTTTCAAAACgg                                                                   >  1:603535/1‑151 (MQ=255)
                                                                                      aaCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTTAGGGGGGTTTATGTTCAATTTTGCAGTCA‑GCCCCGAAAGCCTGTTTTCAGGATTTTAGGGGTTTTTTTTTTTTTTTTTCAACACGGGtgtg                                                              >  2:1151191/1‑151 (MQ=255)
                                                                                       aCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTgg                                                             <  1:809536/151‑1 (MQ=255)
                                                                                       aCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTgg                                                             <  2:1363844/151‑1 (MQ=255)
                                                                                       aCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGGTAAGGGGTGTTTTTGTTCAATTTTTCAGGCA‑GCCGCGAAAGCCTGTTTTTAGGATTTTAGTGGGTTTTTTTTTTTTTTTTCAACACGGGTGTgg                                                             >  2:787709/1‑151 (MQ=255)
                                                                                        cATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTT‑CCGGTTAAGGTGTGTTTATGTTTAATTTTTCAGTCA‑GCCGCGCAAGCCTGTTATCAGGATTTCAGGGGTTTTTCTTTTTTTTTTTCAACACGGTTGTGGt                                                            >  2:1390695/1‑151 (MQ=255)
                                                                                        cATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGt                                                            <  1:395337/151‑1 (MQ=255)
                                                                                         aTCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGtt                                                           <  2:807102/151‑1 (MQ=255)
                                                                                         aTCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGtt                                                           <  2:138386/151‑1 (MQ=255)
                                                                                         aTCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGtt                                                           <  2:1195631/151‑1 (MQ=255)
                                                                                         aTCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGtt                                                           <  1:182856/151‑1 (MQ=255)
                                                                                          tCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTtc                                                          <  1:953388/151‑1 (MQ=255)
                                                                                           cagcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGGTTAGGGGTGTTTTTGTTTAATTTTTCAGTCA‑GCCGCGAAAGCCTGTTATCAGGGTTTCAGTGGTTTTTCTTTATTTTTTTCAACACGGGTGTGGTtc                                                          >  1:308525/1‑150 (MQ=255)
                                                                                             gcagTAATTTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGGGGGTTTTTGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTTTCAGGGTTTCTGTGGTTTTTTTTTTTTTTTTTCCACACGGTTGTGGGtcct                                                        >  2:1297149/1‑150 (MQ=255)
                                                                                             gcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTTT‑CCGGTTTAGGGGTGGTTTTGTTCAATTTTTTCGTCA‑GCCGCGAAAGCCTGTTATCAGGGTTTTAGGGGTTTTTTTTTTTTTTTTTCAACACGGTTGTGGGtcct                                                        >  2:1068201/1‑150 (MQ=255)
                                                                                             gcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCCTTCT‑CCGGGTTAGGTGTGGTTATGTTCAATTTTTCAGGCA‑GCCGCGAAAGCCCGTTATCAGGGTTTCTGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTttct                                                        >  1:1470496/1‑150 (MQ=255)
                                                                                             gcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTT‑CCCGTTATGGGGGGGTTTTGTTCAATTTTGCAGGCA‑GGCCCCGAAGGCTGTTATCCGGATTTTTGTGGGTTTTCTTTTTTTTTTTCCACACGGTTTTGGGtccc                                                        >  1:1297357/1‑149 (MQ=255)
                                                                                             gcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑TCGGGTAAGGGGTGGTTATGTTCAATTTTGCCGTCC‑GCCGCGAAAGCCTGTTATCAGGGGTTTAGTGGGTTTTCTTTTTTTTTTGTAACACGGGTGTGGTtcct                                                        >  1:1396913/1‑150 (MQ=255)
                                                                                             gcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGGTAAGGGGGGTTTATGTTCAATTTTTCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCTGTGGGTTTTCTTTTTTTTTTTCAACACGGTTTTGGGtcct                                                        >  2:1418743/1‑150 (MQ=255)
                                                                                                  aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGGATTTC‑CCGGGTAAGGGGGGGTTATTTTCAATTTTGCAGTCA‑GCCGGGAAAACCTGTTATCAGGATTTTAGTGGGTTTTTTTTTTTTTTTGCAACACGGTTGTGGTTCCTCCTACg                                                  >  2:505850/1‑151 (MQ=255)
                                                                                                  aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTT‑CCGGTTTAGGGGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTTTCAGGATTTCAGTGGGTTTTTTTTTTTTTTTTCAACACGGGTGTGGGTCCTCCTACg                                                  >  2:1451561/1‑151 (MQ=255)
                                                                                                  aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGGGTGTTTATGTTCAATTTTTCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGGTTTTTTTTTTTTTTTTCAACACGGTTGTGGGTTCTCCTACg                                                  >  2:342433/1‑151 (MQ=255)
                                                                                                  aaTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGGTAAGGGGTGTTTATGTTCAATTTTTCAGTCA‑GCCGCGAAAGCCTGTTATCAGGGTTTTAGTGGGTTTTTTTTTTTTTTTGCAACACGGTTGTGGGTCCTCCTACg                                                  >  2:1487105/1‑151 (MQ=255)
                                                                                                   aTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGc                                                 <  1:1053661/151‑1 (MQ=255)
                                                                                                   aTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGc                                                 <  1:749336/151‑1 (MQ=255)
                                                                                                   aTCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTGTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGc                                                 <  1:1390005/151‑1 (MQ=255)
                                                                                                    tCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTT‑CCGGTTTAGGGGTGTTTTTTTTCAATTTTTCAGTCA‑GCCGCGAAAGCCTGTTTTCAGGGTTTCAGTGGTTTTTTTTTTTTTTTTTCAACACGGTTGTGGTTTCTCCTACg                                                  >  2:631243/1‑149 (MQ=255)
                                                                                                         ccccTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGGGTTTATGTTCCATTTTTCAGGCA‑GCCGCCAAAGCCTGTTCTTAGGGTTTTCGGGGTTTTTCTTTTTTTTTTGCAACACGGGTGTGGGTCCTCCTTCCCTACttt                                           >  1:12427/1‑151 (MQ=255)
                                                                                                         ccccTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGGTAAGGGGTGTTTTTGTTCAATTTTTCAGTCA‑GCCGCGAAAGCCTGTTTTCAGGGTTTCAGTGGTTTTTTTTTTTTTTTTTCAACCCGGGTGTGGTTCCTCCTACGCTACttt                                           >  2:672431/1‑151 (MQ=255)
                                                                                                          cccTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTc                                          <  2:419839/151‑1 (MQ=255)
                                                                                                             tGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGc                                       <  1:85844/151‑1 (MQ=255)
                                                                                                              gCAATGCCCGTCCTTGCGGCGGGCCTTCT‑CCGGGTAAGGGGGGTTTTTGTTTAATTTTGCAGGCA‑GCCCCGAAAACCTCTTTTTCGGGTTTTCGTGGTTTTTTTTTCTTTTTTGCCACACGGTTGTTGTTCCTCCTTCGCttctttctt                                        >  1:1437107/1‑148 (MQ=255)
                                                                                                              gCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGGGTGTTTATGTTCCATTTTGCAGTCA‑GCCGCGAAAGCCTGTTTTCAGGGTTTTAGTGGTTTTTCTTTTTTTTTTTCAACACGGTTGTGGTTCCTCCTCCGCTCCTTTCTg                                        >  2:748779/1‑149 (MQ=255)
                                                                                                                        tCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGc                            <  1:1389965/151‑1 (MQ=255)
                                                                                                                                 cggGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGt                    <  2:254776/150‑1 (MQ=255)
                                                                                                                                 cgcgCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTACAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGt                    <  1:384738/147‑1 (MQ=255)
                                                                                                                                    gCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGc                 <  1:824647/150‑1 (MQ=255)
                                                                                                                                               gTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACg        >  1:95924/1‑149 (MQ=255)
                                                                                                                                               gTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGc       <  2:459148/150‑1 (MQ=255)
                                                                                                                                                ttAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCt      <  2:1005807/150‑1 (MQ=255)
                                                                                                                                                ttAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCt      <  1:525683/150‑1 (MQ=255)
                                                                                                                                                   aGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACg  >  2:247733/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TATCACGAAAACCATCTTACGTGGTGATGTGATTTACGACATTGAACAGGGCTTCCCTGTTGCGCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCT‑CCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCA‑GCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACACGGTTGTGGTTCCTCCTACGCTACTTTCTGCTTTTCAGTTGCCGCAAAGTAGCCTGCTTACGCTCACG  >  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: ‑