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

CATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTC  >  NC_000913/3001015‑3001313                                                                                                                                                              |                                                                                                                                                      cATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCaa                                                                                                                                                       <  1:296850/150‑1 (MQ=255) cATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCaa                                                                                                                                                       >  4:534462/1‑150 (MQ=255)  aTGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAAcc                                                                                                                                                     <  4:349899/151‑1 (MQ=255)    gCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAAccca                                                                                                                                                   <  1:1090502/151‑2 (MQ=255)    gCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCg                                                                                                                                                   >  2:30804/1‑151 (MQ=255)     cGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGt                                                                                                                                                  <  2:1081937/151‑1 (MQ=255)     cGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGt                                                                                                                                                  <  2:1104915/151‑1 (MQ=255)      gTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTa                                                                                                                                                 <  2:472185/151‑1 (MQ=255)       tACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTAc                                                                                                                                                >  2:779715/1‑151 (MQ=255)         cGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTAc                                                                                                                                                <  1:309589/149‑1 (MQ=255)          gAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTAccc                                                                                                                                              >  2:199892/1‑150 (MQ=255)            gTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGcc                                                                                                                                           >  3:215726/1‑151 (MQ=255)               atcaAACAATTTGTCGGTGCCGGTTTTGGTAATAAACAGGATGTTCTGGAAGGGCCAAGGGCGGAGTT‑CCTGACCAGCATGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCCCGCACCCCGTACCCGcc                                                                                                                                           <  1:757589/148‑1 (MQ=255)                tcaAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGc                                                                                                                                       >  3:131838/1‑151 (MQ=255)                  aaaCCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTTCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTAACCGCCACGCtt                                                                                                                                     >  2:889309/1‑151 (MQ=255)                  aaaCCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCtt                                                                                                                                     >  2:1070762/1‑151 (MQ=255)                  aaaCCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGAACAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCtt                                                                                                                                     >  4:134344/1‑151 (MQ=255)                       aTTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTAcc                                                                                                                                >  2:477198/1‑151 (MQ=255)                           gTCGGTGGGGGTTGTGTAAAAAAACAGGGTGTAGTGAAAGAGCCAAGGGCGCCTTC‑CGTGACCAGTAAGCTGGGCGGC‑TTCCGGG‑TAAAGTTTCCCC‑TTAGCG‑G‑TAAAGAG‑TTTTT‑CCTCGCACCCCGTACCCGCCACGCTTTTAc                                                                                                                                 <  1:1127101/146‑1 (MQ=255)                              ggTGGCGGTTTTGGGAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGGTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACg                                                                                                                         >  4:27881/1‑151 (MQ=255)                                tGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGG‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACggg                                                                                                                       >  3:552958/1‑151 (MQ=255)                                 ggtggTTTTGAAAAAAAACGGTGTGTACGGGAAGGCCAAAGGGGGGCTTT‑CCTGCCCACCAGGGTTGGC‑GCTTTCCGGTAGAACTTTC‑CCC‑TTCGCC‑G‑TAAGGGG‑TGTTC‑CCCCGCACCCCGTTCCCGCCACGCTTTTACCATTGACGGGc                                                                                                                      <  4:531428/148‑1 (MQ=255)                                 ggCTGTTTTGGAAAAAAACAGGTTGTAGTGGAAGGGCAAGTGGCGTCTTT‑CCTGACCAGCATGTTGGGCGGCTTTCGGGT‑AAAATTTC‑CCC‑TTACCG‑G‑TGAAGGG‑TTTTC‑CCCCGCACCCCTTACCCGCCACGCTTTTACCATTGACGGGc                                                                                                                      <  1:484781/151‑1 (MQ=255)                                    ggTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCaa                                                                                                                    >  1:469538/1‑150 (MQ=255)                                     gTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAAt                                                                                                                  >  3:197275/1‑151 (MQ=255)                                      ttttGAAAAAAAACAGGATGTCGGGGAGGCGCTATGGGCGTCTTT‑CCTCACCAGCAGTGTGGGGGGCATTCGGGT‑‑AAAGTTT‑CCCTTTAGCC‑G‑TGAAGGG‑TTTTC‑CCCCGCACCCCTTACCCGCCACGCTTTTACCATTGACGGGCa                                                                                                                     <  1:238502/147‑1 (MQ=255)                                         tGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGc                                                                                                              >  2:197249/1‑151 (MQ=255)                                          gttaataaACAGGATTACCGGGAAGAGCCAAGGGCGCCATT‑CCTGACCAGCAAGTTGGGCGGCTTTCCGGT‑AAAATTTC‑CCC‑TTAGCC‑G‑TGAAGAG‑TTTTC‑CCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATggg                                                                                                               <  3:211984/147‑1 (MQ=255)                                          gaaaaaaaaCAGTATGCACTGAAAGAGCAAGTGGCGATATC‑CCTCCCCAGAAAGCTGGGGGGCTTCCGGGT‑AAAATTTT‑CCT‑TTGGCA‑G‑TAAAGAT‑TTTTC‑CCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGc                                                                                                              <  2:219589/144‑1 (MQ=255)                                             aataaACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTg                                                                                                           >  2:565453/1‑150 (MQ=255)                                             aaaaaaCGGGTTGTAGGGAAAGAGCAAGTGGCGGCATT‑CCTGACCAGCAAGCTGGGCGGCATTCGGGT‑AAAAGTTT‑CCC‑TTGCCC‑G‑TGAAGAT‑TTTTC‑CCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATggg                                                                                                               <  2:1468813/143‑1 (MQ=255)                                              ataaACAGTATTTACTGGAAGACCCAAGGCCGGCATT‑CCTCCCCACCAGGCTTGGCGCCTTTCCGTT‑AAAATTTT‑CCC‑TTCGCC‑T‑TGAATTT‑TTTTC‑CCTCGCACCCCTTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCg                                                                                                             <  2:692139/147‑1 (MQ=255)                                               aaaaCAGGATGTCGGGGAAGAGCAAAGGGCGCCATC‑CCGGACCAGCAAGCTGGGCGCCTTCCGGGT‑AAAATTTC‑CCC‑TTAGCC‑G‑TGAAGGT‑TGTTC‑CCGCGAACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAAc                                                                                                        <  1:604122/150‑1 (MQ=255)                                                aaaCAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAAcc                                                                                                       >  2:547073/1‑151 (MQ=255)                                                aaaCAGGATGTACGGAAGGAGCAAATGGCGGCATT‑CCTGACCAGAAGGCGGGGCGCCTTCCGGGT‑GAAATTTT‑CCC‑TTAGCC‑G‑GGAGGAG‑TTTTT‑CCTCGAACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAAcc                                                                                                       <  2:1074048/151‑1 (MQ=255)                                                    aGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGAAAGCTTGGCGGAATTCCGGG‑GAAAGTTT‑CCC‑TTAGCC‑G‑CGAAGCG‑TGTTT‑CCTCGCCACCCGTACACGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgc                                                                                                     >  1:1549199/1‑149 (MQ=255)                                                       aTGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACgg                                                                                                >  2:1466595/1‑151 (MQ=255)                                                            ataaaaGAG‑CAAAGGCCGCATTCCCTGACCAGCAGTCTGGGCGCCTTCCGGGT‑AAAATTTC‑CCC‑TTGGCC‑G‑GGAAGAG‑TGTTC‑CCCCGCACCCGTTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAt                                                                                           <  2:1087711/147‑1 (MQ=255)                                                             tGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAtt                                                                                          >  4:71035/1‑151 (MQ=255)                                                                aaggcccAAGGGCGGCTTT‑CCTCCCCAGCAGGTTGGGCGGCCTTCGGGT‑AAAATTTT‑CCC‑TTGGCC‑G‑TAAAGGT‑TTTTT‑CCCCGCACCCCTCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACa                                                                                            <  2:493449/141‑1 (MQ=255)                                                                aaGAGCAAGTGGCGGCTTT‑CCTGACCAGCAAGCTGGGGGGCTTCCGGGT‑AAAATTTT‑CCC‑TTAGCC‑G‑GGAAGAG‑TTTTC‑CCTCGCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGaa                                                                                       <  2:448425/151‑1 (MQ=255)                                                                   agCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAgg                                                                                    >  2:128243/1‑151 (MQ=255)                                                                     ccAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGtt                                                                                  >  3:470930/1‑151 (MQ=255)                                                                     caaggggcggcTTT‑CCTGCCCAGAAAGTTGGGCGGC‑TTCCGGG‑TAAATTTT‑CCCTTTAGCG‑G‑GGAGGGT‑TTTTC‑CCTCGCACCCCTTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAgg                                                                                    <  4:552938/144‑1 (MQ=255)                                                                      cgttGGCGTCTTT‑CCTGCCAAGC‑AGGTTGGCGGCATTCCGGG‑TAAAATTT‑CCC‑TTGGGC‑GTTGAGGAG‑TTTTT‑CCTCCCAACCGTCACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGaaa                                                                                      <  3:534479/143‑1 (MQ=255)                                                                      ccatggggggtttcccTGCCCAGCAAGTTTGGCGGCTTCCCGGT‑GAAAGTTT‑CCC‑TGGGCC‑G‑TGAAGTG‑TTTTT‑CCCCGCACCCCCTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGtt                                                                                  <  1:1344003/137‑1 (MQ=255)                                                                       aggggcggcATC‑CCGGACCAGCAAGCTGGGCGGC‑TTCCGGG‑TAAAAGTT‑CCC‑TTGAGC‑G‑GGAAGGGTTTTTT‑CCTCCCACCCCCTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAg                                                                                     <  4:273722/143‑1 (MQ=255)                                                                       aagggcggcTTT‑CGTGCCCGGCAAGCTGGGCGGCATTCCGTA‑AAAATTTC‑CCC‑TTCGCG‑G‑TGAAGAT‑TTTTT‑CCTCCCACCCCCCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGtt                                                                                  <  4:411297/146‑1 (MQ=255)                                                                        aTGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTtat                                                                                >  2:595466/1‑150 (MQ=255)                                                                          ggcggcATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGt                                                                             >  3:218882/1‑151 (MQ=255)                                                                          ggcggcATC‑CCTGCCCAGAAAGCGTGGCGGCATCCGGGT‑AAAAGTTC‑CCC‑TTAGCG‑G‑GAAAGAG‑TGTTC‑CCCCCCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAg                                                                              <  2:595489/150‑1 (MQ=255)                                                                          ggcgccATT‑CCTGACCAGCAAGCTGGGCGGCATTCCGGT‑AAAAGTTT‑CCC‑TTGGCC‑G‑AGAAGAG‑TGTCT‑CCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGt                                                                             <  1:565370/151‑1 (MQ=255)                                                                               cATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTgg                                                                         >  2:1489603/1‑150 (MQ=255)                                                                                     ctgccaaGCAGGCTGGGCGGCTTTCGGGT‑AAATTTTC‑CCC‑TGAGCG‑G‑TGAGGTG‑TGTTC‑CCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGt                                                                     <  1:1005989/143‑1 (MQ=255)                                                                                      tGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCt                                                                  >  2:629846/1‑151 (MQ=255)                                                                                         caagaaaGCTTGGCGGCTTGCCGTA‑AAAATTTC‑CCC‑TCAGCC‑‑‑TGAAAGT‑TGTTTCCCCCACACCCTGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGt                                                                <  2:1488120/145‑1 (MQ=255)                                                                                           agaaggttgggcggcTTTCCGTT‑AAAATTTC‑CCT‑TTAGCG‑G‑GGAGGTG‑TGTTC‑CCCCGCACCCCTTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTg                                                                 <  1:864105/138‑1 (MQ=255)                                                                                           agaaagttgggcggcATTCCGGT‑AAATGTCC‑CCA‑TTGGCC‑T‑GAAAGAG‑TTTTT‑CCGCCCCCCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTg                                                                 <  1:943680/138‑1 (MQ=255)                                                                                             caggCTTGGCGG‑TTTCCGGG‑TGAAGTTT‑CCC‑TTTGGC‑GTTGAAGGT‑GGTTC‑CCCCCCACCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCt                                                              <  1:261727/145‑1 (MQ=255)                                                                                             caggCTTGG‑GGGATTCCCGG‑GAAATTTTCCCC‑TTGGCC‑G‑GGAAGGG‑TTTTT‑CCCCGCACCCCTTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTc                                                               <  3:10379/144‑1 (MQ=255)                                                                                              agttttGGCGACATTCGGGT‑GAAATTTC‑CCC‑GTGGCC‑G‑AGAAGAT‑TTTTC‑CCCCGCCCCACGCACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCt                                                              <  2:1698266/143‑1 (MQ=255)                                                                                              aaGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAac                                                           >  2:533938/1‑150 (MQ=255)                                                                                                gctggggggttttCCGGA‑AAAATTTC‑CCT‑TTAGCG‑G‑AGAAGAG‑TGTTC‑CCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAacac                                                         <  2:1547428/139‑1 (MQ=255)                                                                                                 cTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCg                                                       >  1:1609093/1‑151 (MQ=255)                                                                                                  gtGGCCGCTTTCGGGT‑AAATG‑TT‑CCC‑TTGGCCGG‑TGAAGGT‑TTTTC‑CCTCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAac                                                           <  1:25290/145‑1 (MQ=255)                                                                                                    ggcggcATCCGGGT‑GAAATTTC‑CCC‑TTACCC‑G‑GAAAGGG‑TTTTC‑CCCCCCACCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCg                                                       <  1:1670280/148‑1 (MQ=255)                                                                                                    ggcggcACCCCGGT‑GAATGTTC‑CCC‑TTAGCC‑G‑GGAAGAG‑TTTTC‑CCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCg                                                       <  1:459002/148‑1 (MQ=255)                                                                                                         cATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTAt                                               >  2:1386045/1‑151 (MQ=255)                                                                                                            tgcggttaaatgTTC‑CCT‑TTAGCC‑G‑TGAAGGT‑TTTTC‑CCGCGCCACCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATg                                              <  1:258392/139‑1 (MQ=255)                                                                                                            tCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCt                                                                  >  4:517417/1‑129 (MQ=255)                                                                                                                 t‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtct                                         >  1:776507/1‑149 (MQ=255)                                                                                                                    aaaGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACg                                     >  1:202710/1‑151 (MQ=255)                                                                                                                     aaGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACgg                                    >  4:2602/1‑151 (MQ=255)                                                                                                                      aGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACggg                                   >  3:596663/1‑151 (MQ=255)                                                                                                                      aGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACggg                                   >  2:518158/1‑151 (MQ=255)                                                                                                                         ttttccttggccc‑G‑GGAAGTG‑TGTTC‑CCTCGCACCCTGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACgg                                    <  4:474091/138‑1 (MQ=255)                                                                                                                          t‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACt                               >  2:1736482/1‑151 (MQ=255)                                                                                                                            ccc‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTc                              >  1:560904/1‑151 (MQ=255)                                                                                                                            ccc‑TTAGAC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGGTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCAc                                >  2:1104340/1‑149 (MQ=255)                                                                                                                             cc‑TTAGCG‑G‑TGAAGTG‑TTTTC‑CCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTcc                             <  2:630702/151‑1 (MQ=255)                                                                                                                             cc‑TTAGCC‑G‑TGAAGAG‑GGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTcc                             >  4:545275/1‑151 (MQ=255)                                                                                                                                ttAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCAt                           >  1:301596/1‑151 (MQ=255)                                                                                                                                 tAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATc                          >  1:1178734/1‑151 (MQ=255)                                                                                                                                    cgatctaagggg‑tgtTC‑CCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCAt                           <  3:62870/138‑1 (MQ=255)                                                                                                                                     c‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCtt                      >  1:1615224/1‑151 (MQ=255)                                                                                                                                     c‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCtt                      >  2:386291/1‑151 (MQ=255)                                                                                                                                       gggaaaGTG‑TGTTC‑CCTCGCCCCCCCCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCtt                      <  1:1353299/147‑1 (MQ=255)                                                                                                                                          gAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCAt                           >  4:62868/1‑143 (MQ=255)                                                                                                                                           aggag‑TATTT‑CCACGTACCCCACCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCt                    <  1:128961/147‑1 (MQ=255)                                                                                                                                           aaGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGc                  >  1:100941/1‑151 (MQ=255)                                                                                                                                           aaGAG‑TGTTC‑CCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCt                    <  4:215708/149‑1 (MQ=255)                                                                                                                                                   ttt‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggggggg           >  1:580725/1‑151 (MQ=255)                                                                                                                                                   ttt‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggggggg           >  2:1092470/1‑151 (MQ=255)                                                                                                                                                    ttcccccGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgggggg            <  2:1031950/144‑1 (MQ=255)                                                                                                                                                          cGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATaa      >  1:437793/1‑150 (MQ=255)                                                                                                                                                            cAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGt   >  2:1620479/1‑151 (MQ=255)                                                                                                                                                             aaCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGt   >  2:1526998/1‑150 (MQ=255)                                                                                                                                                             aaCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTc  >  1:1468657/1‑151 (MQ=255)                                                                                                                                                             aaCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTc  >  2:1475848/1‑151 (MQ=255)                                                                                                                                                              |                                                                                                                                                      CATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATT‑CCTGACCAGCAAGCTTGGCGGCATTCCGGT‑GAAAGTTT‑CCC‑TTAGCC‑G‑TGAAGAG‑TGTTT‑CCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTC  >  NC_000913/3001015‑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: ‑