Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,895,556	T→C	19.2%	M352V (ATG→GTG)	adeP ←	adenine permease, high affinity; adenine:H+ symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,895,556	0	T	C	19.2%	106.6 / 3.5	54	M352V (ATG→GTG)	adeP	adenine permease, high affinity; adenine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base T (25/14);  new base C (2/10);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 6.86e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATC  >  NC_000913/3895330‑3895794                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 cACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACCCGTGCCAGGCTGGAAG‑TCATCAGCACCCCAACGTAAACAAACGCACCAGCTGC‑AGCGTAGCCTGGGCA‑CCA‑CCCCGCCAgcg                                                                                                                                                                                                                                       >  1:379072/1‑238 (MQ=255)                cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTCAAATCCCGCCA‑GATCACGCGTGCCACACTGGAAG‑ACACAACGACGCCCACCGAAATCAAGCCACCAACTCC‑AGCGAAGACC‑GGCA‑CCCACCCCCCCCCCGGCGCCacacaaaaaaa                                                                                                                                                                                                                       >  2:567746/1‑229 (MQ=255)                cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCC‑GGCA‑CCATCCCCGCCACCGGCGACCGAAAGATAAc                                                                                                                                                                                                                       >  1:177540/1‑239 (MQ=255)                      cGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGGACAGATTCAATAAGATCCTGCCA‑GTTCACCCGTGCCAAACTGGAAG‑GCATCAACACCCCAACGTAAATCAAGCCACCCCACGC‑AGCGCACCCT‑GGCC‑CCCCCCCCCCCCCCGGCGCCAGAAAAATAAcacacca                                                                                                                                                                                                                 >  2:771979/1‑233 (MQ=255)                                            cGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACCCGTGCCAGACTGGAAG‑TCAACAGCACGCCCACGTAAAACAGCGCACCAACTGC‑CGCGTAGCCC‑GGCC‑CCACCCCCGCCACCGGCGGCCGAAAGATAAACACCAGGAACAGCAGACcaacaaac                                                                                                                                                                                              >  1:830343/1‑234 (MQ=255)                                               aGTAGGAGAAAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCCCGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGaa                                                                                                                                                                                                                              <  1:548609/201‑1 (MQ=255)                                                  aGGAGTAAAAGCCTAGCGCGATACCTTGGGTAATGGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAAAAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGt                                                                                                                                                                                     <  1:267674/238‑1 (MQ=255)                                                         aaaGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCC‑CCATCCCCCCCAGCGGCGCCCGAAAGATAACCAGCAGGAACAGCAGACC‑ACAACACACTGCCGCCAAACCgg                                                                                                                                                                             >  1:221780/1‑239 (MQ=255)                                                           aGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAAATCCTGCCA‑GTTCACGCGGGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCC‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAAAAAA‑CACTGCCGTCAGACCggg                                                                                                                                                                            >  1:580011/1‑237 (MQ=255)                                                              ccgccgcattcccTTTGGTAATGGAGAA‑CTGAACGGCATCCCCCGCGGGGTAAAAAAACCCAGAAAAAGATTCAG‑AAGATCCCGGCATGTTCCCGCGGCCAAGATGGGAAGCTCA‑AACCCCCCCAACTAAAATCAGCCCCCCGGCTGA‑AGGTTAGCCT‑GCCC‑CCACCCCCGCCAGCGGCGCCAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGa                                                                                                                                                                        <  2:22414/228‑1 (MQ=255)                                                                               gaaaTGGAAAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTAAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGc                                                                                                                                                        <  1:771979/236‑1 (MQ=255)                                                                                   cccaaaaGCTGACCGGCAA‑CATCCGGGGGTTAATAAACGCCGGAAAAGATTCAGTAAGACCCTGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGaa                                                                                                                                                   <  1:25483/234‑1 (MQ=255)                                                                                     ggaaatcgtaaCGGCCT‑CACCACGCGGTTAATAAACCCGGGACAAAATTCGGTACGACCCTGCCA‑GTCCAGGCGTGCCAGGCGGGAAG‑TTACCAGCACCCCAACGAAAAACGGGCCCCAAGCTGA‑AGCGTGGCCT‑GCCA‑CCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAaga                                                                                                                                                 <  2:379072/230‑1 (MQ=255)                                                                                     gagaAGCTGAACGGCAC‑CATCACGGCGGTAATAACCGCCGGAACAGATTCAGAAAGATCCGGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCCCGCCAACGTAAATCAGCGCCCCAGCTGC‑AGCGTAGCCG‑GGCA‑CCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAaga                                                                                                                                                 <  1:586843/239‑1 (MQ=255)                                                                                         aGCGGAACGGCAT‑CATCAGGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTCCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGa                                                                                                                                              <  1:141818/238‑1 (MQ=255)                                                                                          tcTGAACGGCAA‑CACCACGGGGTTAAAAAACCCCGGAAAAGCGTCAGAAAGCTCCTGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCCCGCCAACGGAAATCAGCGCACCAGCTGA‑AGCGTGGCCG‑GGCC‑CCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACt                                                                                                                                            <  2:177540/238‑1 (MQ=255)                                                                                                                 gTAAAAAACCCCGGAACAGTTTCAGGAAGACCCTGCCA‑GTCCACGGGTGCCAGACGGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGgaag                                                                                                                       <  1:756073/238‑1 (MQ=255)                                                                                                                   aataaACGCCGGAACAAATTAAGAAAGCCCCCACCA‑GTCCCCGCGTGCCAGGCTGGAAG‑TCCCCAGCCCCCCAAGGGAAACCACGCCACCAGCGGC‑AGCTTGGCCG‑GGCC‑CCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGaagaag                                                                                                                    <  2:221780/239‑1 (MQ=255)                                                                                                                       acgcccGGAACAGATTCAGTAAGATCCGGCAA‑GTCCACGCGTGCCAGACTGGAAG‑TCATCAGCCCGCAAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTAcc                                                                                                                <  1:23932/235‑1 (MQ=255)                                                                                                                           ccGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACAGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGGTa                                                                                                            >  1:145355/1‑239 (MQ=255)                                                                                                                            cGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCCGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATTTAAGCCCTAACGTGAGAAGTACCGATaa                                                                                                           >  1:258089/1‑239 (MQ=255)                                                                                                                                    aTTCAGTAAGCTCCTGCAA‑GTTCACGCGTGCCAGCCTGGAAG‑TCCTCAGCACCCCAACGTAAATCAGCGCCCCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATaa                                                                                                                                                                                                                        <  2:105009/123‑1 (MQ=255)                                                                                                                                     tccAGTAAGATCCGGCCA‑GTCCACGCTGGCCAGACTGGAAG‑TCAACAGCCCGCAAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGaa                                                                                                      <  1:289807/233‑1 (MQ=255)                                                                                                                                        agagaagtCCTGCCA‑TTTCCCCGGCGCCAACCTGGACG‑TGAGCAGCGCGCAAACGAAAAGCAGCGCACCGGCTGC‑GGCGTAGCAT‑GCCA‑CCATCCCCGACAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CGCTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTCCCGATAAACGAACCGGTCa                                                                                               <  2:139392/232‑1 (MQ=255)                                                                                                                                         gTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCAc                                                                                              >  1:28921/1‑239 (MQ=255)                                                                                                                                           aaGCTCCTGCCA‑GTTCACGGGGGCCAGCCTGGAAG‑TCATCAGCACGCCAACGTAAACGAGCGCCCCAGCTGC‑AGCGAAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCAc                                                                                              <  1:784166/237‑1 (MQ=255)                                                                                                                                              aTCCTGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGaa                                                                                                                                                                                                              >  1:105009/1‑123 (MQ=255)                                                                                                                                              aTCCTGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGgaag                                                                                                                       <  1:781965/209‑1 (MQ=255)                                                                                                                                              aTCCTGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCAATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGCAGAACCCGTAAACGAACCGGTCACGGAag                                                                                         >  1:120465/1‑239 (MQ=255)                                                                                                                                               tCCTGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGCTAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAAAGAACTCAAAATAAGCCGTAACGGAAGAAGTACCGCTAAAGGAACCGGTCACGGagaa                                                                                        >  1:187923/1‑236 (MQ=255)                                                                                                                                                ccgcccA‑GTTCACGCGTGCCAGACGGGAAG‑TCCTCAGCCCGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCC‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAagag                                                                                       <  1:503903/235‑1 (MQ=255)                                                                                                                                                  tGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGACCCGGCCACGGAAGAGAt                                                                                     >  1:344723/1‑239 (MQ=255)                                                                                                                                                  tGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAAAGAt                                                                                     >  1:175317/1‑239 (MQ=255)                                                                                                                                                  tGCCA‑GTTAACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCACCGCCAGCGGCGACATAAAGATAACCAGCAGGAACAGCAGCCCAACAAC‑CCCTGCCGTCAGACAGGTAACACCGCCAACCTGTAAGCCTGCAGAGGCCTCAATATAAACCGCAACGGAAGAAGCAACGATAAACGAACCGGGCACGGaaacaca                                                                                     >  2:698140/1‑234 (MQ=255)                                                                                                                                                    ccA‑GTTCCCGCGGGCCAGACTGGAAG‑TCATGAGCACGCCAACGTAAACCAGCGCACCAGGTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATAc                                                                                   <  1:485427/239‑1 (MQ=255)                                                                                                                                                         ttCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCg                                                                                                                                                                     >  2:140170/1‑153 (MQ=255)                                                                                                                                                         ttCACGCGTGCAAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCg                                                                                                                                                                     <  1:140170/153‑1 (MQ=255)                                                                                                                                                          tcccggttGCCAGACTGGAAG‑TCACCGGCACGCAACCGTAAATCAGGGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCg                                                                              <  1:211445/232‑1 (MQ=255)                                                                                                                                                          ccgcgcgTGCCAGTCTGGATG‑TCACCAGCCCCAAAACGTAAATCAGCGCACCAGCTGC‑AGGGAAGCGT‑GGCA‑CAACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTg                                                                                 <  1:393705/233‑1 (MQ=255)                                                                                                                                                                   ccAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCCCGGAAGGGAAACTGTCgcacaacag                                                                      >  1:337856/1‑230 (MQ=255)                                                                                                                                                                                      gCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑ggg                                                  >  1:596341/1‑239 (MQ=255)                                                                                                                                                                                       cACGACAACGTAAATGAGCCCCCCAGCAGC‑AGGCTAGCC‑‑GGCACCCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGa                                                 <  1:12028/239‑1 (MQ=255)                                                                                                                                                                                            cAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACCGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCACTCTAAGCCGTAACGGAAGAAGTACCGCTAAACGAACCGGTCACGGACGAGATACTGTCGCCATACAACGCCTGCTTCATCCCC‑GGGAAttt                                             >  1:96015/1‑238 (MQ=255)                                                                                                                                                                                                    aTCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGCTAAACGAACCGGTCACGGAAGAGATACTGTCGCCATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTtctc                                    >  1:885110/1‑239 (MQ=255)                                                                                                                                                                                                    aTCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGCCGACATACAGCGCCTCCTTCATGCGC‑GGGAATTTCCACTcct                                     >  1:689137/1‑238 (MQ=255)                                                                                                                                                                                                        gcgcACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGCTAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCCGCTTCAAGCGC‑GGGAAATTCCCCTTCCCaac                                 >  2:272930/1‑236 (MQ=255)                                                                                                                                                                                                         cgcACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTTCTCATcc                                >  1:514246/1‑238 (MQ=255)                                                                                                                                                                                                         cgcACAAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCATCAGGAACAGCAGACCACCAAC‑CACTGCCGACAGACCGGTACGACCGCCAACCGCAACGCCGGCAGAGGCCCCACTATAAGCCGGACAGGAAGAACTCCCGACACACGACCCGGCCACGGAAGGAACCCTTGCGCCATACAACGCCCGCTCCCTTCGC‑GGGCATTTCCCCCTCCCAAcc                                >  2:772742/1‑238 (MQ=255)                                                                                                                                                                                                          gcACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGCAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGCCACGGCAGAGATACTGTAGACAAACAACGCCTGCTTCAT‑CGCGGGGAATTTCCCCTCCCCACCCGc                              >  2:688555/1‑239 (MQ=255)                                                                                                                                                                                                          gcACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAAAGATACTGTCGACATACAGCGCCTGCTTCATCCCC‑GGGAATTTCCCCTTCCCATccc                               >  2:48077/1‑237 (MQ=255)                                                                                                                                                                                                          gcACCAG‑TGCAAGCGTAGCCT‑GGCA‑CCACCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTTCTCATCCGc                              <  1:641166/239‑1 (MQ=255)                                                                                                                                                                                                                  ggC‑AGCTTAGCCT‑GGCC‑CCCGCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CCCTGCCGTCAGACCGGTACGACCGCCAACCGGTACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTTCTCATCCGCCAGAc                         <  1:551062/235‑1 (MQ=255)                                                                                                                                                                                                                            gCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCgg             >  1:745884/1‑239 (MQ=255)                                                                                                                                                                                                                               t‑GGCA‑CCATCCCCGCCAGCGGCGAGAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGACGTCAGGCCGGTACGACCGCCAACCGATTCGCCGGAAGAGCAATCAATGTAAGCAGCAACGAGAGAAGTACCGATAACCGAACCGGCCACGGCAGAGAAACATTCGCGACACAACGCTTACATCAGGCCC‑CGGAATATCACATCCTCCAACGCCAGCACGGCTGTATCggga           >  2:67087/1‑236 (MQ=255)                                                                                                                                                                                                                                         tCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGATATACTGTCGACATACAGCGCCGGCTTCATGCGC‑GGGCATTTCCCCTTCCCATCCGCCAGACCTGCTTTATCGGTCACGCCAATc  >  1:134618/1‑239 (MQ=255)                                                                                                                                                                                                                                         tCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATc  >  1:528918/1‑239 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 CACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGGAAG‑TCATCAGCACGCCAACGTAAATCAGCGCACCAGCTGC‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCGACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATC  >  NC_000913/3895330‑3895794

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