Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,244,539	A→G	20.2%	N179N (AAT→AAC)	malF ←	maltose transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,244,539	0	A	G	20.2%	83.2 / 2.1	46	N179N (AAT→AAC)	malF	maltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (11/24);  new base G (7/3);  total (19/27)
Fisher's exact test for biased strand distribution p-value = 6.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGGGCTTAGCTTTTCATCACCCCAGTTGCCGTCGGCGGTAATGGACTGGTAAAAGCCAATTTGGTTATTCGGACGATATTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTC  >  NC_000913/4244315‑4244757                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  ccGGGGCTTAGCTTTTCATCACCCCAGTTGCCGTCGGCGGTAATGGACTGGTAAAAGCCAATTTGGTTATTCGGACGATATTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCGCTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGAtt                                                                                                                                                                                                                                 <  1:398536/227‑1 (MQ=255)          aGCTTTTCATCCCCCCAATTGCCGTCGGCGGTAATGGACTGGTAAAAGCCAATTTGGTTATTCGGACGATATTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATc                                                                                                                                                                                                                       <  1:243809/228‑1 (MQ=255)                                                   gTAAAAGCCAATTTGGTTATTCGGACGATATTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TTGCACTCCACGCCCTACGGTTTTTGGGAAACACCGCCCGG‑ATTCGCGCGTCGGCCTCGGGGTTgggggtccg                                                                                                                                                                             >  1:119630/1‑222 (MQ=255)                                                                                atTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGGTGCCCTCCGGCAGAATGGCGGGAA‑TGTCACACCGCGACCGTCGGGTTTGGGTAATTA‑CCTGCCGAATTCTCG‑GTCCGCCCTCCGGCTGGGCGGCGCTCTC‑TTTTCGTTTCA‑Gttctttttgctct                                                                                                                                              >  2:254145/1‑219 (MQ=255)                                                                                 tttCGCGCCGGTTTGGTTTTTCGGCAACGTGCCGTCTCCGTCGAGTGTGTAGAGCGGCTGGGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑gccgcc                                                                                                                                           <  1:325467/228‑1 (MQ=255)                                                                                         tcgtctGGTTATTGGTCAACGTGCCGTCACCGGCGGGTGTGTAGAGCGGTTGCGTGCCAGAAAACTGGGGGAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTaa                                                                                                                                   <  2:162197/225‑1 (MQ=255)                                                                                          ggtctGTTTTTTCGTCAACGTGCCGTCACCTTAGAGTGTGGAGAGCGGCTGCGTGCCAGAAAAATGGCGCCGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTaa                                                                                                                                   <  2:309617/225‑1 (MQ=255)                                                                                          ggtcggATTATGCGGCAACGTGCCGTCCCCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTaa                                                                                                                                   <  1:609507/222‑1 (MQ=255)                                                                                            tctGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTcgcgcg                                                                                                                                                                                                    <  1:347903/163‑1 (MQ=255)                                                                                                attattCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGGGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGGCGTGTC‑TTTTAGTTGCC‑GTTTTTGCTT‑GTCGCCCACTTTtaaacggcc                                                                                                                            >  1:471101/1‑223 (MQ=255)                                                                                                     tCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGGAA‑TGTCCCTCAGCGCCTCACGATTTTGGGGAAACC‑CGGGCGG‑ATTTGAGCGTTACCCCTCGGGCTGGGGGGGCGGTTC‑TTTCTGTTTGC‑GTTTTTTCTCTGCCCCCAACTTTAAAAggccggggg                                                                                                                        >  2:409552/1‑220 (MQ=255)                                                                                                     tCGTCAACGTGCCGTCACCGTCGAGTGGGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCCGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGCAA‑TGGGCATCAGCGGCTGCCGCTTCTTGGTAATCA‑CACCCAG‑GTTTGACCGTTCGCCTCTGGGCTGGGCGGGGGTTTC‑TTTCTGTGGTG‑GGTTTTGTGC‑GCCGCCCAATTTaaaataggcgggg                                                                                                                        >  2:419347/1‑218 (MQ=255)                                                                                                      cGTCAAGGGCCCGGCACCGGCGAGTGTGTAGAGCGGCTGCGGGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGaga‑g                                                                                                                     <  1:283412/228‑1 (MQ=255)                                                                                                          aaCGTGCCGTCACCGTCGAGTGTGTAGAGCGGCGGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACCCAAACCCTGACGATTCTCGGGAATCA‑CCCGCAG‑ATTCGCCCGGTCGCCTTCGGGCTGGGGGGGCGTTTC‑TTTCAGTTTCG‑GTGTTTTTTC‑TCCCCCCAATTTAAACACGCCGGGGA‑GGTa‑a                                                                                                                >  2:129894/1‑228 (MQ=255)                                                                                                            cGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTTGGCGGGCGTTTC‑TTTCCGTTGCA‑GGTTTTGCTC‑GCCCCCCAATTTAAAAGCGTCGGAGA‑GGGa‑act                                                                                                              >  1:169738/1‑226 (MQ=255)                                                                                                                     aCCGTCGAGTGTGTAGAGCGGCTGGGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGt                                                                                                     <  1:409552/228‑1 (MQ=255)                                                                                                                        gTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑ATTCCCTCAGCGCCTGACGATTCTGGGGAAACA‑CCCGCCG‑AGTCGCGCGGTCTCCTTCTGGCTGGGGGGGCGTTTC‑TTTCAGTTGCC‑GTTTTTTCTC‑GCCGCCCACTTTTAAAGCGTCGGGGAGGGTA‑TTTTTTGCCGGttt                                                                                                   >  2:546048/1‑228 (MQ=255)                                                                                                                         tCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGCAA‑AGGCACTCATCGCCTGACGGTTCTGCGTAATCA‑CGCGCCA‑CTTCGGGCGCTGGCCTTCGGGCCGGCGGGTCGTTTC‑TTTCAGTTGCAGGTTTTCGCCC‑CCGCCCACATTTAAAA‑CGCCGAGGG‑GGGA‑AATTTTATCGGTTTCg                                                                                                 >  2:36642/1‑228 (MQ=255)                                                                                                                              tgtgtAGAGCGGCTGCGTGCCAGAAAACTGGCCCAGGGAGCTCATCATCACTTTGTTGCCATCAGGCAGAATGGCGGGAA‑TGTCACCAAGCGCCCGCCCGTTCTTGGTAAAAA‑CCCCCAG‑GTTAGCCCCTTCGCCGTGCGCGGGGGCGGGCGTTCC‑TTTCAGTTGTA‑GTCTTTTTCC‑GCGGCCCAATATAATAGGGTGGGGGG‑GGGA‑ATTTTTGCCGGTTTTGCCggg                                                                                            >  2:495382/1‑226 (MQ=255)                                                                                                                                tgtAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAACGCCTGACGAGTCTTGGGAAACC‑CGCGCCG‑ATTCGCGCGGTGGCCCTCTGGGGGGGCGGTCGTTTCTTTTCAGTTCCG‑GTTTTTCCCC‑G‑CGCCACATTTTAAAACGTGCGGGG‑GGTG‑ATTTTTTGCGGTGTCGCCgctgc                                                                                          >  2:470992/1‑224 (MQ=255)                                                                                                                                      gCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCCCTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAgcg                                                                                    <  1:470992/228‑1 (MQ=255)                                                                                                                                         gCTGCGTGCCAGAAAACTGGGGCAGGGAGCTTATCATCCCTTTGGTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGcc                                                                                  <  2:119630/227‑1 (MQ=255)                                                                                                                                          cTGCGTGCCAGAAAACTGGGCCAGGGGGTCCATCATCACTTTGTTGCCATCCGGAAGAATGGCGGTAA‑TGTCCCTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAg                                                                                <  2:147579/228‑1 (MQ=255)                                                                                                                                          cTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGGTTCTGGGTAATCA‑CCGGCAG‑ATTCGCGCGTTTGCCTTCGGGCTGGGCGGGCGTTTC‑TTTCAGTTTCG‑GTTTTTGTTC‑GCCGCCCCATTTTAAAG‑GTCGCGGAGGGGA‑AATTTTTCCGGTTTTCCCCTCCCCGAgcgcacg                                                                                >  2:93457/1‑225 (MQ=255)                                                                                                                                          cTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGGCACTCAGCGACCGACGATTCTGGGGCATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCCGGCTGGGCGGTCGTTCT‑TTTCCGATGCA‑GTATTTCCTG‑GCCGCCCAAATTACAAAAGTCGGAGG‑TGGA‑CTTTTTGCAGCTTTCCCCGCCGCAGGGCACCa                                                                                 >  2:341364/1‑227 (MQ=255)                                                                                                                                          cTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑AGGCACTCAGCGCCTCACGGATCTGGGGAATTA‑CCCGCAG‑ATTCGCGCGTGCCCCTTTCGGCTGGGGGGGCGGTTT‑TTTTAGTTGCA‑TTGTTTGTCC‑GCGCCCACATTTGAAAGCGCTGGTGA‑GGGT‑AATTTTGTCCGTTTTGCCCTcccctaggcgcc                                                                                 >  2:363223/1‑216 (MQ=255)                                                                                                                                                            gcgcAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑AGTCCCTCAGCGCCCGACGATTCTGGGTAATCA‑CGCGCCG‑ATTCGCCCGTTTCCCTTCCGGCTGGGCGGCCGTTTA‑TTTCCGTTGCA‑GGTTTTTCTC‑TCCGCCCAATTTTAAAGCGGCGCGGG‑GGGA‑ATTTTTGCCGGTTTCGCCGTCTCTGAGCGCCAGTTGCCCCCCATcacccgc                                                              >  2:435027/1‑228 (MQ=255)                                                                                                                                                             ggcAGGGGGTTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCCCTCGGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCGTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCGC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATcgcccgc                                                              <  2:437065/226‑1 (MQ=255)                                                                                                                                                             cgcgGGGAGCTCATCATCTCTTTGTTGCCCTCCGGCAGAATGGCGGTAA‑TGTCGCTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATcgcccgc                                                              <  1:72352/227‑1 (MQ=255)                                                                                                                                                                 tggtgCTCATCATCAGTTTGGTGCCATCAGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGAGTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTGCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGt                                                         <  1:641745/224‑1 (MQ=255)                                                                                                                                                                        catcatccCTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCCCTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGAcc                                                  <  1:592803/228‑1 (MQ=255)                                                                                                                                                                         atcctcaCTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGGTTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGAcc                                                  <  1:234762/227‑1 (MQ=255)                                                                                                                                                                         atcatccCTTTGTTTCCATTCGGCAGAATGGCGGTAA‑TGGCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGAcc                                                  <  1:163197/227‑1 (MQ=255)                                                                                                                                                                          tcattaCTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCa                                                 <  1:546048/227‑1 (MQ=255)                                                                                                                                                                          tcatcaCTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGGCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCa                                                 <  1:262759/227‑1 (MQ=255)                                                                                                                                                                                cTTTGTTGCCCTCCGGCAGAATGGCGGTAATTGTC‑CTCAGCGCCTGGCGGTTCTGGGTAATCT‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTta                                           <  1:522038/227‑1 (MQ=255)                                                                                                                                                                                 tttGGTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTtat                                          <  1:5957/227‑1 (MQ=255)                                                                                                                                                                                   tgttgCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCCG‑ATTCGCGCGGTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGAA‑GTTTTTGCTC‑TCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCCGTTTCCCCTCATCGCCCGCCGGGGAAAAACCAGAGTtatacg                                       >  1:382842/1‑226 (MQ=255)                                                                                                                                                                                   tgttgCCATCCGGCAGAATGGCGGTAA‑TATCACTCAGCGCCTGACGGTTCTGGGGAATCA‑CCCGCAG‑GTTCGCGCGTTCGCCCTTCGGCGTGGCGGGCGTTTT‑TTTTTGTTGCA‑GTTTTTGGTC‑GCCGCCCCATTTAAAAACGGCAGTGG‑AGTAGATTTTTTCGGTTTTCCCATTCCCTTGACGCCATTTGCACCTACTCGCCCCCCGGGTACAGCCAAAAGtaaa                                          >  2:638550/1‑223 (MQ=255)                                                                                                                                                                                    gttgCCATCCGGCAGAATGGCGGTAA‑TGTCCCTCAGCGCCTGGCGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTa                                                        <  1:531706/210‑1 (MQ=255)                                                                                                                                                                                            ccGGGAGAATGGCGGTAA‑TTTCACTCAGCGCCTGACGGTTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTg                              <  1:476264/228‑1 (MQ=255)                                                                                                                                                                                                           tAA‑TGTCACTCAGCCCCTGACGATTCGGGGTAATCA‑CGCGCAG‑ATTCGTGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACGCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGa                <  1:84658/227‑1 (MQ=255)                                                                                                                                                                                                                gTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTa            <  1:526079/227‑1 (MQ=255)                                                                                                                                                                                                                   acccagccacTGACGATGCTGGGTAATCG‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTaac          <  1:254145/217‑1 (MQ=255)                                                                                                                                                                                                                         cgcCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCGGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATTGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCCGCTTGCCCGGAGAGCTCCAACAACACTTc  >  2:271107/1‑228 (MQ=255)                                                                                                                                                                                                                         cgcCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATCGGTTTTGCCTGCTTGCCAGGAGCGGTCTAACAACACTTc  >  1:46829/1‑228 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  CCGGGGCTTAGCTTTTCATCACCCCAGTTGCCGTCGGCGGTAATGGACTGGTAAAAGCCAATTTGGTTATTCGGACGATATTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAATCA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGTA‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTC  >  NC_000913/4244315‑4244757

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