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

CGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGG  >  NC_000913/4237909‑4238199                                                                                                                                                |                                                                                                                                                         cGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGGGCGCGGGGTATCGGCGGTA‑CGGGGGTG‑GTTGATGC‑GGGATATCt‑tgt                                                                                                                                                >  1:195381/1‑150 (MQ=255)     gAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGGTG‑TTTTATGC‑GGGATATCt‑tttttgc                                                                                                                                            >  2:921939/1‑147 (MQ=255)       aGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGGGCGCGGGATATCGGCGGGA‑CGGGGGTG‑GTTGGTGC‑GGGATATTt‑tgtt                                                                                                                                               >  2:143433/1‑145 (MQ=255)         cGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGGGCGCGGGATATCGGCGGGA‑AGGGGGGG‑TTTGATGG‑GGGATATCt‑tttt                                                                                                                                               >  2:31897/1‑143 (MQ=255)          gTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCaaa                                                                                                                                      <  2:56518/151‑1 (MQ=255)          gTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGGTG‑GTTGATGC‑GGGATATTT‑TTTTTGCAGCaa                                                                                                                                       >  2:492600/1‑150 (MQ=255)          gTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGGGCGCGGGATATCGGCGGGA‑CGGGGGGG‑TTTGGTGG‑GGGAAATTt‑tgtttgccg                                                                                                                                          >  2:609847/1‑142 (MQ=255)             ctctGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGGTG‑TTTTGTGC‑GGGATATGT‑TGTTGGCAGCaaa                                                                                                                                      >  2:954761/1‑148 (MQ=255)                  gCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTgggg                                                                                                                             <  2:913416/151‑1 (MQ=255)                  gCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTgggg                                                                                                                             <  1:746496/151‑1 (MQ=255)                     aTCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGcc                                                                                                                          <  2:694239/151‑1 (MQ=255)                      tCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGGTG‑TTTTTTGG‑GGGATATCT‑TGTTTCCCGCAAA‑ACCTggggg                                                                                                                            >  2:461573/1‑148 (MQ=255)                        aaaCGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGGTG‑TTTTTTGC‑GGGAAATTT‑TTTTGTCCGCAAA‑ACCCggggg                                                                                                                            >  2:192757/1‑146 (MQ=255)                          aCGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCg                                                                                                                     <  2:587959/151‑1 (MQ=255)                            gTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGAt                                                                                                                   <  2:302316/151‑1 (MQ=255)                                   gataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGGTATCGGCGGTT‑CCGGGCTG‑CTTGATGC‑GGGGAATTT‑TTTTTCCTCCAAA‑GCCCGCGGGGCGGTCGGtttt                                                                                                                >  1:726393/1‑147 (MQ=255)                                    ataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGGTG‑TTTGATGC‑GGGATATTT‑TGTTGCCCGCAAA‑ACCCGGGGGGCGGTCGGttttttgc                                                                                                            >  2:813229/1‑146 (MQ=255)                                     taaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGGGCGCGGGATATCGGCGGGA‑CCGGGGTG‑GTTTTGGC‑GGGAAATTT‑TTTTTGCCGCAAAGGCCGGGGGGCGGTTGGTtttttcccc                                                                                                           >  2:707162/1‑146 (MQ=255)                                        agCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGGGCGCGGGGGATCGGGGGGA‑CGGGGGCG‑GTTGTTGG‑GGGATATCT‑TTTTGGCCGCCAA‑GCCTGGGGGGCGGTTGGTTTTTCGCCCg                                                                                                         >  1:700639/1‑149 (MQ=255)                                         gCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTTGTGCGCGGGATATCGGGGGGAGCGGGGC‑G‑GTTGTTGC‑GGGAAATCT‑TTTTGTCAGCAAA‑ACCCGGGGGGCGGTCCATTTTACGCTCgg                                                                                                        >  1:313207/1‑149 (MQ=255)                                          ccTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTTGTGGGCGGGGTATCGGCGGGA‑CGGGGGGG‑GTTTCTGC‑GGGGGAATT‑TTTTGCCCGCCAA‑GACTGGGGGGCGGTTGGtttttc                                                                                                              >  1:1006206/1‑140 (MQ=255)                                             tCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCctg                                                                                                    <  2:324830/149‑1 (MQ=255)                                             tCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCctg                                                                                                    <  1:666827/149‑1 (MQ=255)                                              cGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGGTG‑TTTTGTGC‑GGGATATCT‑TGTTGGCCGCAAA‑GCCCGGGGGGCGGTCGGTTTTTTGCCCGGCctggg                                                                                                  >  2:371718/1‑150 (MQ=255)                                                 tctcAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGCTG‑TTTTGTGC‑GGGAAATTT‑TTTTTCCCGCAAA‑ACCCGGGGGGCGGTCGATTTTTCGCTCCGCccggggtgg                                                                                              >  2:957673/1‑151 (MQ=255)                                                        cTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGGTT‑TTTTGTGC‑GGGAAATTT‑TTTTTCCCGCAAA‑ACCTGGGGGGCGGTCCGTTTTAACCCCGGCCTGGGCTGGGGGGAt                                                                                        >  2:228476/1‑150 (MQ=255)                                                            gTTTGTGGGAAGAGAGTGACTGGCTGCCGCAAGTGGCGGTTGGGGCGCGGGATAACGGCGGGA‑CGGGGGGG‑GTTGGGGC‑GGGAAATATATGTGGGCAG‑AAA‑GACAGGGGGGGGGGCGCTTTTACGCTCCGCCTGGGGTgggggggttt                                                                                      >  1:291191/1‑143 (MQ=255)                                                              ttGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGGTG‑GTTTATGC‑GGGATATCT‑TTTTGGCCGCAAA‑AACCGGGGGGCGCTTGGTTTTTCGCCCCGCCTGGGGTGGGGGGATTTtgg                                                                                   >  2:167420/1‑149 (MQ=255)                                                                gtgGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTTGTGCGCGGGAAATCGGCGGTA‑TGGGGCCGTTTTGTTGC‑GG‑ATATCT‑TTTTGTCCGCAAA‑GGCTTGGGGCCGGTCGTTTTTTTCCCTGGCCCGGGCAGGGGGGTTTTTGGCAcc                                                                               >  2:737219/1‑151 (MQ=255)                                                                gtgGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGCTG‑TTTTGTGC‑GGGAAATCT‑TGTTTCCCGCAAA‑ACCTGGGGGGCGTTTGATTTTTCGCTCGGCCCGGGGTGGGGGGATTTTGGCAcc                                                                               >  2:1014634/1‑151 (MQ=255)                                                                  gggAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCa                                                                              <  1:743437/150‑1 (MQ=255)                                                                   ggAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGc                                                                            <  2:357955/151‑1 (MQ=255)                                                                    gAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCg                                                                           <  2:691008/151‑1 (MQ=255)                                                                        agagTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGGTATCGGCGGGG‑CGGGGCCC‑GTTGGAGC‑GGGATATTT‑TTTTTGCAGCAAAGGCCGGGGGGCGGTTGGTTTTTA‑GCTCCGCCCGGGGCGGGGGGATTTggggcacagagg                                                                          >  1:172619/1‑139 (MQ=255)                                                                         gagTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAAt                                                                      <  1:903967/151‑1 (MQ=255)                                                                          agTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑AGGGGCTG‑TTTTGTGC‑GGGATATCT‑TTTTGCCCGCCAA‑ACCTGGGGGGCGGTCGGTTTTTCGCCCGGGCTGGGGTGGGGGGATTTTGGGACCAGCGGGaa                                                                       >  2:181656/1‑149 (MQ=255)                                                                           gTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑AGGGGGTG‑TTTTTTGCGGGAAATTTT‑TTTTGCCGCCAAA‑ACCGGGGGGCGGTT‑GATTTTTTGCTCGGGCTGGGGTGGGGGGGTTTTGGGACCCGCGGGaa                                                                       >  2:456318/1‑148 (MQ=255)                                                                             tACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAAtgtg                                                                   <  1:812158/150‑1 (MQ=255)                                                                              aCTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATGGGCGGGA‑GCGGGGTG‑GTTGATGG‑GGGATATTT‑TGTTGGCAGGGAA‑ACCTCGGGGGCGGGTGGTTTTTCCCACGGGCCCGGGTGGGGGGGGTTGGGCAcccccgggactgt                                                                    >  1:169584/1‑137 (MQ=255)                                                                               cTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCCG‑TTTGGTGC‑GGGATATCT‑TGTTGCCAGCAAA‑ACCTGGGGGCCGGTCGGTTTTTCGCCCGGGCTGGGCTGGGGGGATTTTGGCACCAGCGGGAATGTGaa                                                                 >  1:251969/1‑150 (MQ=255)                                                                                 ggCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGaaa                                                                <  2:897349/149‑1 (MQ=255)                                                                                     gCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATcc                                                           <  2:815215/150‑1 (MQ=255)                                                                                     gCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATc                                                            <  2:53989/149‑1 (MQ=255)                                                                                       cGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGct                                                        <  2:930454/151‑1 (MQ=255)                                                                                       cGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGGTG‑GTTTATGC‑GGGATATCT‑TTTTTGCCGCAAA‑AACTGGGGGGCGGTTGATTTTTCGCTCGGGCTGGGGTGGGGGGATTTTGGCACCCGCGGGAAAGTGAAAAAAcccc                                                         >  2:475427/1‑148 (MQ=255)                                                                                              ggCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTc                                                  <  2:409360/150‑1 (MQ=255)                                                                                              ggCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTc                                                  <  2:423445/150‑1 (MQ=255)                                                                                                 ggTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGcc                                              <  1:1007863/151‑1 (MQ=255)                                                                                                  gTTGGTGCGCGGGATATCGGCGGTA‑CGGGGGTG‑TTTGATGC‑GGGATATCT‑TGTTGCCAGCCAA‑GCCTGGGGGCCGTTCGATTTTACGCCCGGCCTGGGCTGGGGGTATTTGGGCACCCGCGGTAATGTGAAAAATCCGCTCTGTTCAGcc                                              >  1:754834/1‑150 (MQ=255)                                                                                                  gTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCa                                             <  1:486901/151‑1 (MQ=255)                                                                                                     ggTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGt                                           <  2:236212/150‑1 (MQ=255)                                                                                                       tGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCCAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGGATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGAt                                        >  2:330209/1‑151 (MQ=255)                                                                                                        gcgcGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGAt                                        <  2:707087/150‑1 (MQ=255)                                                                                                        gcgcGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGAt                                        <  2:502114/150‑1 (MQ=255)                                                                                                            gggATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAata                                   <  2:390873/151‑1 (MQ=255)                                                                                                                    ggcggTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATc                           >  1:529376/1‑151 (MQ=255)                                                                                                                    ggcggTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATc                           >  2:63704/1‑151 (MQ=255)                                                                                                                     gcggTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATcg                          <  2:160876/151‑1 (MQ=255)                                                                                                                              gggCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACa                  <  1:592270/151‑1 (MQ=255)                                                                                                                              gggCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACa                  <  1:907325/151‑1 (MQ=255)                                                                                                                                 cTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGc                >  2:412262/1‑150 (MQ=255)                                                                                                                                   g‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCt               <  1:1019207/149‑1 (MQ=255)                                                                                                                                     tttGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTAc             >  1:844875/1‑150 (MQ=255)                                                                                                                                      ttGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGAGTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCGATCGCGATAACAGCTa              <  1:66454/148‑1 (MQ=255)                                                                                                                                        gATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACaaa          <  2:234705/150‑1 (MQ=255)                                                                                                                                         aTGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACa        <  1:338816/151‑1 (MQ=255)                                                                                                                                          tGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAAc         >  1:1003275/1‑149 (MQ=255)                                                                                                                                                aaTATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCggg  <  1:192649/151‑1 (MQ=255)                                                                                                                                                aaTATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCggg  >  2:180407/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                         CGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTG‑TTTGATGC‑GGAATATCT‑TGTTGCCAGCAAA‑GCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGG  >  NC_000913/4237909‑4238199

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