Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,001,164	0	A	C	43.2%	124.4 / 18.1	74	T274P (ACC→CCC)	xdhA	xanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base A (36/6);  new base C (0/32);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 3.15e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTC  >  NC_000913/3001020‑3001313                                                                                                                                                         |                                                                                                                                                         gTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑Ta                                                                                                                                                 <  1:1366581/151‑1 (MQ=255)    cGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TAccc                                                                                                                                              <  2:304269/151‑1 (MQ=255)    cGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGGGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TAccc                                                                                                                                              >  2:1452153/1‑151 (MQ=255)     gAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TAccca                                                                                                                                             <  2:536064/151‑2 (MQ=255)        tcatcaAACCATTTGTCGGTGGCGGTTTTGGTAATAAACCGGATGAAATGGAAGAACCCA‑TGGCGGCAATCCAGGACCA‑CAAGCATGGCGGCTTTTCCGGGTAA‑GTTTTCC‑TT‑AA‑TCGTGGAGTGTGATT‑CCTCCCAAA‑CCC‑G‑AAcc                                                                                                                                               >  1:1293316/1‑147 (MQ=255)          atcaAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCAc                                                                                                                                         >  1:706977/1‑150 (MQ=255)          atcaAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCAc                                                                                                                                         >  1:488395/1‑150 (MQ=255)                ccATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCtttt                                                                                                                                   >  1:361341/1‑150 (MQ=255)                 cATTTGTCGGTGGCGGGTTTGGTAATAAACAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGGGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTAc                                                                                                                                 >  1:1295881/1‑151 (MQ=255)                      gttggTGGGGTTTTTGAAAAAAAACAGGTTGTACGGGAAG‑GCCAAGTGGCGGCATTCG‑GAACCGGCAGGTTTGGCGAAATTCCGGGAAAA‑TTTCCCC‑TT‑TC‑CGGTAAAGGTTGTTT‑CCGCGC‑AC‑CCC‑G‑TACCCGCCACGCTTTTACCAtt                                                                                                                             <  1:447484/147‑1 (MQ=255)                         gggggggtttttGGAAAAAAACAGGATGTAGGGAAAGGGCAAG‑TGGCGCTTTTCC‑TGCCCAACAGGTTGGGGGCCTTTCGGGTGAAA‑TTTCCCT‑TT‑AG‑CCGGAAAGGTTTTTT‑CCCCCC‑AC‑CCC‑G‑TCCCCGCCACGCTTTTACCATTGACg                                                                                                                         <  2:940034/143‑1 (MQ=255)                         gggggCGTTTTTGGAAAAAAACAGGATGTACGGAAAGAGCAAA‑GGGCGGCATTCC‑TGACCAGCAAGCTGGGCGGCTTTCGGGAAAAA‑TTTTCCT‑TT‑GG‑CCGTAAAGAGTTTTC‑CCCCGC‑AC‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACg                                                                                                                         <  1:1397778/148‑1 (MQ=255)                             gCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGc                                                                                                                      >  1:591939/1‑150 (MQ=255)                                  tttGGTAATAAACAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATgg                                                                                                                >  2:1267846/1‑151 (MQ=255)                                  tttGGAAAAAAACAGGATGTCGGGAAAGCCCAAA‑GGGGGGCATCCC‑GAACCAGAAAGCTGGGCGCCTTCCGGGAAAAA‑GTTCCCC‑TT‑GG‑CGGTAAAGGGTTTTC‑CCCCGC‑AC‑CCC‑G‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATg                                                                                                                 <  1:1141811/150‑1 (MQ=255)                                    tGGTAATAAACAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGc                                                                                                              >  1:277791/1‑151 (MQ=255)                                        aaaaaaCAGGGTGTAGGGGAAGAGCAAA‑GGTCGCCCTTCC‑GGCCCAGAAAGTGGGCCTGCA‑TCCGGGTAAA‑GTTTTCC‑TTCGG‑ACGAAAAGAGGGTTC‑CCTCGC‑AC‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGc                                                                                                              <  2:184291/144‑1 (MQ=255)                                        aaaaaaCAGGATGTAGTGGAGGAGCAAG‑GGGCGGCTTCCC‑TGACCAGCAAGTTGGGCGCCATTCCGGTAAAA‑TTTCCCC‑TT‑CG‑CCGTGAAGAGTTTTC‑CCCCCC‑AC‑CCC‑G‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGt                                                                                                            <  1:715376/146‑1 (MQ=255)                                           aaaCGGGATGTGGTGAAAGGGCAAA‑GGGGGGCTTTCC‑AGCCCAGAAGGTTGGGGGGCCCTCGGGTGAAT‑TTTTCCC‑TT‑AG‑CGGTGAAGTTTGTTC‑CCCCGC‑AC‑CCT‑G‑TTCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAAcc                                                                                                       <  1:715810/151‑1 (MQ=255)                                             aCAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgc                                                                                                     >  2:632197/1‑151 (MQ=255)                                              cAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgcg                                                                                                    >  2:918783/1‑151 (MQ=255)                                              cAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgcg                                                                                                    >  1:1046347/1‑151 (MQ=255)                                              cAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgcg                                                                                                    >  1:1441357/1‑151 (MQ=255)                                                   tGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACgg                                                                                                >  2:334399/1‑150 (MQ=255)                                                     tACGGGAAGAGCCAG‑GGGCGGCATCCC‑GGACCAGAAGGCTTGGCGGCTTTCGGGTAAAA‑TTTCCCT‑TT‑GG‑CCGGGAAGGGTTTTC‑CCCCGC‑AC‑CCC‑T‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGa                                                                                               <  2:654187/149‑1 (MQ=255)                                                     tACGGGAAGAGCCAA‑GGGCGGCATTCC‑TGACCAGCAAGTTTGGCGGCATTCGGGTAAAA‑TTTCCCT‑TT‑GG‑CCGTGAAGATTTTTT‑CCTCGC‑AC‑CCC‑G‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACg                                                                                                 <  1:702324/147‑1 (MQ=255)                                                       cTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGCAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAt                                                                                           >  1:375473/1‑151 (MQ=255)                                                       cTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAt                                                                                           >  2:418386/1‑151 (MQ=255)                                                           aaGAGCCAG‑GGGCGGCTTTCC‑GGACAAGCAAG‑GTGGGGGGAATCCGGGTAAA‑ATTTTCC‑TT‑AGCCCGTGAGGAGTTTTC‑CCCCGA‑AC‑CCC‑G‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAtt                                                                                          <  1:1128551/148‑1 (MQ=255)                                                           aaGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGa                                                                                        >  2:170758/1‑150 (MQ=255)                                                           aaGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGa                                                                                        >  2:553515/1‑150 (MQ=255)                                                                ccAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGtt                                                                                  >  1:94811/1‑151 (MQ=255)                                                                 caaagggcggcATTCC‑TGACCAGCAAGCTGGGCGGCATTCGGGTAAAA‑GTTTCCC‑TA‑AG‑CCGTGAAGAGTTTTC‑CCCCGC‑AC‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGaa                                                                                       <  1:70946/141‑1 (MQ=255)                                                                 cAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGtt                                                                                  >  2:1147989/1‑150 (MQ=255)                                                                     tGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTtat                                                                                >  1:1294990/1‑149 (MQ=255)                                                                      ggcggcTTCCC‑TCACCAGAACGCTGGGCGGCTTCCGGGTAAAA‑GTTCCCC‑TT‑AG‑CGGGGAAGAGT‑TTT‑CCCCCC‑CA‑CCC‑GTTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGt                                                                             <  2:828085/151‑1 (MQ=255)                                                                      ggcggcATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAg                                                                              >  2:497645/1‑150 (MQ=255)                                                                         ggcATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTg                                                                          >  1:1358883/1‑151 (MQ=255)                                                                         gcctttCC‑TGACCAGCAAGTTGGGCGGCATCCCGGTGAAA‑GTTCCCC‑TT‑AG‑CCGGAAAGAGTTTTT‑CCCCGC‑AC‑CCC‑G‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGaaa                                                                                      <  2:292957/135‑1 (MQ=255)                                                                             ttCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATg                                                                      >  1:305550/1‑151 (MQ=255)                                                                             tccc‑TGCCCAGCAAGCTGGGCGCCATCCCGGTAAAA‑TTTTCCC‑TC‑GG‑CGGTAAAGAGTTTTT‑CCCCCC‑AC‑CCC‑T‑GCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTgg                                                                         <  1:1087696/146‑1 (MQ=255)                                                                               cc‑TGCCCAGCAAGCTGGGGGACATCCCGGT‑AAA‑ATTTCCC‑TT‑TGCCCATGAGGGGTGTTC‑CCCCGC‑GC‑CCC‑G‑ATCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTg                                                                          <  1:883123/145‑1 (MQ=255)                                                                               cc‑TGACCAGAAAGTTGGGCGGCATTCGGGGGAAA‑TTTCCCC‑TT‑GG‑CGGTGAAGTGTTTTC‑CCTCGC‑AC‑CCC‑G‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGt                                                                     <  2:366315/150‑1 (MQ=255)                                                                                c‑TGCCCAGAAAGTTGGGCGGCTTCCCGGTGAAT‑TTTCCCC‑TG‑GG‑CCGGGAAGGGTGTTC‑CCCCGC‑AC‑CCC‑G‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGtt                                                                    <  2:286774/150‑1 (MQ=255)                                                                                   gACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTg                                                                 >  1:495316/1‑151 (MQ=255)                                                                                      ccagaaGGTTGGCGGGATTCCCGTTAAATGTTCCCC‑TA‑AG‑CCGGAAAGTGTTTTC‑CCCCCC‑AC‑CCC‑G‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCt                                                                  <  2:229203/144‑1 (MQ=255)                                                                                       aGCAGGTTTGGCGGCTTCCGGGTAAAA‑TTTCCCC‑TT‑AG‑CCGTGAAGAGTTTTT‑CCTCGC‑AC‑CCC‑G‑AACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGt                                                                <  1:809223/148‑1 (MQ=255)                                                                                       aGCAAGTTGGGCGGC‑TTCCGGGTAAA‑ATTTCCCTTT‑AG‑CGGTAAAGGTTTTTC‑CCTCCC‑AC‑CCC‑T‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCt                                                              <  2:56324/150‑1 (MQ=255)                                                                                          atgttgggcggcTTTCGGGT‑AAA‑TTTTCCCTTT‑AG‑CCGTGAAGATTGTTC‑CCCCGC‑AC‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAaca                                                          <  1:208666/145‑1 (MQ=255)                                                                                          aaGCTGGGCGGCTGTCCGGTAAAA‑TTTCCCC‑TT‑AC‑CCGGGAGGTGTTTTT‑CCCCGC‑AC‑CCC‑G‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAac                                                           <  1:845642/150‑1 (MQ=255)                                                                                              tgggggtcatccgggTGAAT‑TTTTCCT‑TT‑GG‑CGGTGAAGTGTTTTC‑CCTCAC‑AC‑CCT‑G‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACAcc                                                        <  1:296892/136‑1 (MQ=255)                                                                                               tGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGgc                                                     >  1:1071887/1‑151 (MQ=255)                                                                                                gggcgctttCCGGGAAAA‑TTTTCCC‑TT‑AG‑GCGGGAAAGGTTTTT‑CCCCCC‑CC‑CCC‑C‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCgg                                                      <  2:903525/142‑1 (MQ=255)                                                                                                  cggcATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCt                                                  >  2:847130/1‑151 (MQ=255)                                                                                                     cATGCCGGTGAAC‑TTTCCCA‑TT‑AG‑CCG‑GAAGGTTTTTTCCCCCGC‑AC‑CCG‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTa                                                <  1:163663/150‑1 (MQ=255)                                                                                                     cATACCGGAGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGCCTGGATGTTCTGTCTAACACCGGCGCTTAt                                               >  2:266228/1‑151 (MQ=255)                                                                                                         cgggTAAAA‑TTTCCCC‑TT‑AG‑CGGTGAAG‑GTTTTT‑CCCCCC‑AACCCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCa                                            <  1:1245561/148‑1 (MQ=255)                                                                                                            ggtaaattttcccc‑TT‑AG‑CCGTAAGGTGT‑TTT‑CCCCCC‑AA‑CCC‑GCTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtct                                         <  1:556689/140‑1 (MQ=255)                                                                                                               aaa‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACg                                     >  2:143083/1‑151 (MQ=255)                                                                                                               aaa‑GTTCCCC‑TT‑GG‑CCGTGAAGGGTGTTC‑CCCCGC‑AC‑CCC‑G‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtctc                                        <  1:985365/148‑1 (MQ=255)                                                                                                                 aattttCCC‑TT‑AG‑CCGTGAAGTGTTTTT‑CCTCGC‑AC‑CCC‑T‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCa                                       <  2:706897/145‑1 (MQ=255)                                                                                                                        cc‑TT‑AG‑CCGTAAAGAGTTTTC‑CCCCGC‑AC‑CCC‑G‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTcc                             <  1:444631/151‑1 (MQ=255)                                                                                                                           tt‑AG‑CCGTGAAGTTTGTTC‑CCACCC‑AC‑CCC‑G‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTcc                             <  2:1264972/149‑1 (MQ=255)                                                                                                                            tcaG‑CCGTAAAGAGTTTTT‑CC‑CCC‑AC‑CCCGG‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCa                            <  2:361292/148‑1 (MQ=255)                                                                                                                              tggcggggaagtttgtTC‑CCCCGC‑AC‑CCC‑G‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCAt                           <  1:356532/137‑1 (MQ=255)                                                                                                                                 ccGTGAAGGTTTTTC‑CCCCCC‑AC‑CCC‑G‑TCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCg                         <  1:695792/149‑1 (MQ=255)                                                                                                                                  cGTGAAGGGTGTTC‑CCCCGC‑AC‑CCC‑G‑GCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCt                       <  2:487616/150‑1 (MQ=255)                                                                                                                                  cGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCtt                      >  1:868495/1‑151 (MQ=255)                                                                                                                                    tGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCt                    >  1:124837/1‑151 (MQ=255)                                                                                                                                      aaGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGc                  >  1:76204/1‑151 (MQ=255)                                                                                                                                             ttt‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggggggg           >  2:510847/1‑151 (MQ=255)                                                                                                                                                  cTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAAt        >  2:320788/1‑150 (MQ=255)                                                                                                                                                     gC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATaa      >  1:74164/1‑149 (MQ=255)                                                                                                                                                     gC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATaa      >  2:27651/1‑149 (MQ=255)                                                                                                                                                        aa‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTc  >  1:136540/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                         GTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAA‑TGGCGGCATTCC‑TGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTTCCC‑TT‑AG‑CCGTGAAGAGTGTTT‑CCTCGC‑AA‑CCC‑G‑TACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTC  >  NC_000913/3001020‑3001313

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