Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,238,049	0	A	G	41.8%	77.3 / 18.1	55	E139G (GAA→GGA)	yjbH	DUF940 family extracellular polysaccharide protein
Reads supporting (aligned to +/- strand):  ref base A (7/25);  new base G (23/0);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.13e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCAGCGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAA  >  NC_000913/4237903‑4238190                                                                                                                                                      |                                                                                                                                                   cagcagCGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAatat                                                                                                                                                <  2:640735/151‑1 (MQ=255)  agcagcGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGGGCGCGGGATATCGGC‑GGGA‑CGGGGGTGGTTTGTTGGGGaaa                                                                                                                                                  >  2:575254/1‑148 (MQ=255)      gcGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGA‑CGGGGGG‑GGTTTGTGCGGGAAATCtttt                                                                                                                                           >  2:229696/1‑149 (MQ=255)       cGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGGTGGTGCGCGGGGTATCGGCGGGTA‑GGGGGTT‑GTTTGTTGGGTAATTTTttgt                                                                                                                                           >  2:606378/1‑151 (MQ=255)         tCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGA‑CGGGGGT‑GTTTTTTGCGGGATATCttttt                                                                                                                                          >  2:215081/1‑149 (MQ=255)             aGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑c                                                                                                                                   <  1:589254/151‑1 (MQ=255)              gCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑a                                                                                                                                 <  1:940501/151‑1 (MQ=255)                gTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGGTATCGGC‑GGGA‑CGGGGGT‑GTTTTTTTCGGGATATAttttt                                                                                                                                          >  2:980917/1‑142 (MQ=255)                gTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑aaa                                                                                                                               <  1:846534/151‑1 (MQ=255)                gTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGGGGGCGGGATATCTGCGGGTT‑CGGGG‑T‑GTTTTGTGGGGGATAACTTTTTGTCAG‑C‑aa                                                                                                                                >  2:623544/1‑150 (MQ=255)                   ctctGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGA‑CGGGGGG‑GTTTTTTGCGGGATATTTTGTTGTCCG‑C‑aaaaac                                                                                                                            >  2:955633/1‑148 (MQ=255)                    tctGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGA‑CGGGGGC‑GGTTTATTGGGGATATTTTTTTTCCCG‑C‑aaaa                                                                                                                              >  2:896325/1‑147 (MQ=255)                       ggCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑Tgg                                                                                                                        <  2:904183/150‑1 (MQ=255)                               aaCGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGGGCGCGGGATATCGGC‑GGGA‑CGGGGGT‑GGTTTGTGGGGGAAATTTTTGTGGCAG‑C‑AAAACC‑Tggggg                                                                                                                     >  1:131676/1‑145 (MQ=255)                                 cGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGa                                                                                                             <  1:210907/151‑1 (MQ=255)                                 cGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGA‑CGGGGGC‑GTTTGGTGCGGGATATTTTTTTGTCCG‑C‑AAAACC‑CGGGGGGCGTTCg                                                                                                              >  2:97737/1‑150 (MQ=255)                                   tataAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGA‑CGGGGGT‑GTTTTGTGCGGGATATTTTTTTGCCAG‑C‑AAAACC‑CGGGGGGCggttggtt                                                                                                           >  2:874024/1‑144 (MQ=255)                                   tataAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGA‑CGGGGGC‑GGTTGTTGGGGGAAATTTTTTTGGCAGCC‑AAAACC‑GGGGGGCGGttggtt                                                                                                            >  1:126909/1‑147 (MQ=255)                                        agataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGA‑CGGGGGT‑GTTTTATTCGGGATATATTGTTTTCCA‑C‑AAAACC‑CCGGGGCCGGTCGAattt                                                                                                         >  2:360213/1‑148 (MQ=255)                                          ataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCt                                                                                                    <  1:236364/151‑1 (MQ=255)                                          ataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGA‑CGGGGGT‑GTTTGGTGCGGGAAATTTTTTTTTCAG‑CGAAAGCC‑CGGGGG‑CGTTCGGtttt                                                                                                         >  1:21338/1‑146 (MQ=255)                                           taaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGA‑CGGGGGC‑GGTTTTTGGGGGAAATTTTGTTTGCCG‑C‑AAAACC‑TGGGGGGCGGTTGGttttttg                                                                                                      >  2:532755/1‑145 (MQ=255)                                            aaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGGT‑GTTTTATGGGGGATATCTTTTTGCCCG‑C‑CAAGAC‑TGGGGGGCGTTTGGTTTTATGccc                                                                                                   >  1:197454/1‑148 (MQ=255)                                            aaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑g                                                                                                 <  2:752472/151‑1 (MQ=255)                                                 cTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGA‑CGGGGCT‑GTTTTGTGGGGGAAATTTTTTTGGCAG‑G‑AAAGCC‑CGGGGGGCGGTCGATTTTTCGCTC‑gggc                                                                                              >  2:370950/1‑147 (MQ=255)                                                          cAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCctgggc                                                                                         <  1:599459/145‑1 (MQ=255)                                                            aaCTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTa                                                                                <  2:557052/151‑1 (MQ=255)                                                                gagTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGACCGGGGCG‑GTTTG‑TGGGGGATATTTTTTTTCCCG‑C‑CAAGCCAGGGGGGCGGTTTGGTTTTCCGCCC‑GGCCGGGGCg‑ggggg                                                                                  >  1:599163/3‑146 (MQ=255)                                                                gCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTg                                                                            <  1:731871/151‑1 (MQ=255)                                                                   tttGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGGGCGCGGGATATCGGC‑GGTA‑CGGGGGT‑GGTTTTTGTGGGAAATCTTTTTTCCCG‑C‑AAAAAC‑CGGGGGGGCTTCGGTTTTTTGCCC‑GGGCTGGGGTGGGGGGGATTTggg                                                                          >  2:759963/1‑151 (MQ=255)                                                                       tgGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGA‑CGGGGGT‑GTTTGATGGGGGAAATCTTGTTGCCAG‑C‑AAAACC‑TGGGGGCCGTTCGGTTTTTCGCTC‑GGCCCGGGCT‑GGGGGGATTTTGGCAcc                                                                      >  2:335247/1‑150 (MQ=255)                                                                         ggAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGc                                                                   <  1:607485/151‑1 (MQ=255)                                                                          gAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGGTATCGGC‑GGGA‑AGGGGGC‑GTTTTGTGCGGGATATTTTGTTGGCCG‑G‑AAAACC‑TGGGGGGCGGTCGGTTTTTATCCC‑GGCCCGGGGT‑GGGGGGATTTTGGGAc                                                                       >  1:833471/1‑146 (MQ=255)                                                                               gagTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTTATGGGGGAAATCTTGTTGCCAG‑C‑AAAACC‑CGGGGGGCGGTTGGTTTTTAGCTC‑GGCCTGGGGT‑GGGGGGATTTTGGGACCCGCGGGaa                                                              >  2:624283/1‑150 (MQ=255)                                                                                       ggCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGA‑CGGGGCT‑GTTTTATGCGGGATATTTTTTTGGCAG‑G‑AAAAAC‑CGGGGGGCGGTCGATTTTTCGCCC‑GGGCTGGGGT‑GGGGGGATTTGGGCAACAGGGGGAATGTGgaaa                                                      >  2:938870/1‑150 (MQ=255)                                                                                           gCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATc                                                   <  2:108330/149‑1 (MQ=255)                                                                                             cGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑AAGGGGC‑GTTTTAAGGGGGATATTTTGTTGTCAA‑C‑CAAACC‑TGGGGGGCGGTGGATTTTAAGCTCGGGCCCGGGGT‑GGGGGTTTTTTGG‑ACCAGCGGGAATGTGAAAAAACCg                                                 >  1:873814/1‑149 (MQ=255)                                                                                             cGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGT‑CGGGGGC‑GGTTGGTGGGGGATATTTTTTTGGCCG‑C‑AAAAAC‑CGGGGGGCGGTCGGGTTTTCGCGC‑GGCCCCGGGG‑GGGGGGGTTTTGGGACCCGCGGGGAAGTGGAAAAAcc                                                  >  1:277952/1‑148 (MQ=255)                                                                                                aaGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGTGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTg                                            <  1:932271/151‑1 (MQ=255)                                                                                                  gTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGtt                                          <  2:40473/151‑1 (MQ=255)                                                                                                    ggCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCa                                        <  2:771380/151‑1 (MQ=255)                                                                                                    ggCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCa                                        <  1:123413/151‑1 (MQ=255)                                                                                                       ggTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGGTGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGGAATGTGAAAAATCCGCTCTGTTCAGcc                                     >  2:166515/1‑151 (MQ=255)                                                                                                         ttGGTGCGCGGGATATCGGC‑GGTA‑CCGGGGT‑GTTTGATGCGGAATATTTTGTTTCCAG‑C‑CAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGGGAAAAATCCGCTCTGGTCAGCCa                                    >  2:821161/1‑150 (MQ=255)                                                                                                                gcgGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATaaa                            <  2:227645/151‑1 (MQ=255)                                                                                                                   ggATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAatat                         <  1:423362/151‑1 (MQ=255)                                                                                                                   ggATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGGATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAata                          >  2:252351/1‑150 (MQ=255)                                                                                                                    gATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATAtt                        <  1:951850/151‑1 (MQ=255)                                                                                                                          ggc‑ggTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATc                  <  2:749018/151‑1 (MQ=255)                                                                                                                                 a‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATa            <  1:705295/151‑1 (MQ=255)                                                                                                                                 a‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATa            <  2:562086/151‑1 (MQ=255)                                                                                                                                   cGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATaa           <  1:968846/151‑1 (MQ=255)                                                                                                                                   cGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATa            >  1:865449/1‑150 (MQ=255)                                                                                                                                      ggCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAAc          >  2:444052/1‑149 (MQ=255)                                                                                                                                           gTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAAGAGCTAc    >  2:576281/1‑151 (MQ=255)                                                                                                                                              tGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACaa  <  2:86509/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CAGCAGCGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTA‑CGGGGCT‑GTTTGATGCGGAATATCTTGTTGCCAG‑C‑AAAGCC‑TGGGGGCCGTTCGATTTTACGCTC‑GGCCTGGGCT‑GGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAA  >  NC_000913/4237903‑4238190

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