Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,275,602	T→C	18.1%	V108A (GTC→GCC)	garD →	D‑galactarate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,275,602	0	T	C	18.1%	147.4 / 4.1	71	V108A (GTC→GCC)	garD	D‑galactarate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (51/4);  new base C (2/13);  total (54/17)
Fisher's exact test for biased strand distribution p-value = 4.60e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGGCCTGAAAGCAGGAACGCGTTTTCCGGATGGGCTGGAATTAATTGAACATATTCCCCA‑GGGGCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCG  >  NC_000913/3275370‑3275837                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             aaTGGCCTGAAAGCAGGAACGCGTTTTCCGGATGGGCTGGAATTAATTGAACATATTCCCCA‑GGGGCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCCCTGGCAACCAAAGTCCCgg                                                                                                                                                                                                                                         >  1:31729/1‑238 (MQ=255)   tGGCCTGAAAGCAGGAACGCGTTTTCCGGATGGGCTGGAATTAATTGAACATATTCCCCA‑GGGGCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGaa                                                                                                                                                                                                                                       >  1:71597/1‑238 (MQ=255)          aaaGCAGGAACGCGTTTTCCGGATGGGCTGGAATTAATTGAACATATTCCCCA‑GGGGCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑c                                                                                                                                                                                                                               >  1:1101925/1‑238 (MQ=255)               aGGAACGCGTTTTCCGGATGGGCTGGAATTAATTGAACATATTCCCCA‑GGGGCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑Ccgccgc                                                                                                                                                                                                                         >  1:75796/1‑239 (MQ=255)                    cgcgTTTTCCGGATGGGCTGGAATTAATTGAACATATTCCCCA‑GGGGCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAC‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑CCACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑ACGCCGCGGGaa                                                                                                                                                                                                                    >  1:1234858/1‑239 (MQ=255)                    cgcgTTTTCCGGATGGGCTGGAATTAATTGAACATATTCCCCA‑GGGGCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACCCACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCCCTGGaa                                                                                                                                                                                                                    >  1:101193/1‑239 (MQ=255)                          gttcggggtgggttgAATTAATTGAAAAATTTCCCGAGGGGGAAAAAAGTCGCTGGGGTGCACGTCGCGGCAAGTGGTAAAATTATTGTTAATGGCG‑AGGGGTTGGTTCGCGGGGGGCGTG‑CAACCCCGCGGGGGAGGTGAACCGAC‑AAACAA‑TGGGTAGTACTCCCGGAAGCCCCGCCGT‑C‑ACACCCGCTGCCACTGGCAACAAAGGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATa                                                                                                                                                                                                               <  2:565029/225‑1 (MQ=255)                            ccGGATGGGCTGGAATTAATTGAACATATTCCCCA‑GGGGCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑AAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAACCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGCACCCTTC‑CCGCCGCTGGAAGCATaca                                                                                                                                                                                                             >  1:587770/1‑238 (MQ=255)                                              attGAACATATTCCCCA‑GGGGCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCaa                                                                                                                                                                                          >  1:427438/1‑239 (MQ=255)                                                  aaCATATTCCCCA‑GGGGCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATCAACCTTTGAGGGCTATCGCAATGcc                                                                                                                                                                                      >  1:222987/1‑239 (MQ=255)                                                       atTCCCCA‑GGGGCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAC‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACCCACGCTGCCACTGGCAACCACAGCCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGAGCTATCCCAATGCCGATgg                                                                                                                                                                                 >  1:408229/1‑239 (MQ=255)                                                             cA‑GGGGCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCg                                                                                                                                                                                                                                          >  1:76305/1‑177 (MQ=255)                                                             cA‑GGGGCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCGGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCg                                                                                                                                                                                                                                          <  2:76305/177‑1 (MQ=255)                                                                ggggCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTgg                                                                                                                                                                         >  1:954160/1‑239 (MQ=255)                                                                ggggCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGCAAGCGTgg                                                                                                                                                                         >  1:315757/1‑239 (MQ=255)                                                                ggggCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCCCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCgggg                                                                                                                                                                         >  1:1351884/1‑239 (MQ=255)                                                                ggggCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCCCCGCTGGAAGGACAAACCTTTGAGGGCTACCCAAATACCGAGGGCAGCgggg                                                                                                                                                                         >  2:138801/1‑239 (MQ=255)                                                                 gggCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGCATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCCGCCACTGGCCACCAAAGTCCCGGACCCCTTA‑CCGCCCTCCGACGGATCCCCCTTTTGGGGCCATCGCAATCCCGAGGGCCGCgggg                                                                                                                                                                         >  2:157327/1‑238 (MQ=255)                                                                 gggCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACCCCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTgg                                                                                                                                                                         >  1:380057/1‑238 (MQ=255)                                                                 gggCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTggg                                                                                                                                                                        >  1:338201/1‑239 (MQ=255)                                                                 gggCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCCCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCgggg                                                                                                                                                                         >  1:912551/1‑238 (MQ=255)                                                                 gggCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑AAACACGCTGCCACAGGCAAACAAAGTCCCGGAACCCTTA‑CCCCCGCTGGAAGGATAAACCCTTGAGGGCTATCGAAATggccaagggcggggggg                                                                                                                                                                        >  2:983508/1‑223 (MQ=255)                                                                    cATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACCCGCTGCCACTGGCAACCAAAGTCCCCGAACCCTTA‑CCGCCGCTGGAAGGATACCCCTTTGAGGGCTATCGAAATGCCGATGCCCGCGTGGGCa                                                                                                                                                                      >  1:783139/1‑238 (MQ=255)                                                                             gCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACCCCTTTGAGGGCCATCGCAATGCCGATGGCGGCGGGGGGACCAAAAAcat                                                                                                                                                            >  2:596930/1‑237 (MQ=255)                                                                             gCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATAACCCTTTGAGGGCTATCCCAATGCCGATGGCAGCGGGGGAACCAAAAAcc                                                                                                                                                             >  1:799915/1‑238 (MQ=255)                                                                              cATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGTTGAATCAC‑T‑GGTTGTACTACCCGAAGCGCCGCCGT‑T‑AGACACGTTGCCACTGGCCACCAAAGTCCCGGCACCCTTA‑CCGCCGCCGGAAGGATACCCCTTTGAGGGCTATCGCAATGCCGATTGcgcaggaggtagcacaaacct                                                                                                                                                            >  2:209113/1‑218 (MQ=255)                                                                              cATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGCTGCCAACGGGGGCACCAAAAActc                                                                                                                                                            >  1:1352550/1‑236 (MQ=255)                                                                                     ggacttctcGGGTAATGGTGAAATTTTTCGGTGTGGGGATGGGTTTGGTT‑ACGGCG‑GTGTGCCCACCCCCAGGGGAGAGTGGAGCCGAGAAACAA‑A‑TGTTGTCCCTCCCGAAGCCCCCCCCT‑TAACCCCCTCTGCCATGGCCAACCAACGCCCCGGCACCCTAA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGt                                                                                                                                                     <  2:1322423/231‑1 (MQ=255)                                                                                           tCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGGTCGACGAATCAA‑T‑GGTTGTACTACCTGAACCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGCCCCGGAACCCTTA‑CCGCCGCTGGCAGCATACACCTTTGGGGGCTCTCGAAAAGCCGATGGCAGCGGGGGCACCAAAAACCTgactggtaccaac                                                                                                                                               >  2:396512/1‑226 (MQ=255)                                                                                           tCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGCGGGCACCAAAAACCGGCTCGGAATcacc                                                                                                                                               >  1:1148128/1‑238 (MQ=255)                                                                                           tCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCAACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGCGGGCTATCGCAATGCCCATGGCAGCGGGGGCACCACAAACCTGCCCGGTATcacc                                                                                                                                               >  1:15662/1‑238 (MQ=255)                                                                                                     ggTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTcc                                                                                                                                     >  1:1113530/1‑238 (MQ=255)                                                                                                                tCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGACCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCCCAATGCCGATGCCAGCGTGGGAACCAAAAACCTGCCCGGGATCCCAACCAGCG‑CCCCAGTGTgggacg                                                                                                                        >  2:15544/1‑235 (MQ=255)                                                                                                                     aTGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTCCTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCACAAACCTGCTCGGTATCACAACCAGCGTCCACTGGGTGGCAGGCgtg                                                                                                                    >  1:189212/1‑239 (MQ=255)                                                                                                                       ggCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCaagggggcg                                                                                                                  >  1:475034/1‑231 (MQ=255)                                                                                                                           aaGTGATTGTTT‑CCGCCGTGCGTG‑CACTCCCCGGCGGAAGCGGGACCGACAAATAAA‑T‑GGTTGTCCTACCGGAAGCGCCTCCGT‑T‑ACACACGCTGCCACGGGAAACCAAAGTCCCGGACCCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGa                                                                                                               <  1:330114/238‑1 (MQ=255)                                                                                                                            aGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGa                                                                                                                                                                                                      >  1:223263/1‑150 (MQ=255)                                                                                                                             agtgggtgtttgcGCCGTGCGTG‑ACCCCCCACGCGGAAGCTGGGAAGAAGTAACAATT‑GGTGGTACACCGGGACGGCCCGCGGT‑C‑ACCCGCGCTGCCACTGGCAACAAAGGTCCCGAACCCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGa                                                                                                               <  2:799915/226‑1 (MQ=255)                                                                                                                                      t‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGCTGGCAGCGTGGGCCCCAAAAACCTGCCCGGTATCCCCCCCAGCGCCCCCGGTGTGGCCggggggggggcctcgatagaaaa                                                                                                   >  2:963007/1‑218 (MQ=255)                                                                                                                                         cGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGCCTATGTAGTacaa                                                                                                  >  1:1143189/1‑238 (MQ=255)                                                                                                                                               ccgtg‑cAATCCCACGCGGAAGCTGGATCGACGAGTCAA‑T‑GGGTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTACCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGCGACACCTTTGATGGCAGTGGCAATGCCAAGGGCAGCATGGGTACCAAAAACCGGCTCGGGATAACCACACGCGCCCACTGTGGGGAGGGCGgggggggcacggaagaaaaaagagatg                                                                                           >  2:904381/2‑218 (MQ=255)                                                                                                                                                cgtg‑cAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGAATCACCACCAGCGTCCACTGTGGGGCAGGCGGGGGGGACTATGTAGGAAAAACCATTGa                                                                                          >  1:748585/1‑239 (MQ=255)                                                                                                                                                             cgcgGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTTAGGGCTATCGCAATGCCGATGGCAGCGGGGGCACCAAAAACCTGCTCGGTATACCCACCAGCGTCCACTTTGTGGCAGGCGGGGTGGACTATGTGGGAAAAATCATTGAAcccggcctcc                                                                               >  2:966943/1‑232 (MQ=255)                                                                                                                                                             cgcgGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCCCGGTATCCCCACCAGCGTCCACTGTGTGGCAGGCGGGGTGGGCTAGTGGGTAAAAATCATTTAACCCGATTTGCt                                                                              >  2:896506/1‑239 (MQ=255)                                                                                                                                                             cgcgGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCCCGGTATCACCCCCAGCGTCCACTGTGTGGCAGGCGTGGGGGACTATGAAGTAAAAACCATTGAACCCGATcctct                                                                              >  1:661551/1‑235 (MQ=255)                                                                                                                                                                  aaGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTCTGAAGGCTATCGCAATGCCGCTGGCGGCGCGGGGACCAAAAACCTCCTCCGTATCCCCACCAGCGTCCACGTCAAGGCAGGCGTGGCGGACAATGTAATAAAAATCATTGAAATCGATCGGCCACCGa                                                                         >  2:1018355/1‑239 (MQ=255)                                                                                                                                                                        gATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTAACGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGGGGGGCAGGCGGGGGGGCCCATCTAGTAAAAAACATTGAAcccggacttccacccaaaaac                                                                    >  2:740988/1‑221 (MQ=255)                                                                                                                                                                          tCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCCAAAACCTTCTCGGTATCACCACCACCGTCCACTGTGTGGCAGGGGTGGGGGACTAGTGAGAAAAAAACATTGAACGCAGTCTGCCACCGAAAAACCCg                                                                 >  2:1002963/1‑239 (MQ=255)                                                                                                                                                                                             gTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTCAAAATCATTGAACGCGATCTGCTCCCGAAATACCCGAACGTCGATGGCgtggt                                                >  1:898654/1‑239 (MQ=255)                                                                                                                                                                                             gTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAAt                                                                                                >  1:967110/1‑191 (MQ=255)                                                                                                                                                                                                  cgcGGACGCGCCCCCCTCT‑ACACCCGTTGCCCGTGCAAACAAAAGCCCGGGATCCCGTA‑CCGCCGCTGGAAGCATACACCTTTGAGGGCAGTCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCGCCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAACCTACCGAACGTCGATGGCGTGGTgggg                                            <  2:222987/237‑1 (MQ=255)                                                                                                                                                                                                   ccGGAAGCGCCGCCGT‑T‑CCACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGGCTATGTAGTAAAAATCATTGAACGCGATCTGCTCCTGGAATACCCGAACGGCGATGGCGTGGCgggggg                                          >  2:259433/1‑237 (MQ=255)                                                                                                                                                                                                          cgccgcCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCCCTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAAAACCCAAACGTCGATggggggggggggtgggacata                                   >  2:494877/1‑230 (MQ=255)                                                                                                                                                                                                             cgccgT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGGGGGGCTGAATCATTTg                                >  1:1128196/1‑239 (MQ=255)                                                                                                                                                                                                             ccccgT‑T‑CCACCCGCCTCCCCTGGCAACCAAAGCCCCGGCACCTTTC‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTg                                <  1:983508/237‑1 (MQ=255)                                                                                                                                                                                                               cgcttttccccccccTGCCTCGGCCACCCAAAGCCCCGGACCCTTTC‑CCGCCGCTGGAAGGATACCCCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTg                                <  2:427438/225‑1 (MQ=255)                                                                                                                                                                                                               cctt‑t‑ACCCAGGCGGCCAATGCCAGCTAAAGCCCCGAGACCTGTT‑CCGGCGCTGGAAGGATACACCTTTGAGGGCTATCTCAATGCCGATGTCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTa                              <  1:1018355/236‑1 (MQ=255)                                                                                                                                                                                                               ccgT‑T‑ACCCCCGCTCCCACTGGCAACCAAAGCCCCGCACCCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTa                              <  1:1082778/239‑1 (MQ=255)                                                                                                                                                                                                                gatgtaaaaaaggaaaagtcTGGCAACCAAGGCCCCGTACCCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTg                                <  1:199207/219‑1 (MQ=255)                                                                                                                                                                                                                 ctttcacaccgccggccccgggCAACCAAAGCCCCGGACCCCTTC‑CCGCCGCTGGAAGGATACCCCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAATCCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTa                              <  1:1087697/219‑1 (MQ=255)                                                                                                                                                                                                                 ccttcccccccccccgcccTGGCACACAAGGCCCCGGACCCCTTC‑CCGCCGCGGGAAGGATACCCCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACTAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTa                              <  2:912551/221‑1 (MQ=255)                                                                                                                                                                                                                 ccttcaccccccctccccttGGCAACCAACCCCCCGGACCCCTTC‑CCGCCGCTGGAAGGATTCCCCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTa                              <  2:730051/220‑1 (MQ=255)                                                                                                                                                                                                                 cctccaccccctgcccccgtcgCAACCAACCCCCCGACACCCCTC‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTa                              <  2:1143189/218‑1 (MQ=255)                                                                                                                                                                                                                       acgcctagccccgtGGCAAACAAA‑CCCAGGAAGCCGTACCCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTtg                        <  1:904381/226‑1 (MQ=255)                                                                                                                                                                                                                            gCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGGGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAAAGTCGATGGCGGGGGGGGGCTGAACCATTTGGACGGTTgggggg                   >  2:1110833/1‑237 (MQ=255)                                                                                                                                                                                                                                 cgacggacacaaaaGCCCCGGACCCCTTA‑CCGCCGCTGGAAGGATACCCCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGtggcgtggcg              <  1:138801/228‑1 (MQ=255)                                                                                                                                                                                                                                      gCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGGTTaa         >  1:536241/1‑239 (MQ=255)                                                                                                                                                                                                                                        aaCCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGCCCGTGGCGGGGTGGGGCCGAATCATTTTTACGGTTGTGGCGGGGGGGTTaaaa       >  2:983855/1‑237 (MQ=255)                                                                                                                                                                                                                                        aaCCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGCGTGGCAGGCGAGGTGGACTATGTAGTAAAAATAATTGAACGCGCTCTGCTACCGCAATACCCGAACGTCGATGGCGGGGGGGGGGTGAAACATTTGTACGGGTgtgggggggggggaaaaa       >  2:305046/1‑225 (MQ=255)                                                                                                                                                                                                                                          ccAAAGTCCCGGAACCCTTA‑ACGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACGATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAAAACCCGAACTCCCACGGCGTGTGGGGGTTGAATCCTTTGTCAGGTTgtagggtgggggtataggca     >  2:228210/1‑228 (MQ=255)                                                                                                                                                                                                                                             aaGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGGGGGGGACTATGTAGTAAAAATCATTGAACGCGGTCTGCTACCGAAATACCCGAACGTCGAAGGCGTGGGGGGGCTGAAGCATTTGTACGGtttggggtgtcctttaaaaaaaaac  >  2:199207/1‑216 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             AATGGCCTGAAAGCAGGAACGCGTTTTCCGGATGGGCTGGAATTAATTGAACATATTCCCCA‑GGGGCATAAAGTCGCATTGCTGGACATTCCGGCTAATGGTGAAATTATTCGTTATGGCGAAGTGATTGGTT‑ACGCCGTGCGTG‑CAATCCCACGCGGAAGCTGGATCGACGAATCAA‑T‑GGTTGTACTACCGGAAGCGCCGCCGT‑T‑ACACACGCTGCCACTGGCAACCAAAGTCCCGGAACCCTTA‑CCGCCGCTGGAAGGATACACCTTTGAGGGCTATCGCAATGCCGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCG  >  NC_000913/3275370‑3275837

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