Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,278,682	C→T	A68A (GCC→GCT)	ghxP →	guanine/hypoxanthine transporter GhxP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,278,682	0	C	T	96.4%	179.6 / NA	57	A68A (GCC→GCT)	ghxP	guanine/hypoxanthine transporter GhxP
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (28/25);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGGTTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGG‑TTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGG  >  NC_000913/4278445‑4278927                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      tGACGAAAATTAAACCCTAGGGGGATTTTTTCTTATGTCTGCGCCATAAGCGCGTACCGGCGGTTCCCTCGACGCCT‑GG‑TTTAAAATTTCACA‑TCG‑TGGAAGC‑TCTGTCCGTCAGGAAGGGGT‑TGGC‑GGGTTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGTT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTc                                                                                                                                                                                                                                            <  2:259379/249‑1 (MQ=255)     gAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGGTTAACAACGTTTCTTGCGATGGTCTACTC‑GTTCATCGTCGTTCCAGGTATGT‑GGGGT‑AAAGCAGGCTTCCCGTCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTTGACTCGGTTCTATCGt                                                                                                                                                                                                                                      <  1:88102/251‑1 (MQ=255)             aaCTCTCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑GGGGAAC‑ACTTTCCGTCAGGAAGGGGG‑TGCC‑GGGTTAACAACGTTTCTGGCGGTGGGCTACTC‑GGACATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCACGCCGGCGGCCGCTTTCGTTGCAACCTTTGTGG‑TTGCTGGGCTCGGTTCTATCGCGATggggc                                                                                                                                                                                                                              >  1:353759/1‑249 (MQ=255)               ctctCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGGCGCCT‑GGGTTTAAAATT‑CACACACG‑TG‑AGGC‑CCTGTCCGCCCGGAAGTGGT‑TGCC‑GGGGTTACACAGTTTCTGGCGATGGGCTACTC‑GGTCATCGTCGGTCCCGGTCTGT‑TGGGT‑AAAGCAAGCGTCCCCCCTGGGGCAGCTTTCGTTTCCACCTGTCTGG‑TTGTCTGACTCGGGTCtttccgcattggtctg                                                                                                                                                                                                                            >  1:353945/1‑236 (MQ=255)                   cAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTTACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCCGGAAGGGGT‑TGCCCGGGTTA‑CAACGGTTCTGGCGGTGGTCTCCTC‑GGTCATCGTCGTTCCCGGGATGT‑TGGGG‑AAAGCAGGCTTCCCCCCTGCGGGACCCTTCGTTGCCACCTTTCTGG‑GTGTTGGGCTCGGGTCTATCTTTATGGGGCTTTggg                                                                                                                                                                                                                        >  2:50005/1‑251 (MQ=255)                         aTGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTAACTCGAAGCCC‑GG‑TTTAAAATTTCACA‑ACC‑TGGAAGC‑ACTGTCCGGCAGGGAGTGGG‑TGCC‑GGGGTAACAACGTTTCTTGCGATGGTCTACTC‑GGTCATCGTCGTTTCAGGGATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGGCGCTTTCGGTGCAACCTGTCTGG‑TTGCTGGGCTCGGGTCTTTCGTGATGGGGCTGTGGGGTAATc                                                                                                                                                                                                                  >  1:146123/1‑251 (MQ=255)                                 ttATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTTGACCCCT‑AG‑ATTAAAATTTTACA‑ACGAGGGAGGC‑A‑TGTCCGTCAGGAAGCGGT‑TGCC‑GGGGTAAAAACGCTTCTTGGGCTGGTCTTAAC‑GGGCATAGTCGTTCCAGGGATGT‑TGGGG‑AAAGAAGGCTTCCCCCCTGCGGGCGCATTCGTTGCAAACCGGCTGG‑TTGCTGGAGTCGGGTCTATCGTGGTGGGGCTGTGGGGTAATCCGCCGTTg                                                                                                                                                                                                          >  1:264854/1‑251 (MQ=255)                                   aTGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGAAGCCCCGG‑TTTAAA‑TTTCACA‑ACG‑AGGAAGC‑ACTGTACGTCAGGCAGTGGT‑TGCC‑CGGTTTACAACGTTTCTTGCGGTGGTATACTC‑GGCCATCGTCGTTCCACGAATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCCACACGTCTGGGTTGT‑GGGCTCGGTTCTATCGTGATGGGTCTGTGGGCTTATCTTCCGTTGGc                                                                                                                                                                                                        >  1:385898/1‑251 (MQ=255)                                             cATCAGCGCGTACCGGCGGTTCACCCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑CCTGTCCGTCAGGAAGTGTT‑TGCC‑GGGTTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTgc                                                                                                                                                                                              <  2:125884/251‑1 (MQ=255)                                                     ggTACCGGGGGTTCACTCGACGCCT‑GG‑TTAAAAATTTCACA‑CTG‑TGGAACC‑ACTGTCCGTCAGGAAGGGGT‑TGCC‑GGGTTAAAAACTGTCCTGGCGTTGTTCTTCTC‑GGGCATCGTCGTTCCAGGTCGGT‑TGGGT‑AAAGCGGTCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCGGTCTGG‑TTGCTGGACTCGGTTCTGTCGTGATGGGTCTGTGTGCTAATCTGCCGTTGGCGATTGGTTGCGCCAtctc                                                                                                                                                                                      <  2:353945/250‑1 (MQ=255)                                                      gtacctgtgtttCACTCGACGCGT‑GT‑TTTAAAATTTCAAA‑GCG‑TGGAAGC‑ACTGTCCGTCAGGAGGTGGT‑TGCC‑GGGTTAACAACGTTTCTGGCGAGGGTCTACCC‑GGTAATCGTCGTTCCGGGTATGT‑TGGGT‑AAGGCGGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTcc                                                                                                                                                                                     <  2:189486/241‑1 (MQ=255)                                                         ccGGCGGTTACCTCGACGCCT‑GG‑TTTAAACTTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGGTTACCAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTg                                                                                                                                                                                  >  2:337485/1‑251 (MQ=255)                                                            gcggTTCACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGGTTAACAACGTTTCTGGCGAGGGTCTACTG‑GGTCATCTTCGTTCCAGGTATGT‑TGGTT‑AAAGCAGTCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGGGATGGGTTGCGCCATCTCCCTGACa                                                                                                                                                                               <  1:126676/251‑1 (MQ=255)                                                             cggTTCACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCCGGAAGTGGT‑TGCC‑GGGCTAACAACGTTTCTGGCGATGGTCCACTC‑GGTCATCGTAGTTCCCGGCATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGCCTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGAc                                                                                                                                                                                >  1:204576/1‑249 (MQ=255)                                                              ggttacctcggacccT‑GG‑TTTAAAATTTACCA‑ACG‑TGGAACCCAC‑GGCCGGCGGGAAG‑GGTGTGCC‑GGGTTAACAACGTTTCTGGCGATGGTCTACCC‑GGTCATCGTCTTTCCAGGTATTT‑TGGGT‑AAACCAGGCTTCCCCCCTTCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTTCTGGACTCGGTTCTTTCGGGATTGGTCTGTGGGCTAATCTTCCGTTTGCGATTGGTTGCGCCATCTCCCTGACAGc                                                                                                                                                                             <  2:248668/238‑1 (MQ=255)                                                                  cACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGGTTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGttt                                                                                                                                                                         >  1:19152/1‑251 (MQ=255)                                                                     ccggcGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGGGTT‑TGCC‑GTGTTAACAACGTTTCGGGCGATGGTCTACTC‑GGTCATCGTCGTGCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTGCTATCGTGATGGGTCTGTGGTCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTAcc                                                                                                                                                                      <  1:172818/247‑1 (MQ=255)                                                                      cGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGGTTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTAc                                                                                                                                                                       <  1:294820/249‑1 (MQ=255)                                                                      cGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGGTTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCg                                                                                                                                                                     >  1:121579/1‑251 (MQ=255)                                                                       gACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGACGAGGC‑TGCC‑GGGTTAACAACGTTTCTGGAGATGGTATACTC‑GGTCATCGTCGTTCCAGGAATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑GTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGc                                                                                                                                                                    >  1:37233/1‑251 (MQ=255)                                                                         cGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCAGTAAGGCAGTGGT‑TGCC‑GGGTTCACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑ACAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTg                                                                                                                                                                                                          >  1:207389/1‑211 (MQ=255)                                                                             t‑GG‑TTTAAAATTTCACC‑ACG‑TGGAAGC‑ACTGTCCGTCCGGAAGTGGT‑TGCC‑GGGTTAACATCGTTTCTGGCGATGGTCTACTC‑GGTCATCGCCGTTCCAGGCATGT‑TGGGT‑AAAGAAGGCCTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTc                                                                                                                                                                >  1:225096/1‑249 (MQ=255)                                                                               gg‑TTTAACATTTCACA‑ACG‑TGTAAGC‑ACTGTCCGTCAGGAAGCGGT‑TACC‑GAATTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCTTCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGTCCGCTTTCGAGGCAACCTGTCTGG‑TTGCTGTACTCGGTTCTCTCGAGGTGTGACAGTGGGCTCATCTGCCGTTGGCGATCGGTGACGCATTCTTCCTGTCAGCGTTTACCGCAttaatg                                                                                                                                                             >  2:20078/1‑247 (MQ=255)                                                                               gg‑TTAAAAATTTCGCA‑ACG‑TGGAAGC‑ACTGTCCGTCAGAAAGTGGT‑TGCC‑GGGTTAACAAGGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGTTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGc                                                                                                                                                             <  1:119871/251‑1 (MQ=255)                                                                                g‑TTTAAAATTCCACA‑CCG‑TGGAAGC‑ACAGTCCGTAAGGGAGTGGT‑TGCC‑GGCCTAACAACGTTACTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGACTGCGGCGGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGCCAGCGTTTACCGCATTCCGcc                                                                                                                                                            >  2:58050/1‑251 (MQ=255)                                                                                           tCACA‑ACG‑TGGAAGA‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGGTTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTggg                                                                                                                                                  >  1:116562/1‑251 (MQ=255)                                                                                            caca‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGTT‑TGCC‑GGTTTACCAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTggg                                                                                                                                                  <  1:79131/250‑1 (MQ=255)                                                                                            caca‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGGTTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCctggggctgggg                                                                                                                                                 >  2:406572/1‑251 (MQ=255)                                                                                                 aCG‑TGGAGGC‑ACTGTCCGTCTGAGAGTGTG‑TGCC‑GGTTTAACAACGTTTCGTGCGATGGTCTACGC‑GGTAATCGGCGTTCCAGGAATGT‑TGGGT‑AAAGCAGGCTTCCCGCCGGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACGCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGTGGCAAc                                                                                                                                             <  2:183265/251‑1 (MQ=255)                                                                                                  cG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGGTTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACa                                                                                                                                            <  1:50005/251‑1 (MQ=255)                                                                                                     gggAACC‑ACTGTCCGTCAGGAAGGGGT‑GGCC‑GGTTTAACAACGTTTCTGGCGATGGTC‑GCTCCGGTCATCGGCGTTCCAGGTA‑GTGTGGGT‑AAGGCAGGTTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑GTGCTGGACTCGGTGCGATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTTGGGCAACata                                                                                                                                          <  2:390727/250‑1 (MQ=255)                                                                                                         aGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGGTTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATAtt                                                                                                                                        >  1:184515/1‑249 (MQ=255)                                                                                                                   cGTCAGGAAGTGGT‑TGCC‑GGGTTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTCTGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGt                                                                                                                             >  1:390073/1‑251 (MQ=255)                                                                                                                            gTGGT‑TGCC‑GGGTTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCa                                                                                                                                                               >  1:159854/1‑208 (MQ=255)                                                                                                                                    cc‑GGGTTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCgtg                                                                                                             >  2:214290/1‑251 (MQ=255)                                                                                                                                         gTTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑ATGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTcc                                                                                                         >  1:302421/1‑251 (MQ=255)                                                                                                                                         gTTAACAACGTTTCTGGCGATGGTCGTCTC‑GGTCATCGTCGTTCCAGGTTTGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTTTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTcc                                                                                                         <  1:224754/251‑1 (MQ=255)                                                                                                                                            aaccacGTTTCTGGCGAGGGTCTACTC‑GGCCAGCGCCGTTCTAGGTATGT‑TGGTG‑ACAGTAGGTTTCCCGCCTGCGGCAGCTTGCGTTGCACCCTGTCTGG‑TTGCGGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTGGGCGATTGGTTGCGCCATCTCCCATACAGCGTTTACCGCATTCAGCCCGGTGCTGGGGCAACATAGTATCGTGCCTGTCGCGCTGGGTGCCGTGTTACTGa                                                                                                      <  2:413690/247‑1 (MQ=255)                                                                                                                                                     tCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCCTCTCCCTGACAGCGTTTACCGCATTCa                                                                                                                                                               >  2:351904/1‑185 (MQ=255)                                                                                                                                                     tCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCa                                                                                                                                                               <  1:351904/185‑1 (MQ=255)                                                                                                                                                           gATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGttt                                                                                        >  2:212986/1‑250 (MQ=255)                                                                                                                                                                      c‑GGTCATCGTCGTTACAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCCACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGGTGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGGTTACCGCATTCAGCCTGGCGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGGGTTCCTGATGGGGGTGCTGTTTCCGGTAATTTCt                                                                            >  2:311905/1‑251 (MQ=255)                                                                                                                                                                             tcgtcgTTCCAGGTATGT‑GGGTT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCGGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCAc                                                                       <  1:247186/250‑1 (MQ=255)                                                                                                                                                                                  gtTCCAGGTGT‑T‑TGGGTTAAAGCTGGTTCCCCGCCTGCGGCAGCTTGCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGGTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGTTGCTGGGGCAACAGATTAGCGTACCTGGCGCGCTGGGTGCCGTGTGCCTGATGGGGGTGCTTTTTACGATAATTTCTGCCACGGGTa                                                                  <  2:385898/250‑1 (MQ=255)                                                                                                                                                                                           aTGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTg                                                        >  1:303182/1‑251 (MQ=255)                                                                                                                                                                                                         gcggggttCCCGCCTGCGGCAGCGTTCGTGGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGGGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTAGCCGTAGCTGGATTTTGCGCaa                                            <  2:353759/245‑1 (MQ=255)                                                                                                                                                                                                                  ccGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTc                                    >  1:62536/1‑250 (MQ=255)                                                                                                                                                                                                                           gCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGgc                          <  1:236223/251‑1 (MQ=255)                                                                                                                                                                                                                           gCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGgc                          >  1:142252/1‑251 (MQ=255)                                                                                                                                                                                                                           gCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGgc                          <  1:406572/251‑1 (MQ=255)                                                                                                                                                                                                                                    gTTGCAACCTGTCTGG‑CTACCGGACTCGGTTCTATCGTGCTGGGTCTGTGGGCTAATCTGCCGTCGGCGACTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTAGCGCTGGGGGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCTAAGGGTATCCGGAGCTGGATTTTGTGCAACTTGCCTCACGGTGTGGCgcccggcac                 >  2:268528/1‑251 (MQ=255)                                                                                                                                                                                                                                      tGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGGTGGGCCTGTGGGCTACTCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGAATTCAGCCTGGTGATGGGGCAACATCTTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGcacggcacgg               >  1:375137/1‑251 (MQ=255)                                                                                                                                                                                                                                        cAACCTGTCTGG‑TTGCTGGACTCGGTGCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACgggg             <  1:308734/251‑1 (MQ=255)                                                                                                                                                                                                                                          aCCTGTCTGG‑TTGCTGGACTCGGTTCGATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTTATGGGTGTGCTGTTTACGGTAAGTTCTGCCACGGGTATCCGTAGCTGGATTGTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGTGGAt           <  2:264854/251‑1 (MQ=255)                                                                                                                                                                                                                                              gTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGt       <  1:421542/251‑1 (MQ=255)                                                                                                                                                                                                                                                  gg‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATc    <  1:346292/250‑1 (MQ=255)                                                                                                                                                                                                                                                   g‑TTGCTGGAATCGGTTCCATCGTGATGGGCCTATGGACTCATCTGCCGATGGCGATTGGCTGCGCCATCTCCCTGCCAGCGTTTAACGCATTCAGCCTGGTGCTGGGGCCACATATTAGCGTACCAGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGAGCTGTTTACGGTAATTTCTGCCCCGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGCTTGGTATCgg  >  2:25299/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      TGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGGTTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGG‑TTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGCGCACGGCACGGGGATTGGTATCGG  >  NC_000913/4278445‑4278927

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