Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,244,538	A→G	22.1%	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,538	0	A	G	22.1%	84.0 / 3.8	51	N179N (AAT→AAC)	malF	maltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (13/25);  new base G (9/3);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGGGCTTAGCTTTTCATCACCCCAGTTGCCGTCGGCGGTAATGGACTGGTAAAAGCCAATTTGGTTATTCGGACGATATTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTC  >  NC_000913/4244314‑4244756                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                    ccGGGGCTTAGCTTTTCATCACCCCAGTTGCCGTCGGCGGTAATGGACTGGTAAAAGCCAATTTGGTTATTCGGACGATATTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCGCTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGAtt                                                                                                                                                                                                                                   <  2:398536/227‑1 (MQ=255)          aGCTTTTCATCCCCCCAATTGCCGTCGGCGGTAATGGACTGGTAAAAGCCAATTTGGTTATTCGGACGATATTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑c                                                                                                                                                                                                                        <  2:243809/228‑1 (MQ=255)                                                   gTAAAAGCCAATTTGGTTATTCGGACGATATTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TTGCACTCCACGCCCTACGGTTTTTGGGAAA‑CACCGCCCGG‑ATTCGCGCGTCGGCCTCGGGGTTgggggtccg                                                                                                                                                                              >  2:119630/1‑222 (MQ=255)                                                                                atTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGGTGCCCTCCGGCAGAATGGCGGGAA‑TGTCACACCGCGACCGTCGGGTTTGGGTAAT‑TA‑CCTGCCGAATTCTCG‑GTCCGCCCTCCGGCTGGGCGGCGCTCTC‑TTTTCGTTTCA‑Gttctttttgctct                                                                                                                                               >  1:254145/1‑219 (MQ=255)                                                                                 tttCGCGCCGGTTTGGTTTTTCGGCAACGTGCCGTCTCCGTCGAGTGTGTAGAGCGGCTGGGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑gccgcc                                                                                                                                            <  2:325467/228‑1 (MQ=255)                                                                                         tcgtctGGTTATTGGTCAACGTGCCGTCACCGGCGGGTGTGTAGAGCGGTTGCGTGCCAGAAAACTGGGGGAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTaa                                                                                                                                    <  1:162197/225‑1 (MQ=255)                                                                                          ggtctGTTTTTTCGTCAACGTGCCGTCACCTTAGAGTGTGGAGAGCGGCTGCGTGCCAGAAAAATGGCGCCGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTaa                                                                                                                                    <  1:309617/225‑1 (MQ=255)                                                                                          ggtcggATTATGCGGCAACGTGCCGTCCCCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTaa                                                                                                                                    <  2:609507/222‑1 (MQ=255)                                                                                            tctGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTcgcgcg                                                                                                                                                                                                     <  2:347903/163‑1 (MQ=255)                                                                                                attattCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGGGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGGCGTGTC‑TTTTAGTTGCC‑GTTTTTGCTT‑GTCGCCCACTTTtaaacggcc                                                                                                                             >  2:471101/1‑223 (MQ=255)                                                                                                     tCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGGAA‑TGTCCCTCAGCGCCTCACGATTTTGGGGAAA‑CC‑CGGGCGG‑ATTTGAGCGTTACCCCTCGGGCTGGGGGGGCGGTTC‑TTTCTGTTTGC‑GTTTTTTCTCTGCCCCCAACTTTAAAAggccggggg                                                                                                                         >  1:409552/1‑220 (MQ=255)                                                                                                     tCGTCAACGTGCCGTCACCGTCGAGTGGGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCCGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGCAA‑TGGGCATCAGCGGCTGCCGCTTCTTGGTAAT‑CA‑CACCCAG‑GTTTGACCGTTCGCCTCTGGGCTGGGCGGGGGTTTC‑TTTCTGTGGTG‑GGTTTTGTGC‑GCCGCCCAATTTaaaataggcgggg                                                                                                                         >  1:419347/1‑218 (MQ=255)                                                                                                      cGTCAAGGGCCCGGCACCGGCGAGTGTGTAGAGCGGCTGCGGGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGaga‑g                                                                                                                      <  2:283412/228‑1 (MQ=255)                                                                                                          aaCGTGCCGTCACCGTCGAGTGTGTAGAGCGGCGGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACCCAAACCCTGACGATTCTCGGGAAT‑CA‑CCCGCAG‑ATTCGCCCGGTCGCCTTCGGGCTGGGGGGGCGTTTC‑TTTCAGTTTCG‑GTGTTTTTTC‑TCCCCCCAATTTAAACACGCCGGGGA‑GGT‑a‑a                                                                                                                >  1:129894/1‑228 (MQ=255)                                                                                                            cGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTTGGCGGGCGTTTC‑TTTCCGTTGCA‑GGTTTTGCTC‑GCCCCCCAATTTAAAAGCGTCGGAGA‑GGG‑a‑act                                                                                                              >  2:169738/1‑226 (MQ=255)                                                                                                                     aCCGTCGAGTGTGTAGAGCGGCTGGGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGt                                                                                                     <  2:409552/228‑1 (MQ=255)                                                                                                                        gTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑ATTCCCTCAGCGCCTGACGATTCTGGGGAAA‑CA‑CCCGCCG‑AGTCGCGCGGTCTCCTTCTGGCTGGGGGGGCGTTTC‑TTTCAGTTGCC‑GTTTTTTCTC‑GCCGCCCACTTTTAAAGCGTCGGGGAGGGT‑A‑TTTTTTGCCGGttt                                                                                                   >  1:546048/1‑228 (MQ=255)                                                                                                                         tCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGCAA‑AGGCACTCATCGCCTGACGGTTCTGCGTAAT‑CA‑CGCGCCA‑CTTCGGGCGCTGGCCTTCGGGCCGGCGGGTCGTTTC‑TTTCAGTTGCAGGTTTTCGCCC‑CCGCCCACATTTAAAA‑CGCCGAGGG‑GGG‑A‑AATTTTATCGGTTTCg                                                                                                 >  1:36642/1‑228 (MQ=255)                                                                                                                              tgtgtAGAGCGGCTGCGTGCCAGAAAACTGGCCCAGGGAGCTCATCATCACTTTGTTGCCATCAGGCAGAATGGCGGGAA‑TGTCACCAAGCGCCCGCCCGTTCTTGGTAAA‑AA‑CCCCCAG‑GTTAGCCCCTTCGCCGTGCGCGGGGGCGGGCGTTCC‑TTTCAGTTGTA‑GTCTTTTTCC‑GCGGCCCAATATAATAGGGTGGGGGG‑GGG‑A‑ATTTTTGCCGGTTTTGCCggg                                                                                            >  1:495382/1‑226 (MQ=255)                                                                                                                                tgtAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAACGCCTGACGAGTCTTGGGAAA‑CC‑CGCGCCG‑ATTCGCGCGGTGGCCCTCTGGGGGGGCGGTCGTTTCTTTTCAGTTCCG‑GTTTTTCCCC‑G‑CGCCACATTTTAAAACGTGCGGGG‑GGT‑G‑ATTTTTTGCGGTGTCGCCgctgc                                                                                          >  1:470992/1‑224 (MQ=255)                                                                                                                                tgtAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGGCACTCCGCGCCCTGCGAGTCTGCGTCATCCA‑GCCGCGC‑TTTCGC‑CGTTCGGCTTCTGGGGTGGCGGGCTTTTC‑TTTCAGTTGCA‑GGTTTTTTCC‑TCCGCCAAATTTTAAAACGTCGGGGA‑GGG‑G‑Attttttcgtggtttcccgtccc                                                                                          >  1:473758/1‑211 (MQ=255)                                                                                                                                      gCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCCCTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAgcg                                                                                    <  2:470992/228‑1 (MQ=255)                                                                                                                                       cGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGGCACTCAGCGCCTTACGGTTTTGGGTAAA‑CA‑CGCCCAG‑ATTTGCGCGTTCGCCTTCGGGCTGGGCGGGCCGTTC‑TTTCCGTTGCA‑CTTTTTGTGC‑GCCGCCCAATTTTAAACAGTCGGGGG‑GGTGA‑ATTTTTTCTGGttcgccctcccgtcgcg                                                                                     >  1:650148/1‑212 (MQ=255)                                                                                                                                         gCTGCGTGCCAGAAAACTGGGGCAGGGAGCTTATCATCCCTTTGGTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGcc                                                                                  <  1:119630/227‑1 (MQ=255)                                                                                                                                          cTGCGTGCCAGAAAACTGGGCCAGGGGGTCCATCATCACTTTGTTGCCATCCGGAAGAATGGCGGTAA‑TGTCCCTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAg                                                                                <  1:147579/228‑1 (MQ=255)                                                                                                                                          cTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGGTTCTGGGTAAT‑CA‑CCGGCAG‑ATTCGCGCGTTTGCCTTCGGGCTGGGCGGGCGTTTC‑TTTCAGTTTCG‑GTTTTTGTTC‑GCCGCCCCATTTTAAAG‑GTCGCGGAGGGG‑A‑AATTTTTCCGGTTTTCCCCTCCCCGAgcgcacg                                                                                >  1:93457/1‑225 (MQ=255)                                                                                                                                          cTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGGCACTCAGCGACCGACGATTCTGGGGCAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCCGGCTGGGCGGTCGTTCT‑TTTCCGATGCA‑GTATTTCCTG‑GCCGCCCAAATTACAAAAGTCGGAGG‑TGG‑A‑CTTTTTGCAGCTTTCCCCGCCGCAGGGCACCa                                                                                 >  1:341364/1‑227 (MQ=255)                                                                                                                                          cTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑AGGCACTCAGCGCCTCACGGATCTGGGGAAT‑TA‑CCCGCAG‑ATTCGCGCGTGCCCCTTTCGGCTGGGGGGGCGGTTT‑TTTTAGTTGCA‑TTGTTTGTCC‑GCGCCCACATTTGAAAGCGCTGGTGA‑GGG‑T‑AATTTTGTCCGTTTTGCCCTcccctaggcgcc                                                                                 >  1:363223/1‑216 (MQ=255)                                                                                                                                                    gAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGGAA‑TGTCACTCAGCGCCTGACGGTTCTGGGTAAT‑CA‑CGCGCGG‑ATTCGCGCGTTCGCCCTCGGGCCGGGCGGGCGTTTC‑TTTTCGTTGCT‑GTTTTTGTTC‑GCCCCCCAATTTTAAAACGTCGGTGA‑GGG‑T‑ATTTTTTCCGGTTTTCCCGTCCCttgagcccactgcccccacc                                                                      >  1:413110/1‑209 (MQ=255)                                                                                                                                                            gcgcAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑AGTCCCTCAGCGCCCGACGATTCTGGGTAAT‑CA‑CGCGCCG‑ATTCGCCCGTTTCCCTTCCGGCTGGGCGGCCGTTTA‑TTTCCGTTGCA‑GGTTTTTCTC‑TCCGCCCAATTTTAAAGCGGCGCGGG‑GGG‑A‑ATTTTTGCCGGTTTCGCCGTCTCTGAGCGCCAGTTGCCCCCCATcacccgc                                                              >  1:435027/1‑228 (MQ=255)                                                                                                                                                             ggcAGGGGGTTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCCCTCGGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCGTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCGC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATcgcccgc                                                              <  1:437065/226‑1 (MQ=255)                                                                                                                                                             cgcgGGGAGCTCATCATCTCTTTGTTGCCCTCCGGCAGAATGGCGGTAA‑TGTCGCTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATcgcccgc                                                              <  2:72352/227‑1 (MQ=255)                                                                                                                                                                 tggtgCTCATCATCAGTTTGGTGCCATCAGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGAGTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTGCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGt                                                         <  2:641745/224‑1 (MQ=255)                                                                                                                                                                        catcatccCTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCCCTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGAcc                                                  <  2:592803/228‑1 (MQ=255)                                                                                                                                                                         atcctcaCTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGGTTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGAcc                                                  <  2:234762/227‑1 (MQ=255)                                                                                                                                                                         atcatccCTTTGTTTCCATTCGGCAGAATGGCGGTAA‑TGGCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGAcc                                                  <  2:163197/227‑1 (MQ=255)                                                                                                                                                                          ttagctgtttgtgccaacccGGCAGAATGGGGGAAA‑TTGCCATCAGCGCCTGCCGATTCTGTGGAAT‑CT‑CGCCCAT‑ATTCTGGCGTTAGCCTTGGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTCTTTGCTC‑GCCGGCAAATTTAAAAGCGTCGGAGA‑TGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCa                                                 <  1:86198/209‑1 (MQ=255)                                                                                                                                                                          tcattaCTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCa                                                 <  2:546048/227‑1 (MQ=255)                                                                                                                                                                          tcatcaCTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGGCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCa                                                 <  2:262759/227‑1 (MQ=255)                                                                                                                                                                                cTTTGTTGCCCTCCGGCAGAATGGCGGTAATTGTC‑CTCAGCGCCTGGCGGTTCTGGGTAAT‑CT‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTta                                           <  2:522038/227‑1 (MQ=255)                                                                                                                                                                                 tttGGTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTtat                                          <  2:5957/227‑1 (MQ=255)                                                                                                                                                                                   tgttgCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCCG‑ATTCGCGCGGTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGAA‑GTTTTTGCTC‑TCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCCGTTTCCCCTCATCGCCCGCCGGGGAAAAACCAGAGTtatacg                                       >  2:382842/1‑226 (MQ=255)                                                                                                                                                                                   tgttgCCATCCGGCAGAATGGCGGTAA‑TATCACTCAGCGCCTGACGGTTCTGGGGAAT‑CA‑CCCGCAG‑GTTCGCGCGTTCGCCCTTCGGCGTGGCGGGCGTTTT‑TTTTTGTTGCA‑GTTTTTGGTC‑GCCGCCCCATTTAAAAACGGCAGTGG‑AGT‑AGATTTTTTCGGTTTTCCCATTCCCTTGACGCCATTTGCACCTACTCGCCCCCCGGGTACAGCCAAAAGtaaa                                          >  1:638550/1‑223 (MQ=255)                                                                                                                                                                                    gttgCCATCCGGCAGAATGGCGGTAA‑TGTCCCTCAGCGCCTGGCGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTa                                                        <  2:531706/210‑1 (MQ=255)                                                                                                                                                                                            ccGGGAGAATGGCGGTAA‑TTTCACTCAGCGCCTGACGGTTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTg                              <  2:476264/228‑1 (MQ=255)                                                                                                                                                                                                          ggaa‑TGTCACTCCGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCCGGCTGGGCGGTCGTTTC‑TTTTAGTTGCA‑GTTTTTGCTT‑GCCGCCCAATTTAAAAGCGTGGGAAG‑GGT‑A‑ATTTTTTCCGGTTTCTCCGCCGCGGCGGACCCATTGCCGCCCATTGCCCCCCGGGGTAAGACCCAAGAGATAGGGTTTTCCCGcctccacgggggg                 >  1:106080/3‑215 (MQ=255)                                                                                                                                                                                                           tAA‑TGTCACTCAGCCCCTGACGATTCGGGGTAAT‑CA‑CGCGCAG‑ATTCGTGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACGCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGa                <  2:84658/227‑1 (MQ=255)                                                                                                                                                                                                                gTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTa            <  2:526079/227‑1 (MQ=255)                                                                                                                                                                                                                   acccagccacTGACGATGCTGGGTAAT‑CG‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTaac          <  2:254145/217‑1 (MQ=255)                                                                                                                                                                                                                         cgcCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCGGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATTGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCCGCTTGCCCGGAGAGCTCCAACAACACTTc  >  1:271107/1‑228 (MQ=255)                                                                                                                                                                                                                         cgcCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATCGGTTTTGCCTGCTTGCCAGGAGCGGTCTAACAACACTTc  >  2:46829/1‑228 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                    CCGGGGCTTAGCTTTTCATCACCCCAGTTGCCGTCGGCGGTAATGGACTGGTAAAAGCCAATTTGGTTATTCGGACGATATTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTAACAACACTTC  >  NC_000913/4244314‑4244756

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