Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	2,152,038	G→C	16.7%	A362G (GCG→GGG)	yegI ←	protein kinase YegI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	2,152,038	0	G	C	16.7%	124.5 / 2.6	66	A362G (GCG→GGG)	yegI	protein kinase YegI
Reads supporting (aligned to +/- strand):  ref base G (34/19);  new base C (3/8);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 4.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTTGTACGCTTTGCTTGTCAGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGC  >  NC_007779/2151804‑2152272                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 gCTTTGTACGCTTTGCTTGTCAGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGGTGCAATGCTGGCTGCG‑CTACTGGCGGCa                                                                                                                                                                                                                                               >  1:1764329/1‑238 (MQ=255)      gTACGCTTTGCTTGTCAGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCGGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑a                                                                                                                                                                                                                                        <  1:2043331/239‑1 (MQ=255)         cGCTTTGCTTGTCAGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGAAGCCATCACC‑aaa                                                                                                                                                                                                                                      >  1:1706235/1‑236 (MQ=255)            tttGCTTGTCAGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGccc                                                                                                                                                                                                                                  <  1:465229/239‑1 (MQ=255)            tttGCTTGTCAGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGG‑TTACTGACGCCATCACC‑ATCGcc                                                                                                                                                                                                                                   >  1:45031/1‑238 (MQ=255)                cTTGTCAGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACAAAAATAATATAACGCGGTTCTGCCTGAAGGCAGAAAGGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGCCCAG‑a                                                                                                                                                                                                                              <  2:1705656/238‑1 (MQ=255)                     cAGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑GTAATGACGCCATCACC‑CTCCCCCAG‑ACTTTc                                                                                                                                                                                                                         >  1:438270/1‑238 (MQ=255)                        aCTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑GTACCTACGCCCTCAAC‑AACACCCAG‑ACTTTCACCa                                                                                                                                                                                                                     >  1:1679511/1‑239 (MQ=255)                         cTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGc                                                                                                                                                                                                                       <  1:1709752/236‑1 (MQ=255)                               ataatCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑a‑aa                                                                                                                                                                                                            <  1:1503090/239‑1 (MQ=255)                                  atCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑GCACCGGACCCATCACC‑AAAGCCCCA‑ACATTCCCCAGGAC‑a‑caac                                                                                                                                                                                                          >  1:1563244/1‑237 (MQ=255)                                  acccAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑a                                                                                                                                                                                                                                        <  1:800186/209‑1 (MQ=255)                                       gATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGATCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑GCACCGGCGCCATCCCC‑CACGCCCAG‑ACCTTACCCCACCC‑A‑AAA‑ATcccc                                                                                                                                                                                                    >  1:1846269/1‑237 (MQ=255)                                             cgccAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑GTACTTAAGCCCTCACC‑AACGACCAC‑ACTTTCGCCCGCAC‑ACAAA‑TCCccgcgcggg                                                                                                                                                                                              >  2:1743051/1‑234 (MQ=255)                                             cgccAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑GTACTGACGCCATCACC‑AACGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCggtgg                                                                                                                                                                                             >  1:1631616/1‑239 (MQ=255)                                                                     tataCGCGGGTGTGCCCGGAGGCCACAAAAAAGGAGGCAATAAGGCTGA‑AGTGGAATCCCAAGGAGAAATATTTTTTGCGCGCGGAAAAAGCCATAAAGGAAGGGCCC‑TGCCAGTCGGTTGGAAATGATCGGTTAATGCCAGTTGCAATGCGGGCGGTG‑CACCTGACGCCATCCCC‑ATCCCCCAG‑ACTTTCGCAACCAA‑A‑AAA‑ATCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCa                                                                                                                                                                      <  2:675454/236‑1 (MQ=255)                                                                       accgcgGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGAGGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATa                                                                                                                                                                   <  1:1135002/237‑1 (MQ=255)                                                                       aaCGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGG‑CTACTGACGCCCTCACC‑ATCGCCCAG‑ACCTTCGCCCGCAC‑A‑AAA‑TTCCCGCCGGTGGTCATGACCTCTTCGCCGCGCTCAATa                                                                                                                                                                   >  1:1163583/1‑239 (MQ=255)                                                                         cgcgGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAGGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CCACTGACCCCCTCAAC‑AACGCCCAC‑ACTTTCCCCCGCCC‑A‑AAATTTCCCCCCGGCGG‑AATGGACCCTTCGCCgcggtccaataa                                                                                                                                                                 >  2:142631/1‑228 (MQ=255)                                                                         cgcgGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGG‑GTTCTGACGCCCTCACCAACCCCCCAA‑ATTTTCCCCA‑CACAA‑AAA‑TTCCCCCCGGGGGGAATAACCTCTTCCCCG‑GATCaaaaaa                                                                                                                                                                 >  2:1671510/1‑239 (MQ=255)                                                                         cgcgGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGG‑CTACTGACGCCATCACC‑AACGCCCAG‑ACTTTTGCCAGCCC‑A‑AAA‑TTCCCCCCCGGGGGAATGACCTCTTCCCCGAGATCCATaaa                                                                                                                                                                 >  1:616542/1‑239 (MQ=255)                                                                              ttCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGAG‑CTACTGAAACCCTCACC‑ATCGCCCCG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGCGGGAATGACCTCTTCCCCGAGCTCAATAAAATaca                                                                                                                                                            >  1:98423/1‑237 (MQ=255)                                                                                ttgccgagaGGCCGAAAAGGCGGA‑GCAATAAAGAGGACAGGGGGACCTCAAGACAAAAGGTTTATTCCGCCCGGAAACGGGCAAAAGGGAGGCGCCC‑TGCCGGTCGGTGGAAAATGATCGGTAAATGGCAGTGCCAATGCTGGGGGGG‑CTACTGCCCCCATCCCC‑ATCCCCAAG‑CTTTCCGCAAGCAA‑A‑AAA‑TTCCCGCCGGTAGTAATGACCTCTTCGCCGAGATCAATAAAATAGAt                                                                                                                                                           <  2:616542/231‑1 (MQ=255)                                                                                             gaaaaGCAGAA‑GCGATAATGTTGAAAGTGGGATCTCAAGAAGAATGTTTTTTT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ACCGCCCAG‑CCTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTAt                                                                                                                                             <  1:1323578/237‑1 (MQ=255)                                                                                                               aTGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGGGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATc                                                                                                                                            >  1:40108/1‑223 (MQ=255)                                                                                                                tGCTAACGGTGCGATCTAAAGAAGAATGAGTTATT‑CGCGCCGTACCAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGCCGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑TAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGc                                                                                                                          <  1:1671510/239‑1 (MQ=255)                                                                                                                    gACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAaca                                                                                                                                                       >  1:1407216/1‑207 (MQ=255)                                                                                                                     aCAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGGGGTAATGACCTCTTCGCCGCGATCAATAAAATAGATAACACCTTGATTATc                                                                                                                                            >  1:114699/1‑217 (MQ=255)                                                                                                                         tGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCCTCACC‑ATCGCCCCG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGGGGTAATGACCTCATCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACAACCGGCAgtcagtca                                                                                                                 >  1:497226/1‑239 (MQ=255)                                                                                                                           cGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGG‑CCACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑ATCCCGCCGGGGGTAATGACCTCTTCCCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCa‑cac                                                                                                                                 >  1:313593/1‑221 (MQ=255)                                                                                                                              gctcAAGAAGAATGATTTATT‑CGCGCCGTAACGGGCCAAAGGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAAGGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGCCCAT‑ACTTTCGCCAGCAA‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCgag                                                                                                                                                                          <  1:125308/177‑1 (MQ=255)                                                                                                                               ctcAAGCAGAATGTTGTATG‑CGGGCCGTAACGGGCATAAAGGAGGCGGCC‑GGCCAGTCGGTTGGAAATGATGGGGAAAGGGCAGTTGCAATGCTGGCGGCG‑CCTCGGACGCCATCCCC‑ACCGCCCAG‑ACTTTCGCCAGAAA‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAg                                                                                                           <  2:1788457/239‑1 (MQ=255)                                                                                                                               ctcAAGCAGAATGATGTTTT‑CGCGCCGTAACAGGCTAAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACGGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAg                                                                                                           <  1:1609831/239‑1 (MQ=255)                                                                                                                               acaggcacaaaTGATGGTTG‑CGCGGCGTAACGGGCCAGAAGGAGGGGGCC‑TGCCGGGCGGTTGAAAATGACCGGGTAATGGTAGTTGCAAGGCGGGCGGTG‑CTCAGGACCCCACCACC‑ATCGCCCGT‑ATTTCCGCCACAAA‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAg                                                                                                           <  2:45031/230‑1 (MQ=255)                                                                                                                                   aGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CGACTGACGCCATCACC‑CACGCCCCG‑ACTTTCGCCCGCCC‑A‑AAA‑TTCCCCCCCGGGGTAATGACCCCCTCCCCGCGAACAAAAAAAAAGAAAACCCCCTGATTATACCACCCA‑CCCCACCGGCAAACAGGCAAAAGAGCggg                                                                                                       >  2:1944558/1‑239 (MQ=255)                                                                                                                                   aGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGG‑CTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCCGCAC‑A‑AAA‑GTTCCCCCGGTGGGAAt                                                                                                                                                                                         >  1:2127868/1‑158 (MQ=255)                                                                                                                                        aaTGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGG‑CTTCCGACCCCCTCACC‑CTTCCCCCA‑ACCTTTGCCCGCAC‑A‑AAA‑ATCCCCCCCGGGGGAAATTCCTCTTTCCCCGGATCCAAAAAAAAAATTACACCCTTGTTATCCAACCCA‑CACCCCAGGCCAACAGTCAAAATAACGGGGGaaa                                                                                                  >  2:2149833/1‑239 (MQ=255)                                                                                                                                        aaTGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGG‑CCACCGAAGCCCTCCCC‑CCCCCCCAG‑ACCTTCCCCCGCCC‑A‑AAA‑ATCCCGCCGGGGGTAAAGAACTCCTCCCCCAGAACAATAAAAAAAAAAACACCCTGGTTTACCACCCCACCACCAC‑GGCAGTCAAGCAAAAGAGCGGGGGaaa                                                                                                  >  2:2001766/1‑239 (MQ=255)                                                                                                                                            atgTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGc                                                                                                                                                                                                                                    <  1:846470/110‑1 (MQ=255)                                                                                                                                               atatcccgccggTAACGGGCC‑AAAGGAAGGGGCCTTGCCGGTGGGTGGGAAAGGATCGGGTAAGGGCGGTTGCAATGCGGGCGGTC‑CTACCCAGCCCATCCCC‑ACCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCa                                                                                            <  2:98423/228‑1 (MQ=255)                                                                                                                                                  t‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGAAAACACCTTGATTATCCAGCGCA‑CACCCCGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCAtt                                                                                          >  1:1651264/1‑237 (MQ=255)                                                                                                                                                  t‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGGGGGAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGGAAACATGCAtt                                                                                          >  1:86694/1‑237 (MQ=255)                                                                                                                                                          gTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTT‑CTAACGACGCCATCCCC‑ACCGCCCCG‑ACTTTCGCCCGCAC‑A‑AAA‑TTCCCCCCGGCGGGGATGACCCCCTTGCCCCGGTCAAAAAAAAAAATAACAACTTGGTTTTTCTGCCCC‑CACCCCAGGCCGTCAGTCAAAATAGCCGGGTAAACCAGCAATTCTGAAACg                                                                                 >  2:1782770/1‑239 (MQ=255)                                                                                                                                                          gTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTC‑CTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTTCCGCCGGGGGTAATGACCCCTTCGCCGGGACCCATAAAATa                                                                                                                                                              >  1:832783/1‑163 (MQ=255)                                                                                                                                                                 gcagaaggggaGCGGCC‑TGCCGGTCGGTGGGAAATGACCGGGTAATGGCAGTTGAAATGCGGGCGGTG‑CTCCTGCCCCCATCACC‑ACCCCCCAG‑ACTTCCGCAAGAAA‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACAt                                                                           <  2:86694/228‑1 (MQ=255)                                                                                                                                                                           aGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGGCCTCCGGAGGCC‑TCACC‑AACGCCCCGAACTTTCCCCGCCCC‑A‑AAT‑TTCCCCCGGGGGGTAA‑GACCCCTTCACCGAGACCAAAAAAAAAAATAACACCCTTATTATCCAGCGCC‑AACCCCGGGCAGGCAGACAAATAAACGGGGGAAAAATGTCATGCTCGAACCATACAATTCtttt                                                                    >  2:1619251/1‑234 (MQ=255)                                                                                                                                                                             ggccc‑TCCCAGTGGTTTGAAATTGATCGGGTAAGGGCAGTTGCAAGGCGGGCGGTG‑CTTTTGACGCCACCACC‑ACCCCCCAG‑ACTTCCGCCAGAAA‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTg                                                                <  1:2001766/234‑1 (MQ=255)                                                                                                                                                                                cc‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGGG‑CTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTTGCCAGCAC‑A‑AAA‑TTTCCGCCGGGGGGAATGACCTCCTCGCCGGGATCAATAAAATAGATAACACCCTGATTATCCACCGCC‑CACCCCGGGCAGACAGTCAAAAGAGCGGGGGAAACATGGATTGCTGAAACGATACATTTCttttaatttttgc                                                           >  2:1309119/1‑229 (MQ=255)                                                                                                                                                                                                gAAATGATCGGTTAATGGCTGTTGCAATGCTGGCGGTG‑CTCCGGCCGCCATCCCC‑ATCGCCCAG‑ACTTTCGCAAGCAA‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATaaa                                                                                                                                                                 <  1:1535143/123‑1 (MQ=255)                                                                                                                                                                                                       tCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATTGCCCCG‑ACTTTTGCCCGCAC‑A‑AAA‑ATTCCCCCGGGGGGAAAGACCTCCTCGCCCAGAACAATAAAAAAGATAACACCCTGATTATTCAGCCCC‑CACCACCGGGAATCAGACAAAAGAGCGGGGGAAAACAGCCTTGCCTAAACGATACATTTCTTTAACAGTTGGGGTAGAGAGTCCAGTGCCGCTa                                      >  2:1855237/1‑238 (MQ=255)                                                                                                                                                                                                       tCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCCGCAC‑A‑AAA‑TTCCCCCCGGTGGTAATGACCTCTTCGCCGCGATCCATAAAATAGATAACACCTTGATTATCCAGCGCC‑CACCACGGGCAGTTAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTa                                      >  1:26480/1‑238 (MQ=255)                                                                                                                                                                                                       tCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACACCCTCCCC‑ATCGCCCCG‑ACTTTCGCCACCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCCGGATCCATAAAATAACTAACCCCTTGCTTATc                                                                                                                                            >  1:1428502/1‑137 (MQ=255)                                                                                                                                                                                                       tCGGGTAATGGCAGTTGCAATGCTGGCGGGG‑GTACTTACGCCCACACC‑AACGCCCAG‑ACTTTCGCCCGCAC‑A‑AAA‑TTTCCCCCGGGGGGAAATGCCTCTTCCCCCGGGTCCAAAAAAAAGAGAAACCCTTTAATATCCCACCCA‑CCCCCCGGGCAGTCAGGCCAATGAACGGGGGAAAAATGCCATGGTGAAAAGAAACATTTTTTTAACCGTTGGGGTAAAGAGGCCAGGGcccccc                                      >  2:1432511/1‑234 (MQ=255)                                                                                                                                                                                                           gTAATGCCAGTTGCAATGCTGGCGGTG‑CTACTGACCCCACCCCC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCa                                  <  1:1309119/238‑1 (MQ=255)                                                                                                                                                                                                           gTAAGGGCGGTTGAAAAGCTGGCGGTG‑CTCCGGCCCCCACCACC‑ACCGCCCAG‑ACTTCCGCAAGAAA‑A‑AAA‑TCCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCa                                  <  2:55688/238‑1 (MQ=255)                                                                                                                                                                                                                 gCAGTTGAAATCGTGGCGGTG‑CTCCTGACGCCATCACC‑ATCGCCCAG‑ACTTCCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCttt                            <  2:1966564/238‑1 (MQ=255)                                                                                                                                                                                                                    gTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGAAAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGACAAATGAGCCGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTaac                         >  1:1966564/1‑236 (MQ=255)                                                                                                                                                                                                                                gCGGTG‑CTACTGACGCCATCACC‑AACGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTTCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAAAAAAATAGATAACACCTTGATTATCCAGGGCA‑CACCACGGGCAATCAGTCAAAAGAGCGGGGTAAACATGCAATGCTTAAACGATACATTTCTTTAACTGGTGGCGAAGAGAGTCCCGGGCCGCCACCCCCGCTTTAGCCGTTGGTCGCCCg            >  2:1785579/1‑239 (MQ=255)                                                                                                                                                                                                                                       cTACTGACGCCCTCACC‑AACGCCCAA‑ACTTTTGCCCGCAC‑A‑AAA‑TTTCCGCCGGGGGTAATGACCTCTTCGCCGAGATCAATAAAAAAGATAAAACCTTGATTATCCAGCGCA‑CACCACGGGCAGGCAGTCAAAACAACGGGGGAAACATGCAATGCTGAAACGATACATTTTTTTAACTGTTTGCGGAAAGAGTCCAGTGCCGCCACCCACGCTTTAGACCCCGGACGCCCGGTTgcc      >  2:1194796/1‑239 (MQ=255)                                                                                                                                                                                                                                       cTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCTGGTTgcc      >  1:643389/1‑239 (MQ=255)                                                                                                                                                                                                                                        tACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCCCTTCGCCGAGATCACTCAAATAGATAACACCTTGATTATCCAGCCCC‑CACCCCGGGCCGTCAGTCAAATGAACGGGGTAAACATGCATTGCTGACACGCTCCATTTCTTTAACCGTTGGCGTAGAGAGTCCAGTGCCCCTACCCACGCCTTAGCCGTCGGACGCCCgggtccc      >  2:1241308/1‑233 (MQ=255)                                                                                                                                                                                                                                         ccTGACCCCACCCCC‑ATCGCCCAA‑ACTTCCGCCAACAA‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGTGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTgccac    <  1:1475219/238‑1 (MQ=255)                                                                                                                                                                                                                                           tGACGCCATCACC‑ATCGCCCAG‑ATTTTCGCCAGCAA‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTgccacgc  <  1:142631/239‑1 (MQ=255)                                                                                                                                                                                                                                           tGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGaaa                                                                                   >  2:393357/1‑158 (MQ=255)                                                                                                                                                                                                                                           tGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAA‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGaaa                                                                                   <  1:393357/158‑1 (MQ=255)                                                                                                                                                                                                                                           tGACCCCATCACC‑ACCCCCCAA‑ACTTCCCACAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTgccacgc  <  1:1582833/239‑1 (MQ=255)                                                                                                                                                                                                                                               gCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGcagtcagt                                                                                                                   >  2:996260/1‑122 (MQ=255)                                                                                                                                                                                                                                               gCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGcagtcagt                                                                                                                   <  1:996260/122‑1 (MQ=255)                                                                                                                                                                                                                                               gCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAAc                                                                  >  2:1471746/1‑171 (MQ=255)                                                                                                                                                                                                                                               gCCATCACC‑ATCGCCCAG‑ACTTTCGCAAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAAc                                                                  <  1:1471746/171‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 GCTTTGTACGCTTTGCTTGTCAGACTGCCAATAATCCAGATAGCCGCCAGCACAGGAACCAAAATAATATAACGCGGTTCTGCCTGAAGGCCGCAAAGCAGCA‑GCGATAATGCTGACAGTGCGATCTCAAGCAGAATGATGTATT‑CGCGCCGTAACAGGCCTAAAGGAAGCGGCC‑TGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTG‑CTACTGACGCCATCACC‑ATCGCCCAG‑ACTTTCGCCAGCAC‑A‑AAA‑TTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCA‑CACCACGGGCAGTCAGTCAAATGAGCGGGGTAAACATGCATTGCTGAAACGATACATTTCTTTAACTGTTGGCGTAGAGAGTCCAGTGCCGCTACCCACGCTTTAGCCGTCGGACGCCCGGTTGCCACGC  >  NC_007779/2151804‑2152272

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