Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,278,657	T→C	V60A (GTT→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,657	0	T	C	98.1%	149.6 / ‑6.3	54	V60A (GTT→GCT)	ghxP	guanine/hypoxanthine transporter GhxP
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (29/23);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGG‑TTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGC  >  NC_000913/4278413‑4278903                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                        ttCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGACGCCT‑GG‑GTTAAAATTTCCCA‑ACG‑TGGGAGG‑CCTGTACGTCAGGGAGTGGT‑TGCC‑GGG‑TTAACAACGCTTATGGGGATGGTCTACAC‑GGTCATCGTCGGTAAAGGGATGT‑TGGGG‑AAAGCAGGAGTCCCGCCTGCGGACGTTTTCGt                                                                                                                                                                                                                                                   >  1:36080/1‑251 (MQ=255)                ttGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGACGCCT‑GG‑GTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGG‑TGCC‑GGG‑GTAACAAAGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCCGGGATGT‑TGGGG‑AAAGCAGGCTTTCCGCCCGCGGCAGCCTTCGTTTGAACCTGTCTggg                                                                                                                                                                                                                                    >  1:125884/1‑250 (MQ=255)                  gTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑CCTGTCCGTCAGGAAGTGGT‑TGCC‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCGTTCGTTGCAACCTGTCTGG‑TTg                                                                                                                                                                                                                                 <  1:343569/251‑1 (MQ=255)                            cATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCCACCCCT‑GG‑GTTAAAAATTTAAA‑AAA‑TGGAAGC‑ACCGTCCACCAGGAAGTGGT‑TGCC‑GGGGTTAA‑AAAGATTCTGGCGATGGTCTACTC‑GGTCCTCGTCGTTTCCGGTATGG‑TGGGG‑GAAGGAGGGTTCCCCGCTGCGGCCGGTTTCGTTGTCACCAGTCTGG‑TTGcgtggctcgg                                                                                                                                                                                                                       >  1:893/1‑242 (MQ=255)                                 tGACGAAAATTAAACCCTAGGGGGATTTTTTCTTATGTCTGCGCCATAAGCGCGTACCGGCGGTTCCCTCGACGCCT‑GG‑TTTAAAATTTCACA‑TCG‑TGGAAGC‑TCTGTCCGTCAGGAAGGGGT‑TGGC‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGTT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTc                                                                                                                                                                                                                    <  2:259379/249‑1 (MQ=255)                                     gAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGG‑TTAACAACGTTTCTTGCGATGGTCTACTC‑GTTCATCGTCGTTCCAGGTATGT‑GGGGT‑AAAGCAGGCTTCCCGTCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTTGACTCGGTTCTATCGt                                                                                                                                                                                                              <  1:88102/251‑1 (MQ=255)                                             aaCTCTCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑GGGGAAC‑ACTTTCCGTCAGGAAGGGGG‑TGCC‑GGG‑TTAACAACGTTTCTGGCGGTGGGCTACTC‑GGACATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCACGCCGGCGGCCGCTTTCGTTGCAACCTTTGTGG‑TTGCTGGGCTCGGTTCTATCGCGATggggc                                                                                                                                                                                                      >  1:353759/1‑249 (MQ=255)                                               ctctCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGGCGCCT‑GGGTTTAAAATT‑CACACACG‑TG‑AGGC‑CCTGTCCGCCCGGAAGTGGT‑TGCC‑GGG‑GTTACACAGTTTCTGGCGATGGGCTACTC‑GGTCATCGTCGGTCCCGGTCTGT‑TGGGT‑AAAGCAAGCGTCCCCCCTGGGGCAGCTTTCGTTTCCACCTGTCTGG‑TTGTCTGACTCGGGTCtttccgcattggtctg                                                                                                                                                                                                    >  1:353945/1‑236 (MQ=255)                                                   cAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTTACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCCGGAAGGGGT‑TGCCCGGG‑TTA‑CAACGGTTCTGGCGGTGGTCTCCTC‑GGTCATCGTCGTTCCCGGGATGT‑TGGGG‑AAAGCAGGCTTCCCCCCTGCGGGACCCTTCGTTGCCACCTTTCTGG‑GTGTTGGGCTCGGGTCTATCTTTATGGGGCTTTggg                                                                                                                                                                                                >  2:50005/1‑251 (MQ=255)                                                         aTGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTAACTCGAAGCCC‑GG‑TTTAAAATTTCACA‑ACC‑TGGAAGC‑ACTGTCCGGCAGGGAGTGGG‑TGCC‑GGG‑GTAACAACGTTTCTTGCGATGGTCTACTC‑GGTCATCGTCGTTTCAGGGATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGGCGCTTTCGGTGCAACCTGTCTGG‑TTGCTGGGCTCGGGTCTTTCGTGATGGGGCTGTGGGGTAATc                                                                                                                                                                                          >  1:146123/1‑251 (MQ=255)                                                                 ttATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTTGACCCCT‑AG‑ATTAAAATTTTACA‑ACGAGGGAGGC‑A‑TGTCCGTCAGGAAGCGGT‑TGCC‑GGG‑GTAAAAACGCTTCTTGGGCTGGTCTTAAC‑GGGCATAGTCGTTCCAGGGATGT‑TGGGG‑AAAGAAGGCTTCCCCCCTGCGGGCGCATTCGTTGCAAACCGGCTGG‑TTGCTGGAGTCGGGTCTATCGTGGTGGGGCTGTGGGGTAATCCGCCGTTg                                                                                                                                                                                  >  1:264854/1‑251 (MQ=255)                                                                   aTGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGAAGCCCCGG‑TTTAAA‑TTTCACA‑ACG‑AGGAAGC‑ACTGTACGTCAGGCAGTGGT‑TGCC‑CGG‑TTTACAACGTTTCTTGCGGTGGTATACTC‑GGCCATCGTCGTTCCACGAATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCCACACGTCTGGGTTGT‑GGGCTCGGTTCTATCGTGATGGGTCTGTGGGCTTATCTTCCGTTGGc                                                                                                                                                                                >  1:385898/1‑251 (MQ=255)                                                                             cATCAGCGCGTACCGGCGGTTCACCCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑CCTGTCCGTCAGGAAGTGTT‑TGCC‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTgc                                                                                                                                                                      <  2:125884/251‑1 (MQ=255)                                                                                     ggTACCGGGGGTTCACTCGACGCCT‑GG‑TTAAAAATTTCACA‑CTG‑TGGAACC‑ACTGTCCGTCAGGAAGGGGT‑TGCC‑GGG‑TTAAAAACTGTCCTGGCGTTGTTCTTCTC‑GGGCATCGTCGTTCCAGGTCGGT‑TGGGT‑AAAGCGGTCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCGGTCTGG‑TTGCTGGACTCGGTTCTGTCGTGATGGGTCTGTGTGCTAATCTGCCGTTGGCGATTGGTTGCGCCAtctc                                                                                                                                                              <  2:353945/250‑1 (MQ=255)                                                                                      gtacctgtgtttCACTCGACGCGT‑GT‑TTTAAAATTTCAAA‑GCG‑TGGAAGC‑ACTGTCCGTCAGGAGGTGGT‑TGCC‑GGG‑TTAACAACGTTTCTGGCGAGGGTCTACCC‑GGTAATCGTCGTTCCGGGTATGT‑TGGGT‑AAGGCGGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTcc                                                                                                                                                             <  2:189486/241‑1 (MQ=255)                                                                                         ccGGCGGTTACCTCGACGCCT‑GG‑TTTAAACTTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGG‑TTACCAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTg                                                                                                                                                          >  2:337485/1‑251 (MQ=255)                                                                                            gcggTTCACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGG‑TTAACAACGTTTCTGGCGAGGGTCTACTG‑GGTCATCTTCGTTCCAGGTATGT‑TGGTT‑AAAGCAGTCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGGGATGGGTTGCGCCATCTCCCTGACa                                                                                                                                                       <  1:126676/251‑1 (MQ=255)                                                                                             cggTTCACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCCGGAAGTGGT‑TGCC‑GGG‑CTAACAACGTTTCTGGCGATGGTCCACTC‑GGTCATCGTAGTTCCCGGCATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGCCTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGAc                                                                                                                                                        >  1:204576/1‑249 (MQ=255)                                                                                              ggttacctcggacccT‑GG‑TTTAAAATTTACCA‑ACG‑TGGAACCCAC‑GGCCGGCGGGAAG‑GGTGTGCC‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACCC‑GGTCATCGTCTTTCCAGGTATTT‑TGGGT‑AAACCAGGCTTCCCCCCTTCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTTCTGGACTCGGTTCTTTCGGGATTGGTCTGTGGGCTAATCTTCCGTTTGCGATTGGTTGCGCCATCTCCCTGACAGc                                                                                                                                                     <  2:248668/238‑1 (MQ=255)                                                                                                  cACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGttt                                                                                                                                                 >  1:19152/1‑251 (MQ=255)                                                                                                     ccggcGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGGGTT‑TGCC‑GTG‑TTAACAACGTTTCGGGCGATGGTCTACTC‑GGTCATCGTCGTGCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTGCTATCGTGATGGGTCTGTGGTCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTAcc                                                                                                                                              <  1:172818/247‑1 (MQ=255)                                                                                                      cGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTAc                                                                                                                                               <  1:294820/249‑1 (MQ=255)                                                                                                      cGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCg                                                                                                                                             >  1:121579/1‑251 (MQ=255)                                                                                                       gACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGACGAGGC‑TGCC‑GGG‑TTAACAACGTTTCTGGAGATGGTATACTC‑GGTCATCGTCGTTCCAGGAATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑GTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGc                                                                                                                                            >  1:37233/1‑251 (MQ=255)                                                                                                         cGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCAGTAAGGCAGTGGT‑TGCC‑GGG‑TTCACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑ACAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTg                                                                                                                                                                                  >  1:207389/1‑211 (MQ=255)                                                                                                             t‑GG‑TTTAAAATTTCACC‑ACG‑TGGAAGC‑ACTGTCCGTCCGGAAGTGGT‑TGCC‑GGG‑TTAACATCGTTTCTGGCGATGGTCTACTC‑GGTCATCGCCGTTCCAGGCATGT‑TGGGT‑AAAGAAGGCCTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTc                                                                                                                                        >  1:225096/1‑249 (MQ=255)                                                                                                               gg‑TTTAACATTTCACA‑ACG‑TGTAAGC‑ACTGTCCGTCAGGAAGCGGT‑TACC‑GAA‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCTTCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGTCCGCTTTCGAGGCAACCTGTCTGG‑TTGCTGTACTCGGTTCTCTCGAGGTGTGACAGTGGGCTCATCTGCCGTTGGCGATCGGTGACGCATTCTTCCTGTCAGCGTTTACCGCAttaatg                                                                                                                                     >  2:20078/1‑247 (MQ=255)                                                                                                               gg‑TTAAAAATTTCGCA‑ACG‑TGGAAGC‑ACTGTCCGTCAGAAAGTGGT‑TGCC‑GGG‑TTAACAAGGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGTTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGc                                                                                                                                     <  1:119871/251‑1 (MQ=255)                                                                                                                g‑TTTAAAATTCCACA‑CCG‑TGGAAGC‑ACAGTCCGTAAGGGAGTGGT‑TGCC‑GGC‑CTAACAACGTTACTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGACTGCGGCGGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGCCAGCGTTTACCGCATTCCGcc                                                                                                                                    >  2:58050/1‑251 (MQ=255)                                                                                                                           tCACA‑ACG‑TGGAAGA‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTggg                                                                                                                          >  1:116562/1‑251 (MQ=255)                                                                                                                            caca‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGTT‑TGCC‑GGT‑TTACCAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTggg                                                                                                                          <  1:79131/250‑1 (MQ=255)                                                                                                                            caca‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCctggggctgggg                                                                                                                         >  2:406572/1‑251 (MQ=255)                                                                                                                                 aCG‑TGGAGGC‑ACTGTCCGTCTGAGAGTGTG‑TGCC‑GGT‑TTAACAACGTTTCGTGCGATGGTCTACGC‑GGTAATCGGCGTTCCAGGAATGT‑TGGGT‑AAAGCAGGCTTCCCGCCGGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACGCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGTGGCAAc                                                                                                                     <  2:183265/251‑1 (MQ=255)                                                                                                                                  cG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACa                                                                                                                    <  1:50005/251‑1 (MQ=255)                                                                                                                                     gggAACC‑ACTGTCCGTCAGGAAGGGGT‑GGCC‑GGT‑TTAACAACGTTTCTGGCGATGGTC‑GCTCCGGTCATCGGCGTTCCAGGTA‑GTGTGGGT‑AAGGCAGGTTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑GTGCTGGACTCGGTGCGATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTTGGGCAACata                                                                                                                  <  2:390727/250‑1 (MQ=255)                                                                                                                                         aGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATAtt                                                                                                                >  1:184515/1‑249 (MQ=255)                                                                                                                                                   cGTCAGGAAGTGGT‑TGCC‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTCTGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGt                                                                                                     >  1:390073/1‑251 (MQ=255)                                                                                                                                                            gTGGT‑TGCC‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCa                                                                                                                                       >  1:159854/1‑208 (MQ=255)                                                                                                                                                                    cc‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑TTGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCgtg                                                                                     >  2:214290/1‑251 (MQ=255)                                                                                                                                                                         g‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGCTTTCGTTGCAACCTGTCTGG‑ATGCTGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTcc                                                                                 >  1:302421/1‑251 (MQ=255)                                                                                                                                                                         g‑TTAACAACGTTTCTGGCGATGGTCGTCTC‑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)                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                        TTCGCTGGCAGGGGATTGTCCGCCACGCATCTTGACGAAAATTAAACTCTCAGGGGATGTTTTCTTATGTCTACGCCATCAGCGCGTACCGGCGGTTCACTCGACGCCT‑GG‑TTTAAAATTTCACA‑ACG‑TGGAAGC‑ACTGTCCGTCAGGAAGTGGT‑TGCC‑GGG‑TTAACAACGTTTCTGGCGATGGTCTACTC‑GGTCATCGTCGTTCCAGGTATGT‑TGGGT‑AAAGCAGGCTTCCCGCCTGCGGCAGTTTTCGTTGCAACCTGTCTGG‑TTGCCGGACTCGGTTCTATCGTGATGGGTCTGTGGGCTAATCTGCCGTTGGCGATTGGTTGCGCCATCTCCCTGACAGCGTTTACCGCATTCAGCCTGGTGCTGGGGCAACATATTAGCGTACCTGTCGCGCTGGGTGCCGTGTTCCTGATGGGTGTGCTGTTTACGGTAATTTCTGCCACGGGTATCCGTAGCTGGATTTTGCGCAACTTGCCTCACGGTGTGGC  >  NC_000913/4278413‑4278903

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