Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	388,319	C→G	17.4%	R175R (CGC→CGG)	tauD →	taurine dioxygenase, 2‑oxoglutarate‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	388,319	0	C	G	17.4%	157.6 / 2.9	69	R175R (CGC→CGG)	tauD	taurine dioxygenase, 2‑oxoglutarate‑dependent
Reads supporting (aligned to +/- strand):  ref base C (41/12);  new base G (2/11);  total (43/26)
Fisher's exact test for biased strand distribution p-value = 6.26e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCC  >  NC_000913/388098‑388554                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 gtgACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGCGGTAAGCCCTGGGGCAGCGGGCTTACCGGCCATGAGGCGGCCTATTTTTCCTTCCGCCCGCGTCCGGGTGGGGTGCGGGGGGAGCATGGATTTCCGTAAACCTTTCCC‑GAATA‑ACAATAACGCAAAAC‑CCAGGAGGGACATCCACCCATGCGGGCGGCGGGCGG‑CGaaaa                                                                                                                                                                                                                                >  2:38600/1‑239 (MQ=255)      aTTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGGTACGCCCTTGGCCCGCCGTATTTTGGCCTTTGAGGCGGTCTCTCTTTCCTTCCCCCCGCCGCTGAGTGGGCTGGCGGCGGGGGCT‑GGTTTTCCTTAATCCATTCCGGGCTACCACATACC‑CAAAAC‑CCACGAGGAACAACAACCCTGGCCGGGGGCGGTG‑G‑CGaaaaaaacc                                                                                                                                                                                                                           >  2:372987/1‑234 (MQ=255)                                                                     ttgttctctctGGGCCGGGGGTTTTGCGCCTTATGGGGCGCTTTCTGTCCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGGGGAGCAT‑GATTCCCGTAAATCGTCCCCGGAATA‑CAAATCCCGAAAAAC‑CGAGGAGAAACACCAACGCGGGCGGGAGGGGGTC‑G‑CGAAAAC‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGtgaatg                                                                                                                                                         <  2:1299756/232‑1 (MQ=255)                                                                       agcgcctctGGACCAGGTGTTGTGAGGGCGTTGGGGCGCTCTCTTTTCCCTTCGCCCGGCTGTGGGGGGGGTCGGGCGGGGATCAT‑GTTTCCCGTAAATCGTTCCGGAAAAA‑CAAAGACCGAAAACC‑CGGGGAGGAACAACAACCGTGGCGAGAGGCGGTC‑G‑AGAAAAC‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGtgaatga                                                                                                                                                        <  2:1219549/233‑1 (MQ=255)                                                                               tGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACAACCGGTGGTGCGAACGCATCCGGTGACCGGTAA‑ACCGGCGCTGTTTGTGAATGAAGGCTTtat                                                                                                                                               >  1:256765/1‑238 (MQ=255)                                                                                    cgggggTTTTGGGGCCTTGGGGGCCTTCTTTGTTCCCCTCCGCAAGTGGCTGGGGGGGCTGGGGGGGGAGCTT‑GATTCCCGAAAATCGTCCCCGGAAAA‑CAAATCCCGAAAGAC‑CGAGGAGGAACATAAACGCTGGGGCGGGGCGGTC‑G‑CGAAAAC‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTAcgc                                                                                                                                          <  2:353209/235‑1 (MQ=255)                                                                                      gCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCt                                                                                                                                                                                                                  >  2:810553/1‑166 (MQ=255)                                                                                      gCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCt                                                                                                                                                                                                                  <  1:810553/166‑1 (MQ=255)                                                                                       cGGTATTGCGGCCTATGAGGCGCTCTTTGTTCCCTCCCGCCGGCTGCTGAGTGGGCTGCGGGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATCCCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGt                                                                                                                                                                                       <  1:681528/192‑1 (MQ=255)                                                                                       cGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGACAAA‑CCCGCCGTTGCTAACTCCGGGGGTGCGAAAGCATCCGGGGAGCGGTAA‑AAAGGCGCGGTTTGTGAATGGAGGCTTtaacaggggaa                                                                                                                                       >  2:1279787/1‑230 (MQ=255)                                                                                       cGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGATCGGTAA‑ACAGGCCGTGTTTGTGAATGAAGGCTTTACAAcgttct                                                                                                                                       >  1:1219549/1‑235 (MQ=255)                                                                                           aTTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCCAACGCATCCGGTGAGGGGTAACAAAGGC‑CGGGTTTTGGATGAAGGCTTtactacccacattgg                                                                                                                                   >  2:561400/1‑230 (MQ=255)                                                                                           aTTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGCTGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGGGCTGTTTTGGAATGAAGGCGTtaatatgtgcaatgt                                                                                                                                   >  1:453391/1‑229 (MQ=255)                                                                                             tGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTtgt                                                                                                                                                              >  1:124894/1‑210 (MQ=255)                                                                                                    tATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAg                                                                                                                          >  1:1311401/1‑239 (MQ=255)                                                                                                    tATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAg                                                                                                                          >  1:1159711/1‑239 (MQ=255)                                                                                                      tGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGACAAA‑CCAGCCGTTGCTACATCCGGTGGGGCGAACGCAACCGGTGAGGGGTAA‑ACAGGCGCTGTTTTTGTATGAAGGCTTTACACGGCGAAttgttgttggtggcg                                                                                                                        >  2:609194/1‑233 (MQ=255)                                                                                                                  ttttCCCTCCCGCCGGGGGGTGGGTGGGGTGGGGGCGGAGCAT‑GTTTACCGAAAATCGTACCCGGAACA‑CAAATCCCGAAAAAC‑CGGGGAGGAACACCAACGGGGGCGGGGGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCg                                                                                                           <  2:569888/237‑1 (MQ=255)                                                                                                                    ttCCCTTCCGCCAGCTGCTGAGTGGGCTGCGGGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGAAATA‑CAAATCCCGCAAAAC‑CGAGGAGGAACATCAACGCGGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGaa                                                                                                         <  1:375749/239‑1 (MQ=255)                                                                                                                         ttCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGGCtg                                                                                                    >  1:548005/1‑238 (MQ=255)                                                                                                                          tCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAA‑AGAGCGAAGCCTTg                                                                                                   >  1:732532/1‑239 (MQ=255)                                                                                                                            cGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTTAGCGGGAAAA‑GGAGCGAAGCCTTGt                                                                                                  >  1:569888/1‑238 (MQ=255)                                                                                                                            cGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑GGAGCGAAGCCTTGtt                                                                                                 >  1:123999/1‑239 (MQ=255)                                                                                                                            cGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTAAGCGACTTGTTTATGTGGGCGGGAAAG‑AGAGCGAAGCATtttt                                                                                                 >  2:999807/1‑239 (MQ=255)                                                                                                                                 tcttgtgggTGGGGTGCGGGGGGGACAG‑GTTTTCCGTAAATCGGTCCCGGAATA‑CAAATACCGAAAAAC‑GGAGGAGGAACATCAACGCTGGCGCGGGGCGGGC‑G‑CGAAAAC‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGtt                                                                                            <  2:1159711/231‑1 (MQ=255)                                                                                                                                 gctgctGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAACTCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGGAGGCTTTACTACGCGAATTGTTGAGGTGAGCGGGAAAG‑AGGGCGAAGCCATGTTAAGtt                                                                                            >  2:915643/1‑239 (MQ=255)                                                                                                                                 gctgctGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAggt                                                                                            >  1:1062139/1‑237 (MQ=255)                                                                                                                                               tgcgtgcgGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATCCCGCAAAAC‑CGAGGCGGAACATCAACGCTGGCGCGAGGCGGCC‑G‑CGCAAAA‑CCCGCCGTTCCTACATCCGGTGGTGGGACCGTATCCGGTGAGCGGGAA‑TCCGGCCCCGTTTCTGAATGAACGTTTTACTACGCGCATTTCTGCTGTGTGCGTGCAGC‑AGGGGGAAGACTTGTTAAATTTTTTTTTTGtccca                                                                              >  2:271188/1‑237 (MQ=255)                                                                                                                                               tgcgtgcgGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGGGCGAACGCATCCGGGGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCat                                                                              >  2:1008279/1‑239 (MQ=255)                                                                                                                                                 cgtgcgGAGCCT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGTTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATACGGGGATGCGAACGCATCCGGTGCGCGTGTA‑ACAGGCGCCGTTGGTGAATCAAGGCTTTACTATGGTAATTGTTGAGGTGGGGGGGGAAG‑GAAGGGAGGCGATCTTAAGTTTTTTGTTTGACAatac                                                                            >  2:1242190/1‑238 (MQ=255)                                                                                                                                                  gtgcgGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGGGAATGAAGGCTTTACTACGCGAATTGTTGATGTGGGGGGGAAAG‑AGAGGGAAGCCTTGTTAAGTTTTTTTTTTGCCCAAATc                                                                           >  2:315456/1‑239 (MQ=255)                                                                                                                                                  gtgcgGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGCGAGCGGTAA‑ACAGGCGCTGCTTGTGAATGCAGGCTTTACTACCGGACATCGTGATGTGAGCGGGAAGGAAGAG‑GAAGCCTTGTTAAGTTTTTTTTTTGCCCaaaaa                                                                           >  2:420968/1‑235 (MQ=255)                                                                                                                                                           aT‑GATTCCCGAAATTCGTTCCGGGAATA‑CAAATCCCACAAACC‑GGAGGAGAAACATCAACGCGGGGGCGAGGCGGTC‑G‑CGAAAAC‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAAcc                                                                   <  1:315456/238‑1 (MQ=255)                                                                                                                                                           aT‑GATTCCCGAAAATCGTCCCCGGAAAC‑CAAATCCCGAAAACC‑CGAGGAGGACCATAAACGCTGGCGCGAGGCGGCG‑G‑CGAAAAC‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAAcc                                                                   <  1:322152/238‑1 (MQ=255)                                                                                                                                                                  tCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CAAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCCTCCGGTGAGCGGTAA‑ACCGACGCTGTTTGTGAGTGAAGGCGTTACTACGCGAATTGTTGATGTGAGCGAGAAGG‑AGAGCGAGGCCTTGTTGAGTTTTTTGTGTGCCCAAATCGCCAAACCGGAGttg                                                            >  2:1058005/1‑238 (MQ=255)                                                                                                                                                                     gaaaaTGGTTCCGGGAATA‑CAAATCCCGAAAACC‑CGGGGAGGAACATAACCGCGGGGGGGGGGGGGTC‑G‑CGAAAAC‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAg                                                         <  1:561400/237‑1 (MQ=255)                                                                                                                                                                     gTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTTAGCGGGAAAG‑AGACGGAAGCCTGGTTAATTTTTTTTTTTGCCCATATAACCAAACAGGGGTTTCAg                                                         >  2:84022/1‑239 (MQ=255)                                                                                                                                                                              tCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTgcgctg                                                 >  1:943716/1‑238 (MQ=255)                                                                                                                                                                              tCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGGgcgctg                                                 >  1:1354955/1‑238 (MQ=255)                                                                                                                                                                                ccGGAAAA‑AAAATCCCACAAAAC‑CG‑GGGGGAAAAACAAAGCGGGGCGGAGGGGGGC‑C‑CGAAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTgcgctggcg                                              <  2:1354955/239‑1 (MQ=255)                                                                                                                                                                                     aca‑CAAAGCCCGAAAAAC‑GGGGGAGAAACCACCACTGTGGCGAGGGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGcgctggcgctgg                                          <  1:1206123/236‑1 (MQ=255)                                                                                                                                                                                     aTA‑CAAATACC‑CAAAACGCGGGGAGAACCAGAAACGCTGGCGCGAGGGGGTC‑G‑AAAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCgctggcgctggc                                         <  1:925972/239‑1 (MQ=255)                                                                                                                                                                                     aTA‑AAAATACCGCAACAC‑GGGGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CAAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGcgctggcgctgg                                          <  1:356670/238‑1 (MQ=255)                                                                                                                                                                                             tcccccAAAAC‑ACAGGAGGCACATCAAGGGTGGCGGGAGGCGGTC‑G‑GAAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTCGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCaaa                                  <  1:1058005/234‑1 (MQ=255)                                                                                                                                                                                               ccGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTAc                                                                                                                                                                                                                >  2:1095068/1‑65 (MQ=255)                                                                                                                                                                                               ccGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTAc                                                                                                                                                                                                                <  1:1095068/65‑1 (MQ=255)                                                                                                                                                                                               ccGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGGGCTGGCAACAAAATg                                >  1:431563/1‑239 (MQ=255)                                                                                                                                                                                               ccGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGGGGCCTGGGCCTGGCAACCaaaga                                >  2:661588/1‑237 (MQ=255)                                                                                                                                                                                                  cAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGGGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCCCTGGCCACCAAACGata                             >  1:614804/1‑239 (MQ=255)                                                                                                                                                                                                  cAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGCAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGCCGCTGGCAACCAAATGata                             >  1:217565/1‑239 (MQ=255)                                                                                                                                                                                                  cAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCCAACCGGAGTTTCAGGTGCGTTGGCGCTGGCAACCacacaaaa                             >  1:1037688/1‑232 (MQ=255)                                                                                                                                                                                                  cAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTAATACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGGGTTTCAGATGCGCTGGCGCTGGCAACCAATTGaaa                             >  1:15332/1‑237 (MQ=255)                                                                                                                                                                                                     aaC‑CGGGGGGGAAAACAAACGGTGGGGGGGGGGGGGC‑GCCGAAAAC‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATAtt                           <  2:943716/239‑1 (MQ=255)                                                                                                                                                                                                      cc‑cGGGGGGGAACATCCACGGTGGGGGGGGGGGGGG‑G‑GGAAAAACCCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGGGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTg                          <  2:123999/238‑1 (MQ=255)                                                                                                                                                                                                       c‑cGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGGGCTGGCAACCAAATGATATTGCg                        >  1:774518/1‑239 (MQ=255)                                                                                                                                                                                                                gaACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGCGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGCGTTTCAGGTGCGCTGGCGCTGGCAACCAAATAAAATTTGGATTTGGGa                >  2:538144/1‑239 (MQ=255)                                                                                                                                                                                                                gaACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGa                >  1:946928/1‑239 (MQ=255)                                                                                                                                                                                                                  aCATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATa              >  1:87625/1‑239 (MQ=255)                                                                                                                                                                                                                        cctgtGGCGTGGGGGGGTC‑G‑CAAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACcgcg        <  1:372987/235‑1 (MQ=255)                                                                                                                                                                                                                        aggggggCGGGAGGGGGTC‑G‑CAAAAAC‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACcgcg        <  1:999807/234‑1 (MQ=255)                                                                                                                                                                                                                        accggggCGGGAGGCGGTG‑T‑CAAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTTAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACcgcg        <  1:735352/234‑1 (MQ=255)                                                                                                                                                                                                                        aCGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGGTTTACTACGAGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGGTTTTTGTTTGCCCATATCCCCAAACCGGAGTTTCAGGTACGCTGGCGCTGGCAACCAAATGATAATCCGATTAGGGAGAACcgcg        >  2:244623/1‑239 (MQ=255)                                                                                                                                                                                                                        aCGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGGGACGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGGCGGCTGGCGGTGGCACCCAAAAGATATTGCGATTTGGGAACACcgcg        >  2:343179/1‑239 (MQ=255)                                                                                                                                                                                                                        aCGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCCCCAAACCGGAGTTTCAGGTGCGCTGGGGCTGGCCAACCAATTCTATTGCGGTTTggggacaccccg        >  2:375869/1‑230 (MQ=255)                                                                                                                                                                                                                          gcgggCGCGAGGCGGAC‑G‑CG‑AACA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACcgcg        <  1:38600/233‑1 (MQ=255)                                                                                                                                                                                                                          gCTGGGGGGAGGCGGTC‑G‑AAAAAAC‑CCCGCCGTTGCTACATCCGGTGGGGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACcgcg        <  1:112170/237‑1 (MQ=255)                                                                                                                                                                                                                           tgggCGGGGGGGGGTC‑G‑GAAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACcgcg        <  1:169983/234‑1 (MQ=255)                                                                                                                                                                                                                           cTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCCCCAAACCGGAGTTTCAGTGTCGCTGGGGCTGGCAACCAAATAAAATTGCGATTTGGGTAAAcccggga      >  2:288436/1‑233 (MQ=255)                                                                                                                                                                                                                            tGGCGCGAGGGGGTC‑G‑CAAAAAA‑CCCGCCGTTGCTACATCCGGGTGTGCGAACGCAGCCGGGGAGGGGTAA‑ACTGGCGCTGTTTGTGAAGGAAGGCTTTACTACGCGAATTGTTGATGTGATCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACcgcg        <  1:673876/235‑1 (MQ=255)                                                                                                                                                                                                                              gcgcGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGCATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTTCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGAACCCAAATGATATTGCGATTTTGGAAAACCGCGGGGccc  >  2:1017803/1‑239 (MQ=255)                                                                                                                                                                                                                              gcgcGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGAt                                                                                                                               >  1:1203118/1‑115 (MQ=255)                                                                                                                                                                                                                              gcgcGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGCAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGGGTTTCAGGTGGGCTGGCGGCGGGAACCAAATGGTATTGCGAATTGGGATAAACGCGTGAcca  >  2:1218870/1‑238 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 GTGACATTTATTGAAACGCCACCCGCAGGGGCGATTCTGGCAGCTAAAGAGTTACCTTCGACCGGCGGTGATACGCTCTGGACCAGCGGTATTGCGGCCTATGAGGCGCTCTCTGTTCCCTTCCGCCAGCTGCTGAGTGGGCTGCGTGCGGAGCAT‑GATTTCCGTAAATCGTTCCCGGAATA‑CAAATACCGCAAAAC‑CGAGGAGGAACATCAACGCTGGCGCGAGGCGGTC‑G‑CGAAAAA‑CCCGCCGTTGCTACATCCGGTGGTGCGAACGCATCCGGTGAGCGGTAA‑ACAGGCGCTGTTTGTGAATGAAGGCTTTACTACGCGAATTGTTGATGTGAGCGAGAAAG‑AGAGCGAAGCCTTGTTAAGTTTTTTGTTTGCCCATATCACCAAACCGGAGTTTCAGGTGCGCTGGCGCTGGCAACCAAATGATATTGCGATTTGGGATAACCGCGTGACCC  >  NC_000913/388098‑388554

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: ‑