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

CATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGG‑TAATAAACAGGATGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAA  >  NC_000913/3001015‑3001309                                                                                                                                                             |                                                                                                                                                  cATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGG‑TAATAAACAGGATGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCaa                                                                                                                                                   <  1:566180/150‑1 (MQ=255)  aTGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGG‑TAATAAACAGGATGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAAcc                                                                                                                                                 >  1:1110267/1‑151 (MQ=255)     cGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGG‑TAATAAACAGGATGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCAGTG‑AAGAGTGTTTCCTCGCAACCCg                                                                                                                                               >  1:737484/1‑150 (MQ=255)     cGTACGAGTCATCAAACAATTTGTCGGTGGCGGTTTTGG‑TAATAAACAGGATGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGt                                                                                                                                              <  2:1140244/151‑1 (MQ=255)        aCGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGG‑TAATAAACAGGATGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTAc                                                                                                                                            <  1:975545/150‑1 (MQ=255)                     ccATTTGTCGGTGGCGGTTTTGG‑TAATAAACAGGATGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCtttt                                                                                                                               >  2:1253584/1‑150 (MQ=255)                              ggTGGCGGTTTTGG‑TAAAAAACAGGGTGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCATCAAGGTTGGCGGGATTCCGGTGAAA‑GGTT‑TCC‑TTA‑GCCGTG‑AAGAGCGTTTCCTCGCAACCCGTACCCCCCACGCTTTTTCCATTg                                                                                                                        >  1:612112/1‑148 (MQ=255)                                tGGCGGTTTTGG‑TAATAAACAGGATGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGGGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACggg                                                                                                                   >  1:526798/1‑151 (MQ=255)                                  ggggTTTTGGAAAAAAAACAGGATGT‑CTGGAAAGGGCCAAGGGCGGCATTCCGGACCAGCAAGCTGGGCGGCATTCGGGTAAAT‑GTTC‑CCC‑TTA‑GCGGTG‑AAGAGTTTTTCCCCGCACCCCTCCCCCGCCACGCTTTTACCATTGACGGGc                                                                                                                  <  1:888863/149‑1 (MQ=255)                                           g‑TAATAAACAGGATGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGt                                                                                                        >  2:483572/1‑151 (MQ=255)                                              aataaACAGGATGTACGGG‑AAGAGCCAATGGCGGCATCCCTGACCAGCAGGCTGGGCGGCATTCCGGTAAAA‑TTTC‑CCC‑TTA‑GCCGTA‑AAGAGTGTTTCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGc                                                                                                          <  2:781664/147‑1 (MQ=255)                                               aaaaaCGGGATGTAGTGG‑AAGAGCAATGGGGGGCTTTCCGGACCAGAAGGTTGGGGGGCATTCGGGGAAAA‑TTTT‑CCC‑TTA‑GCCGGA‑AAGAGTGTTCCCTCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGaa                                                                                                     <  2:321161/149‑1 (MQ=255)                                                  aaCAGGATGTACGGA‑AAGGGCAAAGGGCGGCTTCCCTGACAAGCAAGCTTGGCGGAATCCGGGGGAAA‑TTTC‑CCC‑TTA‑GCCGTG‑AAGAGTGTTCCCCCGCACCCCTTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcg                                                                                                  <  1:322489/151‑1 (MQ=255)                                                   aCAGGATGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgc                                                                                                 >  2:124175/1‑151 (MQ=255)                                                      ggATGTCCGGG‑AAGAGCAAAGGGCGGCATCCCTGACCAGAAAGTTGGGCGGCA‑TCCGGGAAAA‑TTTTCCCC‑TTA‑GCCGTG‑AAGAGTTTTCCCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgcg                                                                                                <  1:1151626/149‑1 (MQ=255)                                                         tGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACgg                                                                                            >  2:713458/1‑150 (MQ=255)                                                         tGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACgg                                                                                            >  1:836162/1‑150 (MQ=255)                                                         tGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACgg                                                                                            >  1:1022846/1‑150 (MQ=255)                                                         tGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACgg                                                                                            >  1:1071140/1‑150 (MQ=255)                                                         tGTACGGA‑AAGAGCAAAGGGCGGCTTTCCTGACCAGCAAGTTGGGCGGCATCCGGGTAAAA‑GTTC‑CCC‑TTA‑GCCGGA‑AAGAGTTTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACg                                                                                             <  2:387322/149‑1 (MQ=255)                                                          gTACGGG‑AGGAGCAAGTGGCGCCATCCCTGCCCAGCAAGCTTGGCGGCATTCCGGT‑AAA‑TTTT‑CCC‑CTAGGCCGTG‑AAGTGTTTTCCCCCACACCCCTTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGAc                                                                                              <  1:1130545/147‑1 (MQ=255)                                                            aCTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAAc                                                                                         >  2:948408/1‑150 (MQ=255)                                                              tGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAtt                                                                                      >  2:1120399/1‑151 (MQ=255)                                                              tGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAtt                                                                                      >  2:231043/1‑151 (MQ=255)                                                               gg‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAtt                                                                                      >  1:319550/1‑150 (MQ=255)                                                               gg‑AAGAGCAAAGGGCGCCATTCCTAACCACCAAGCTGGGCGGC‑TTCCGGGAAAA‑ATTT‑CCC‑CTG‑GCCGGGAAAGAGTTTTCCCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACa                                                                                        <  1:158065/148‑1 (MQ=255)                                                                gaaaggcccAAGGGCGCCATCCCTAACCAGAAGGCTGGGCGGC‑TTGCGGGAAAA‑ATTC‑CCCTTTG‑GCGGTG‑AGGAGTTTTCCCCCCCACCCCGGCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAAc                                                                                         <  1:648678/140‑1 (MQ=255)                                                                gaaaGAGCCAGGGGCGGCTTCCCTGACCAGCAGGTTGGGCGGCATCCGGGTAAAT‑TTTC‑CCC‑TTA‑GCGGTA‑AAGGTTTTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTg                                                                                     <  2:1108605/149‑1 (MQ=255)                                                                  aaGGGCAAAGGGGGGCATCCCTGACCAGCAAGTTGGGCGGCTTCCGGGTGAAA‑TTTC‑CCC‑TTG‑GCCGTG‑AAGGGTTTTCCCCCGCACCCGTTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgc                                                                                                 <  2:1143778/137‑1 (MQ=255)                                                                  aaGAGCAAAGGGCGGCTTTCCGACCCGGCACGTTTGGGGGCATTCCGGTAAAA‑GTTC‑CCC‑TTG‑GCCGTA‑AAGAGTTTTCCCCCGCACCCCCTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTg                                                                                     <  2:760070/149‑1 (MQ=255)                                                                         aaTGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTtat                                                                            >  2:1102137/1‑151 (MQ=255)                                                                                   ttCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATg                                                                  >  2:790480/1‑151 (MQ=255)                                                                                   ttCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATg                                                                  >  1:489287/1‑151 (MQ=255)                                                                                          ccAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTc                                                           >  1:135681/1‑151 (MQ=255)                                                                                          ccAGAAAGCTTGGCGGCATTCGGGTGAAA‑TTTC‑CCC‑TTG‑GCCGTA‑AAGATTTTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCt                                                              <  2:683039/148‑1 (MQ=255)                                                                                           cAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGAGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCt                                                          >  1:186479/1‑151 (MQ=255)                                                                                              catgttgggcggcTTCCGGGTAAAA‑GTTC‑CCC‑TTG‑GCGGGG‑AAGGGTTTTCCCTCGCACCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTaa                                                        <  1:786620/143‑1 (MQ=255)                                                                                               aggCTGGGCGGCTTCCCGGT‑AAA‑TTTT‑CCCTTTA‑GCTGTG‑AGGAGTGTTTCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCt                                                          <  1:733332/145‑1 (MQ=255)                                                                                               agcgtgggcggcATCCGGGAGAAG‑TTTC‑CCC‑TGG‑GCCGTG‑AGGGGTTTTCCCCCCCACCCGGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTa                                                         <  1:1178914/142‑1 (MQ=255)                                                                                                ggCTTGGGGGCTTCCGGGTAAAA‑TTTT‑CCC‑TTG‑GCAGAG‑AAGTTTTTTCCCCCCCACCCCTCGCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAacac                                                     <  2:710681/150‑1 (MQ=255)                                                                                                   ttGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑CAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAATGGTTATAGTCTGGATGTTCTGTCTAacac                                                     >  2:667549/1‑148 (MQ=255)                                                                                                     ggcggcTTTCCGGTAAAA‑TTTT‑CCC‑TTG‑GCCGTG‑AAGATTTTTCCCTCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGgcg                                                <  2:873794/151‑1 (MQ=255)                                                                                                     ggcggcATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGgcg                                                >  1:756073/1‑151 (MQ=255)                                                                                                     ggcggc‑TTCCGGGAAAATGTTT‑CCT‑TTC‑GCGGAA‑AAGAGTTTTCCCCCCCCCCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACAcc                                                    <  2:1167529/147‑1 (MQ=255)                                                                                                     ggcgacttgctggTGAAA‑‑TTT‑CCT‑TGAGGCCGTG‑AGGGGTTTTTCCCGCCAACCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCgg                                                  <  2:612120/138‑1 (MQ=255)                                                                                                         ccacccgggTGAAG‑GTTT‑CCT‑TTC‑GCCGGG‑AGGGGTGTCCCCCCGCAGCCGGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCtt                                             <  1:580193/143‑1 (MQ=255)                                                                                                          cATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTAt                                           >  1:105840/1‑151 (MQ=255)                                                                                                                ggTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtct                                     >  2:46353/1‑151 (MQ=255)                                                                                                                 gTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtctc                                    >  1:272883/1‑151 (MQ=255)                                                                                                                 gTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtctc                                    >  2:1121605/1‑151 (MQ=255)                                                                                                                    aaa‑TTTC‑CCT‑TTA‑GCCGTG‑AAGGTTTTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACg                                 <  2:57419/151‑1 (MQ=255)                                                                                                                    aaa‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACg                                 >  2:252211/1‑151 (MQ=255)                                                                                                                    aaa‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACg                                 >  2:752616/1‑151 (MQ=255)                                                                                                                          tt‑CCCTTTA‑GCGGTG‑AAGAGTTTTCCCCCCCACCCCTTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGc                              <  2:1042127/150‑1 (MQ=255)                                                                                                                               c‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTcc                         >  1:118573/1‑150 (MQ=255)                                                                                                                                     agccgtaagggttgtTTCCCCGCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATc                      <  2:609569/138‑1 (MQ=255)                                                                                                                                      ccGTG‑AAGGGTGTTCCCTCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATc                      <  2:1233377/148‑1 (MQ=255)                                                                                                                                         tG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCt                >  2:54763/1‑151 (MQ=255)                                                                                                                                         ggaagggggctttCCACGCACCCCGTAGCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCtt                  <  1:469369/140‑1 (MQ=255)                                                                                                                                              gagTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgg            >  1:573876/1‑151 (MQ=255)                                                                                                                                               agTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggg           >  2:311655/1‑151 (MQ=255)                                                                                                                                                 tgtTTCCTCCCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgg            <  1:680585/148‑1 (MQ=255)                                                                                                                                                   tttCCCCGCACCCCGCACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggggg         <  1:1028062/149‑1 (MQ=255)                                                                                                                                                   ttcccccGCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgggggg        <  2:687745/144‑1 (MQ=255)                                                                                                                                                    ttCCCCGCACCCCGTAACCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgg            <  2:513547/145‑1 (MQ=255)                                                                                                                                                        ccGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATaa  <  2:1100547/150‑1 (MQ=255)                                                                                                                                                             |                                                                                                                                                  CATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGG‑TAATAAACAGGATGTACTGG‑AAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAA‑GTTT‑CCC‑TTA‑GCCGTG‑AAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAA  >  NC_000913/3001015‑3001309

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