Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,244,521	T→G	17.1%	D185A (GAC→GCC)	malF ←	maltose transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,244,521	0	T	G	17.1%	90.4 / 2.4	53	D185A (GAC→GCC)	malF	maltose transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (16/27);  new base G (7/3);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 8.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTGTAACCGGGGCTTAGCTTTTCATCACCCCAGTTGCCGTCGGCGGTAATGGACTGGTAAAAGCCAATTTGGTTATTCGGACGATATTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTA  >  NC_000913/4244308‑4244746                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           ggtctatcgggtttttgtttttCATCACCCGCGGGGGGGACGGAGGTAATGGACTGGTAAAAGCCAATGTGGGTATGCGGGCGAAATTTCACGACGTTCCGATTATTCGTCAACGGTCCGTCCCCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGTCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGa                                                                                                                                                                                                                              <  1:335533/210‑1 (MQ=255)       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)                                                         gTAAAAGCCAATTTGGTTATTCGGACGATATTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTg                                                                                                                                                                                                                               >  2:459220/1‑171 (MQ=255)                                                          tAAAAGCCAATGTGGTTATTCGGACGATATTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGa                                                                                                                                                                                                                              <  1:459220/171‑1 (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)                                                                                               cgcTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCatc                                                                                                                                                                                                                                    >  1:407985/1‑126 (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)                                                                                                gcTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAgcg                                                                                                                                                                                                                                   <  2:407985/128‑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)                                                                                                                                                                                                                       |                                                                                                                                                                                                                                           GTGTAACCGGGGCTTAGCTTTTCATCACCCCAGTTGCCGTCGGCGGTAATGGACTGGTAAAAGCCAATTTGGTTATTCGGACGATATTTCACGCCGCTCTGATTATTCGTCAACGTGCCGTCACCGTCGAGTGTGTAGAGCGGCTGCGTGCCAGAAAACTGGCGCAGGGAGCTCATCATCACTTTGTTGCCATCCGGCAGAATGGCGGTAA‑TGTCACTCAGCGCCTGACGATTCTGGGTAAT‑CA‑CGCGCAG‑ATTCGCGCGTTCGCCTTCGGGCTGGGCGGTCGTTTC‑TTTCAGTTGCA‑GTTTTTGCTC‑GCCGCCAAATTTAAAAGCGTCGGAGA‑GGT‑A‑ATTTTTGCCGGTTTCGCCGTCGCTGAGCGCCAGTTGCCACTCATCGCCCGCCGGGTAAAGACCAAAGTTATAGGTTTTGCCTGCTTGCCAGGAGCGATCTA  >  NC_000913/4244308‑4244746

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