Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,763	C→T	N10N (AAC→AAT)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,763	0	C	T	90.0%	92.8 / NA	41	N10N (AAC→AAT)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (11/25);  minor base G (0/4);  total (12/29)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAACGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGC  >  NC_000913/3525‑4006                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                         tCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCa                                                                                                                                                                                                                                                <  1:246053/249‑1 (MQ=255)                aggcgGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTa                                                                                                                                                                                                                              <  1:621000/251‑1 (MQ=255)                   cgGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGGCCGAGGACTGGGAAACTAAAAAAAA‑ATCTT‑AAATATGAAAAATCACAATGAGCAgggacacctttcgcaaaccgtaaac                                                                                                                                                                                                                            >  1:101763/1‑228 (MQ=255)                   cgGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑cc                                                                                                                                                                                                                          <  1:499070/251‑1 (MQ=255)                     gttgcgGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑a                                                                                                                                                                                                                       <  2:493145/248‑1 (MQ=255)                        gcgGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑Agg                                                                                                                                                                                                                     <  1:145299/250‑1 (MQ=255)                                   cgcgGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CGCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCaaaa                                                                                                                                                                                                         <  1:615222/251‑1 (MQ=255)                                    gcgttagcgtggggtttctcCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAGACCGCCCAGCGTGTTGACGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGCTCATATTTGCCGGCTGGATGCGTCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CTATCTGAAAGACCGCAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAT‑CAC‑AGGGGGTGGGCaaaaa                                                                                                                                                                                                        <  2:528377/235‑1 (MQ=255)                                                 gTATCTCCGGCTCCGGCCCGACTTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGGAA‑CCC‑AGGTGTTGGGCAAAAATCAGGGGCTGttt                                                                                                                                                                                           <  2:593496/251‑1 (MQ=255)                                                       ccggctccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑cc                                                                                                                                                                                     <  2:156231/250‑1 (MQ=255)                                                          gctccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGGGGAATGGAAAAATAAATGAAA‑CACTA‑CAATCTTGAAGATCACAAAAAGCAGAGCAG‑CTTTGCGCAAACCCTAA‑CCC‑AGGGGGTGGGGAAAAATCAAGGGATGGTTTTT‑TCGCAc                                                                                                                                                                                 >  2:410505/1‑251 (MQ=255)                                                             ccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACg                                                                                                                                                                                <  1:224215/249‑1 (MQ=255)                                                             cccgcacttccTTGGTCGTTCTGTGGGCCAAGCCGGATACCGCCGAGCGGGTTGTCGACTGGTTGGGTAAGAACTACCTGCAAAATCGGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CCCTA‑CAATCTGAAAGAGCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGTTTGGGCAAAAATCAGGGGATGTTTTTT‑CCGCACGAc                                                                                                                                                                              <  2:461903/242‑1 (MQ=255)                                                                 cccGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGc                                                                                                                                                                          <  1:119443/251‑1 (MQ=255)                                                                      ccTTGTTCGCTCTGTGTGACAAGCCGGATACCGCC‑AGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAACATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCCCGGGTACTGGAAAAATAAATGAAAACACT‑‑CAATCTTAAAGATCACAATCAGCCGGTCAG‑ATTTCCGCAAGCCGGAA‑CCCCGGGGGTTGGGCACAA‑TCAGGGGCTGGTTTTT‑TCGCACGACCTGCC‑ggcat                                                                                                                                                                   >  1:291432/1‑248 (MQ=255)                                                                              gCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTTGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCct                                                                                                                                                            <  1:596710/251‑1 (MQ=255)                                                                                      ggACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CGCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGaaa                                                                                                                                                    <  2:453054/250‑1 (MQ=255)                                                                                              cGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGTGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATgag                                                                                                                                            <  2:587062/251‑1 (MQ=255)                                                                                                                  gCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTGGCCGGCTGGATACGGCGGGCGAACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTg                                                                                                                         <  2:511198/250‑1 (MQ=255)                                                                                                                          gTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGCTACGGCGGGCGCACGAGTACTGGAAAAATAAAAGAAA‑ACCTAACAA‑CTGAAAGAACAACATAAGGAGGGCAGATTTTGGGCA‑GCCGTAAACCC‑AGGGGTGGGGAAAAAAT‑AGGGGGTGTTTTTTTCCGC‑CCACCTGCCCGGATTTC‑GCCTGACTGGCATTGTTGAGATggttaagcagggttttttcaaccccc                                                                                                                >  2:314670/1‑226 (MQ=255)                                                                                                                             gggTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGGTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAAGTGATGAGATGCGGAAGCTTGATTTTGTCACCCTCAGTg                                                                                                             <  2:415338/251‑1 (MQ=255)                                                                                                                               gTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAAGGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGAGGCTGAAGCTGGATTTTGTCACCCGCAGTGCg                                                                                                           <  2:94810/251‑1 (MQ=255)                                                                                                                                  aGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAg                                                                                                        <  1:408943/251‑1 (MQ=255)                                                                                                                                  aGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAAATAAA‑CTATA‑CAATCTGAAAGATCACAATGAGCAGGTCCG‑CTTTTCCCAAGACGTAA‑CCC‑CGGGGGTGGGCGAAAATCAGGGGCTGTTTTTT‑TCTCCCCACCTGCCCGGATTCCACCCTG‑CTGAAATTGATGAGATGCTTAATCTGGAttttttttccccccgtgcgg                                                                                                          >  1:297698/1‑233 (MQ=255)                                                                                                                                               aaaTCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGGGTACTGGGAAACTAAATGAAA‑ATCTA‑CAATCTGAAAAATTACAATGAGCAGGGCAG‑CTTTTCGCCAGCCGTAA‑ACC‑CGGGGGTGGGCCAAAATCAGGGGGTGTTTTTT‑CCGCACGACCTGCC‑GGGATTTAGCCTGCATGTAATTTATGAGATGCTTGAGCTGGGTTTTTTCACCCGCCGTGCGGCGATCCTCTCGGCGt                                                                                           >  1:67096/1‑251 (MQ=255)                                                                                                                                                       aaGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGAGCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGt                                                                                   <  2:114386/251‑1 (MQ=255)                                                                                                                                                                 tCATATTTGCCGGCTGGATACGGCGGGCGAACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGAGCACAAGGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATc                                                                          <  2:101763/250‑1 (MQ=255)                                                                                                                                                                       ttGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAg                                                                    <  2:485798/250‑1 (MQ=255)                                                                                                                                                                          ccGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGaaa                                                                <  1:494094/251‑1 (MQ=255)                                                                                                                                                                                ggATTCGGCGGGCGCACGAGTACTGGAAAACAAAATGAAA‑CTCTA‑CAATCGGAAAGATCACAAGGATCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGGTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCATGAAATCCtt                                                           <  2:489242/250‑2 (MQ=255)                                                                                                                                                                                      ggcgggcgCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCGGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAgc                                                    <  1:323595/251‑1 (MQ=255)                                                                                                                                                                                               aCGAGGACTGGAAAACTAAATGACA‑CTCTC‑CAATCTGAAAGATAACAATGAGCAGATCAG‑CATTGCGCAAGCCGTAA‑CCC‑AGGGCATGGCCAAAAATAAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCATGACTGAAATTGATGAGATGCTGCAGATGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGATTATTGGTGTTGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTgcgcg                                           >  2:372365/1‑251 (MQ=255)                                                                                                                                                                                                        ggAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTg                                  <  2:535405/251‑1 (MQ=255)                                                                                                                                                                                                                aaaTGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCg                          <  1:410505/251‑1 (MQ=255)                                                                                                                                                                                                                 aaTGAAA‑CTCTA‑CAATCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCA‑GGCATTCAGCCTGACTGACATTGATGAGATGCTGCAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTGGCCTTCCCgg                         >  1:355499/1‑251 (MQ=255)                                                                                                                                                                                                                           ctA‑CAATCTGAAAGATCACAAGGGGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACCGAAATTGATGAGATGCTGAAGCTGGATGTTGGCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGt                  <  2:297839/249‑1 (MQ=255)                                                                                                                                                                                                                                aaTCTGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCa             >  2:352818/1‑250 (MQ=255)                                                                                                                                                                                                                                    tGAAAGATCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGtt        >  1:160289/1‑251 (MQ=255)                                                                                                                                                                                                                                    ggAAAGATCACAATGAGCAGGTCAG‑CTGTGCGCAAGCCGTAT‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGAGCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGtt        <  2:196954/250‑1 (MQ=255)                                                                                                                                                                                                                                          aTCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGc  >  1:233944/1‑251 (MQ=255)                                                                                                                                                                                                                                          aTCACAATGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGc  >  2:129398/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                         TCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAA‑CTCTA‑CAATCTGAAAGATCACAACGAGCAGGTCAG‑CTTTGCGCAAGCCGTAA‑CCC‑AGGGGTTGGGCAAAAATCAGGGGCTGTTTTTT‑CCGCACGACCTGCC‑GGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGC  >  NC_000913/3525‑4006

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