Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,356,201	0	T	G	41.2%	116.1 / 19.1	68	T304P (ACC→CCC)	sapA	antimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (2/38);  new base G (28/0);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 4.36e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCC  >  NC_000913/1356051‑1356351                                                                                                                                                           |                                                                                                                                                          ccTCGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCTGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑Ggt                                                                                                                                                           <  2:605184/151‑1 (MQ=255)        gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑TGG‑GGGG‑tt                                                                                                                                                       >  2:200171/1‑147 (MQ=255)             aGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCa                                                                                                                                               <  2:432815/150‑1 (MQ=255)             aGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAg                                                                                                                                              <  1:180124/151‑1 (MQ=255)                  cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGCGGGGG‑TTAAAtcc                                                                                                                                                 >  1:164415/1‑142 (MQ=255)                    cGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCg                                                                                                                                      <  1:694134/151‑1 (MQ=255)                    cGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑GTT‑TGG‑GGGG‑TTAAATGcc                                                                                                                                                >  2:71779/1‑142 (MQ=255)                     gAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑TGC‑GGGG‑TTAAATGCCAGATA‑GGCg                                                                                                                                      >  1:384134/1‑150 (MQ=255)                      aGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGAc                                                                                                                                    <  2:1090157/151‑1 (MQ=255)                       ggCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGT‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACg                                                                                                                                   <  1:530133/151‑1 (MQ=255)                         cgcgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGtt                                                                                                                                 <  2:723709/151‑1 (MQ=255)                          gcgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTT‑GGG‑GGGG‑TTaaagacccagaggcgg                                                                                                                                       >  2:831072/1‑133 (MQ=255)                            gcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑TGC‑GGGG‑GTTAAAGCCAGaaacggggaggt                                                                                                                                  >  1:488254/1‑137 (MQ=255)                                 aaaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑TGC‑GGGG‑GTAAAAGCCAGAAA‑AGCGACGGTCATcc                                                                                                                            >  1:1152343/1‑148 (MQ=255)                                 aaaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCa                                                                                                                          <  1:151106/150‑1 (MQ=255)                                 aaaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAg                                                                                                                         <  1:45986/151‑1 (MQ=255)                                 aaaaTAGAGGCCGCCGTTTCAGCCGTACCACAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAg                                                                                                                         <  1:704528/151‑1 (MQ=255)                                  aaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAgg                                                                                                                        <  1:374768/151‑1 (MQ=255)                                   aaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGGG‑TTAAATGCCAAAAA‑AGCGACGGTCCTccccag                                                                                                                        >  1:761941/1‑147 (MQ=255)                                     tAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGGTTATTTAGCGGC‑GG‑CTT‑GTG‑GGGG‑TTTAAAGCCAGATA‑TGAGGCAGTCATCCCCgggc                                                                                                                      >  1:255744/1‑148 (MQ=255)                                      agagGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACg                                                                                                                     <  1:487894/150‑1 (MQ=255)                                        agGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑ATT‑GGG‑GGGG‑TTATAAGCCAGAAA‑GGCGGCGCTTATCCCCAGACg                                                                                                                     >  2:1176361/1‑148 (MQ=255)                                              ccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGGG‑TTTAAAGCCAGATA‑GGCGCCGGTCCTCCCCgg                                                                                                                        >  1:97526/1‑139 (MQ=255)                                               cgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGATTT‑GGG‑GGGG‑TTAAA‑GTCAGAAA‑AGCGAAGTTTATCCCCGGGAGCACCGGTaa                                                                                                            >  2:832628/1‑150 (MQ=255)                                                gTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAAc                                                                                                          <  2:784538/151‑1 (MQ=255)                                                  ttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACg                                                                                                         <  1:109477/150‑1 (MQ=255)                                                    cAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGGCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCa                                                                                                       <  2:668643/150‑1 (MQ=255)                                                      gCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAgc                                                                                                    <  1:433703/151‑1 (MQ=255)                                                        cGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAgcgc                                                                                                  <  1:235037/151‑1 (MQ=255)                                                        cGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑GTT‑TGC‑GGGG‑GTTAAAGTCAAAAA‑AGCGGCGTTTATCCCCGGACCCCGCGGTTAAAGCaa                                                                                                      >  2:647054/1‑147 (MQ=255)                                                           aCCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTT‑GGG‑GGGG‑TTAAA‑GACAGATA‑GGCGACGTTCTTCCCCCGAAGCCGCGGTTAAAGCAAACGgggg                                                                                               >  2:942920/1‑151 (MQ=255)                                                           aCCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑TGC‑GGGG‑TTTAAAGCCCGATA‑AGCGGCGGTCATCCCCGGACGCAGCGGTTAAAGCAAACGCggg                                                                                               >  1:821770/1‑151 (MQ=255)                                                             cATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑TGC‑GGGG‑GTTAAAGCCCGATA‑AGCGGCGGTTCTCCCAAGACGCAGCGTTAAACGCAAACGCggggc                                                                                             >  1:924914/1‑149 (MQ=255)                                                             cATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgtc                                                                                             <  2:79778/151‑1 (MQ=255)                                                               taataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGGGTTTAATTGCCA‑ATA‑GGCGGCGTTATTCCCCGGCGGCAGGGCTaaaaaacaccgcgggg                                                                                              >  2:1053535/1‑135 (MQ=255)                                                                   aGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGGG‑GTTAAAGCCAAAAAGGGGGGGGTTCTCCCCGGGACCCGGGGTTAAAC‑CAAACGCggg                                                                                               >  2:1163361/1‑143 (MQ=255)                                                                           tGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATg                                                                               <  2:707261/150‑1 (MQ=255)                                                                              aTCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATa                                                                           <  2:782970/151‑1 (MQ=255)                                                                                cAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑TGC‑GGGG‑GTTAATGCCCGAAA‑AGCCCCGGTCCTCCCCGGACCCCGCGTTTAACGCCAACGCGGGCCCTCACGCCAAAATggta                                                                            >  2:777987/1‑147 (MQ=255)                                                                                 aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGG‑TTT‑GGG‑GGGG‑TAAAAAG‑CAGATA‑AGCGACGGTTATCCCCCGACGCCGCGTTAAAAGCAAGCGGGGGGTGGCCCC‑CAAAAAGGAt                                                                            >  1:982827/1‑147 (MQ=255)                                                                                   gcgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGGG‑TTTAAAGCCAGAAA‑AGGGACGTTCCTCCCCGGGCGGCGCGGTTAACGCAAACGCGGGGCGGCACG‑CCAAAAGGATAgctgg                                                                      >  1:985395/1‑151 (MQ=255)                                                                                        ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAgctggctg                                                                   <  2:659598/149‑1 (MQ=255)                                                                                         gTTATTAATCGCCAGTGCCAGCACATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGGG‑GGGG‑TTAAAGACCA‑CTA‑GGGGGCGTTCATCCCCCGACGCCGCGGTAAAAGCAAACGCGGGGCGGCACC‑CAAAAAAGAGAGcagggtggg                                                                 >  2:273880/1‑149 (MQ=255)                                                                                                   gCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGGG‑GTTAAAGACAGAAA‑GGGGACGGTTATCCCCGGACCCAGCGGTGAAAGCAAACGGGGGGCGGCACG‑CCAAACGGAAAACagggtggcgg                                                               >  2:851681/1‑140 (MQ=255)                                                                                                     cAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                    <  2:124425/151‑1 (MQ=255)                                                                                                      aGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑TGC‑GGGG‑TTAAAAGCCCGATA‑AGCGGAGGTCCTCCCACGACGCCGCGTTAAAAGCAAGCGCGGGGCGGCCCG‑CCAAAAGGGAAGCTGGCTGGGAGCaggccaggcc                                                    >  1:1176670/1‑150 (MQ=255)                                                                                                          ccAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAAc                                               <  1:655666/151‑1 (MQ=255)                                                                                                           cAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCTAGaa                                                <  1:948048/149‑1 (MQ=255)                                                                                                             gcgcATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTc                                            <  2:1213256/151‑1 (MQ=255)                                                                                                             gcgcATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTc                                            <  2:578997/151‑1 (MQ=255)                                                                                                              cgcATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑TGC‑GGGG‑TTTAAAGCCCGATA‑AGCGACGGTTATTCCAGGGCCCAGCGGTAAAAGCAAGCGCGGGGCGTCACC‑CAAAAAGGAAAGATGGCTGGGAGCCGGGCCGGGCAGaa                                                >  1:753325/1‑146 (MQ=255)                                                                                                                 aTGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCa                                         <  2:883705/150‑1 (MQ=255)                                                                                                                    gCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTcc                                     <  1:204591/151‑1 (MQ=255)                                                                                                                    gCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTcc                                     <  2:975924/151‑1 (MQ=255)                                                                                                                    gCGGACAGCGGGATTATTTAGCGGC‑GG‑GTT‑GGC‑GGGG‑GTAAAAAGCGAAAA‑AGCGAGGTTCGTCCCCGGGCGCCGGGGTTAAAGGAAGGGCGGGGTCACCCG‑CCAAAAGGGAAGCTGGCTGGGGGGAGGCCAGGCCACAACGACGCaatt                                      >  1:693019/1‑147 (MQ=255)                                                                                                                      ggACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGGG‑TTAAAAGCCAGAGA‑GGCGACGTTCATCCCAGGAGGCAGCGTTAAAAGCAAACGCGGGGCGGCACG‑CAAAATGGAAAAGTGGGTGGGAGCGGGGCAGGCCAGAAAGGCGCATTtccc                                   >  1:316655/1‑151 (MQ=255)                                                                                                                        aCAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCgg                                 <  1:68881/151‑1 (MQ=255)                                                                                                                         cAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGt                                <  1:198481/151‑1 (MQ=255)                                                                                                                          aGCGGGATTATTTAGCGGC‑GG‑TTT‑GGG‑GGGG‑TTAAATACCAGAAA‑AGCGACGTTGAACCCCGGGAGCCGAGGTGAAAGCAAGCGCCGGGCGCCACC‑CAAAAAAGAAATGTGGCTTGGAGCCGGCCCGGCCAGGAAACCGCATTtccc                                   >  2:261884/1‑147 (MQ=255)                                                                                                                          aGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTc                               <  2:742292/151‑1 (MQ=255)                                                                                                                           gCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCa                              <  1:94384/151‑1 (MQ=255)                                                                                                                             gggATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAgg                            <  1:200179/151‑1 (MQ=255)                                                                                                                                      tAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGa                    <  2:607254/150‑1 (MQ=255)                                                                                                                                       aGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACa                  <  2:108028/151‑1 (MQ=255)                                                                                                                                           gc‑gg‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACg              <  1:808478/151‑1 (MQ=255)                                                                                                                                              gg‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTc            >  2:1086840/1‑151 (MQ=255)                                                                                                                                               g‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTcc           <  2:621301/151‑1 (MQ=255)                                                                                                                                                 ttt‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCg          >  1:148086/1‑151 (MQ=255)                                                                                                                                                  tt‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCg          <  2:283955/150‑1 (MQ=255)                                                                                                                                                           tg‑tTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgcc  >  2:682216/1‑151 (MQ=255)                                                                                                                                                           tg‑tTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgcc  >  2:1128145/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                                          CCTCGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTT‑GGC‑GGTG‑TTAAATGCCAGATA‑GGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACG‑CAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCC  >  NC_000913/1356051‑1356351

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑