Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,895,555	T→C	19.0%	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,555	0	T	C	19.0%	125.3 / 4.3	71	M352V (ATG→GTG)	adeP	adenine permease, high affinity, adenine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base T (32/17);  new base C (3/12);  total (41/30)
Fisher's exact test for biased strand distribution p-value = 2.84e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATC  >  NC_000913/3895329‑3895793                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  cacgggattagagtccggagcgttCCGGTGCCTATTGGCATCACGCGGGAGGAAAGAAAGCCCAGGGCGATCCCTGCGGAAATCGAGAAGCTGACGGGCAA‑CATGACGGCGGTAATAAACGCAGAAACAGATTCAGGAAGAACCTGCGA‑GTTCACGCGTGCCAGGCTGG‑AAG‑TCATCAGCCCGCAAACGTAATTGAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAgcgg                                                                                                                                                                                                                                       <  1:785973/216‑1 (MQ=255) cACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACCCGTGCCAGGCTGG‑AAG‑TCATCAGCACCCCAACGTAAACAAACGCA‑CCAGC‑TG‑C‑AGCGTAGCCTGGGCA‑CCA‑CCCCGCCAgcg                                                                                                                                                                                                                                        >  2:379072/1‑238 (MQ=255)                cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTCAAATCCCGCCA‑GATCACGCGTGCCACACTGG‑AAG‑ACACAACGACGCCCACCGAAATCAAGCCA‑CCAAC‑TC‑C‑AGCGAAGACC‑GGCA‑CCCACCCCCCCCCCGGCG‑CCacacaaaaaaa                                                                                                                                                                                                                       >  1:567746/1‑229 (MQ=255)                cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCC‑GGCA‑CCATCCCCGCCACCGGCG‑ACCGAAAGATAAc                                                                                                                                                                                                                       >  2:177540/1‑239 (MQ=255)                cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTAAGTAAGATCCTGCCA‑GTTCACCCGTGCCAAACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TT‑C‑AGCGTACCCT‑GGCC‑CCACCCCCCCCAccggggccgagaagagaac                                                                                                                                                                                                                        >  2:590289/1‑224 (MQ=255)                cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGAATCAGTAAGATCCTTCCC‑GTTTAACCGCGCCAGAATTG‑AAG‑GCCTACGACCGCCAACCGAAATCAGCGCA‑CCAGC‑TG‑A‑AGCGCAACCT‑GGCA‑CCAACccccccccgggggaaacaaaaaaaaa                                                                                                                                                                                                                        >  1:267674/1‑219 (MQ=255)                cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAAATTCAGTAAGATCCTGCCC‑ATTCACGCGTGCCAAACTGG‑AAA‑GTATCAGCACGCCCACCCAAATCAACCCC‑CCAAC‑TC‑C‑AGCGTAGCCT‑GGCA‑CCCACCCCCCCAgcggcgccgaaaaaaaaac                                                                                                                                                                                                                        >  1:756073/1‑226 (MQ=255)                cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGGAATAAACGCCGGAAACGATTCAATAAGATCCTGCCA‑GTTCACCCGTGCCAGGCTGG‑AAG‑TCATCAGAACCCCAACGCAAATCAGCCGC‑CCAGCGTC‑A‑AGGGTA‑CCT‑GGCA‑CCAACccccccccccgcggcacgaagcaaac                                                                                                                                                                                                                        >  2:590348/1‑219 (MQ=255)                cGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGGAATAAACACCGGAACAGAATCAAGAAAAACCTGCCA‑GCTCACGCGTGCCCGACTGG‑AAG‑TCATACGCACGCCAACCGAAATAAGGCCA‑CCCAC‑TG‑C‑AACGAACCCT‑GGCA‑CCAACCCCCCCCgcggcggacaaaagagaaa                                                                                                                                                                                                                        >  1:691727/1‑226 (MQ=255)                      cGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGGACAGATTCAATAAGATCCTGCCA‑GTTCACCCGTGCCAAACTGG‑AAG‑GCATCAACACCCCAACGTAAATCAAGCCA‑CCCCA‑CG‑C‑AGCGCACCCT‑GGCC‑CCCCCCCCCCCCCCGGCG‑CCAGAAAAATAAcacacca                                                                                                                                                                                                                 >  1:771979/1‑233 (MQ=255)                                 gATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAAAGCCCGAACAGATTCAGTAAAATCCTGCCA‑GGTCACGCGTGCCAGACTGG‑AAG‑TCCTCAGCACGCCAACGCAAATCAGCGCA‑CCCCC‑TG‑C‑AGCCGAGCCT‑GGCC‑CCATCCCCGCCACCCGCG‑ACCAAAAGAAAAAcacaagaaaaaacaaac                                                                                                                                                                                                      >  1:141818/1‑224 (MQ=255)                                  aTTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAAAGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACCCGTGCCAGACTGG‑AAG‑TCATCCGCACGCCAACCTAAATCAGCGCA‑CCAGC‑TGCC‑AGCG‑AGCCT‑TGCA‑CCATCCCCGCCACCGGCG‑GCAGAAAGGTaaaccaccgggcaagcaaacc                                                                                                                                                                                                     >  2:27336/1‑220 (MQ=255)                                            cGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACCCGTGCCAGACTGG‑AAG‑TCAACAGCACGCCCACGTAAAACAGCGCA‑CCAAC‑TG‑C‑CGCGTAGCCC‑GGCC‑CCACCCCCGCCACCGGCG‑GCCGAAAGATAAACACCAGGAACAGCAGACcaacaaac                                                                                                                                                                                              >  2:830343/1‑234 (MQ=255)                                               aGTAGGAGAAAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCCCGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGaa                                                                                                                                                                                                                              <  2:548609/201‑1 (MQ=255)                                                 tAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGGAATAAACGCAGGGACAAAATCAGTAAGATCCTGCCA‑GCTCACACGGGCCAGACTGG‑GAG‑TCATCCACCAGCCAACGTAACTCAAACCA‑CCCAC‑TG‑C‑GCGGTAAACC‑GGGC‑CACTCCCCCCCAACGGCG‑AACGAAAAATAACCACCAGGAAAAGCAAAcacaaaaacaacaccgc                                                                                                                                                                                      >  1:393705/1‑223 (MQ=255)                                                  aGGAGTAAAAGCCTAGCGCGATACCTTGGGTAATGGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAAAAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGt                                                                                                                                                                                     <  2:267674/238‑1 (MQ=255)                                                         aaaGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCC‑CCATCCCCCCCAGCGGCG‑CCCGAAAGATAACCAGCAGGAACAGCAGACC‑ACAACACACTGCCGCCAAACCgg                                                                                                                                                                             >  2:221780/1‑239 (MQ=255)                                                           aGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAAATCCTGCCA‑GTTCACGCGGGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCC‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAAAAAA‑CACTGCCGTCAGACCggg                                                                                                                                                                            >  2:580011/1‑237 (MQ=255)                                                           aGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAAAGCCGGAACAGATTCAGTAAGATACTGCCC‑GTTCACCCGTGCCAGACTGG‑AAG‑GCATCAGCACGCCCACGCAAAACAGCCCACCCCGC‑TC‑C‑A‑CCTAGCCC‑GGCA‑CCATCCCCCCCCGCGGCGCACA‑AAAGAAAAACAGCCGGaacaacaaaacaacaaaacaccgccccaaaccggga                                                                                                                                                                            >  1:586843/1‑219 (MQ=255)                                                              ccgccgcattcccTTTGGTAATGGAGAA‑CTGAACGGCATCCCCCGCGGGGTAAAAAAACCCAGAAAAAGATTCAG‑AAGATCCCGGCATGTTCCCGCGGCCAAGATGGG‑AAGCTCA‑AACCCCCCCAACTAAAATCAGCCCC‑CCGGC‑TG‑A‑AGGTTAGCCT‑GCCC‑CCACCCCCGCCAGCGGCG‑CCAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGa                                                                                                                                                                        <  1:22414/228‑1 (MQ=255)                                                                               gaaaTGGAAAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTAAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGc                                                                                                                                                        <  2:771979/236‑1 (MQ=255)                                                                                   cccaaaaGCTGACCGGCAA‑CATCCGGGGGTTAATAAACGCCGGAAAAGATTCAGTAAGACCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGaa                                                                                                                                                   <  2:25483/234‑1 (MQ=255)                                                                                     ggaaatcgtaaCGGCCT‑CACCACGCGGTTAATAAACCCGGGACAAAATTCGGTACGACCCTGCCA‑GTCCAGGCGTGCCAGGCGGG‑AAG‑TTACCAGCACCCCAACGAAAAACGGGCCC‑CAAGC‑TG‑A‑AGCGTGGCCT‑GCCA‑CCACCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAaga                                                                                                                                                 <  1:379072/230‑1 (MQ=255)                                                                                     gagaAGCTGAACGGCAC‑CATCACGGCGGTAATAACCGCCGGAACAGATTCAGAAAGATCCGGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCCCGCCAACGTAAATCAGCGCC‑CCAGC‑TG‑C‑AGCGTAGCCG‑GGCA‑CCACCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAaga                                                                                                                                                 <  2:586843/239‑1 (MQ=255)                                                                                         aGCGGAACGGCAT‑CATCAGGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTCCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGa                                                                                                                                              <  2:141818/238‑1 (MQ=255)                                                                                          tcTGAACGGCAA‑CACCACGGGGTTAAAAAACCCCGGAAAAGCGTCAGAAAGCTCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCCCGCCAACGGAAATCAGCGCA‑CCAGC‑TG‑A‑AGCGTGGCCG‑GGCC‑CCACCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACt                                                                                                                                            <  1:177540/238‑1 (MQ=255)                                                                                           cgtaaccgccaacctcccggcgaGAATAAACCCCGGAACAGATTCAGAAAGACCCCGCCA‑GTCCACGCGTGCCAGACTGA‑AAG‑CCACCACCACGCCACCGTAAATCGGCGAC‑CCAGC‑TG‑C‑AGCGTGGCCG‑GGCC‑CAATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACt                                                                                                                                            <  1:590289/222‑1 (MQ=255)                                                                                             tgaacggcctcccacggggggTAAAAAACGGCGAAACAGATTCAAAAAACCCCGGCCG‑TTCCACGCGTGCAAGACGGG‑AAG‑TCATCCGCCCCCCAACGTAAAACGGCGCA‑CCAGC‑GG‑C‑AGCGTGGCCG‑GGCC‑CCACCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTc                                                                                                                                           <  1:27336/219‑1 (MQ=255)                                                                                                                 gTAAAAAACCCCGGAACAGTTTCAGGAAGACCCTGCCA‑GTCCACGGGTGCCAGACGGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGgaag                                                                                                                       <  2:756073/238‑1 (MQ=255)                                                                                                                   aataaACGCCGGAACAAATTAAGAAAGCCCCCACCA‑GTCCCCGCGTGCCAGGCTGG‑AAG‑TCCCCAGCCCCCCAAGGGAAACCACGCCA‑CCAGC‑GG‑C‑AGCTTGGCCG‑GGCC‑CCACCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGaagaag                                                                                                                    <  1:221780/239‑1 (MQ=255)                                                                                                                       acgcccGGAACAGATTCAGTAAGATCCGGCAA‑GTCCACGCGTGCCAGACTGG‑AAG‑TCATCAGCCCGCAAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCACCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTAcc                                                                                                                <  2:23932/235‑1 (MQ=255)                                                                                                                           ccGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACAGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGGTa                                                                                                            >  2:145355/1‑239 (MQ=255)                                                                                                                            cGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCCGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATTTAAGCCCTAACGTGAGAAGTACCGATaa                                                                                                           >  2:258089/1‑239 (MQ=255)                                                                                                                                    aTTCAGTAAGCTCCTGCAA‑GTTCACGCGTGCCAGCCTGG‑AAG‑TCCTCAGCACCCCAACGTAAATCAGCGCC‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATaa                                                                                                                                                                                                                        <  1:105009/123‑1 (MQ=255)                                                                                                                                    aTTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCCACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGCAGCCT‑GGCC‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAAACAC‑CACGCCCGTCAAACCGGTACACACACCACCCGCAGCGCCCCACAATGACCCACTATAACCCGTCACGgaagaaacacacaaaaacaaaaag                                                                                                   >  1:381462/1‑221 (MQ=255)                                                                                                                                     tccAGTAAGATCCGGCCA‑GTCCACGCTGGCCAGACTGG‑AAG‑TCAACAGCCCGCAAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGaa                                                                                                      <  2:289807/233‑1 (MQ=255)                                                                                                                                        agagaagtCCTGCCA‑TTTCCCCGGCGCCAACCTGG‑ACG‑TGAGCAGCGCGCAAACGAAAAGCAGCGCA‑CCGGC‑TG‑C‑GGCGTAGCAT‑GCCA‑CCATCCCCGACAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CGCTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTCCCGATAAACGAACCGGTCa                                                                                               <  1:139392/232‑1 (MQ=255)                                                                                                                                         gTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCAc                                                                                              >  2:28921/1‑239 (MQ=255)                                                                                                                                           aaGCTCCTGCCA‑GTTCACGGGGGCCAGCCTGG‑AAG‑TCATCAGCACGCCAACGTAAACGAGCGCC‑CCAGC‑TG‑C‑AGCGAAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCAc                                                                                              <  2:784166/237‑1 (MQ=255)                                                                                                                                            aaaacgtcccccattCACGCGTGGCAGACAGGAAAG‑CCACACCCCCGCAACCGTAAATCAGCCCA‑CCACG‑GG‑G‑AGCGTACCCG‑GGCC‑CCACCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CCCTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCAc                                                                                              <  1:618402/225‑1 (MQ=255)                                                                                                                                             ataccccgcaaggtCACGCGTGCCAGACGGG‑AAG‑CCATCGACAGCCCAACGTAAACCAGCACT‑ACAGC‑TG‑C‑AGCGAAGCCT‑GGCA‑CCATCCCCGACAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGa                                                                                           <  2:162623/226‑1 (MQ=255)                                                                                                                                              aTCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGaa                                                                                                                                                                                                              >  2:105009/1‑123 (MQ=255)                                                                                                                                              aTCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGgaag                                                                                                                       <  2:781965/209‑1 (MQ=255)                                                                                                                                              aTCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCAATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGCAGAACCCGTAAACGAACCGGTCACGGAag                                                                                         >  2:120465/1‑239 (MQ=255)                                                                                                                                               tCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGCTAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAAAGAACTCAAAATAAGCCGTAACGGAAGAAGTACCGCTAAAGGAACCGGTCACGGagaa                                                                                        >  2:187923/1‑236 (MQ=255)                                                                                                                                                ccgcccA‑GTTCACGCGTGCCAGACGGG‑AAG‑TCCTCAGCCCGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCC‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAagag                                                                                       <  2:503903/235‑1 (MQ=255)                                                                                                                                                  tGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGACCCGGCCACGGAAGAGAt                                                                                     >  2:344723/1‑239 (MQ=255)                                                                                                                                                  tGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAAAGAt                                                                                     >  2:175317/1‑239 (MQ=255)                                                                                                                                                  tGCCA‑GTTAACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCACCGCCAGCGGCG‑ACATAAAGATAACCAGCAGGAACAGCAGCCCAACAAC‑CCCTGCCGTCAGACAGGTAACACCGCCAACCTGTAAGCCTGCAGAGGCCTCAATATAAACCGCAACGGAAGAAGCAACGATAAACGAACCGGGCACGGaaacaca                                                                                     >  1:698140/1‑234 (MQ=255)                                                                                                                                                    ccA‑GTTCCCGCGGGCCAGACTGG‑AAG‑TCATGAGCACGCCAACGTAAACCAGCGCA‑CCAGG‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATAc                                                                                   <  2:485427/239‑1 (MQ=255)                                                                                                                                                         ttCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCg                                                                                                                                                                     >  1:140170/1‑153 (MQ=255)                                                                                                                                                         ttCACGCGTGCAAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCg                                                                                                                                                                     <  2:140170/153‑1 (MQ=255)                                                                                                                                                          tcccggttGCCAGACTGG‑AAG‑TCACCGGCACGCAACCGTAAATCAGGGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCACCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCg                                                                              <  2:211445/232‑1 (MQ=255)                                                                                                                                                          ccgcgcgTGCCAGTCTGG‑ATG‑TCACCAGCCCCAAAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGGGAAGCGT‑GGCA‑CAACCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTg                                                                                 <  2:393705/233‑1 (MQ=255)                                                                                                                                                                   ccAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCCCGGAAGGGAAACTGTCgcacaacag                                                                      >  2:337856/1‑230 (MQ=255)                                                                                                                                                                                   atcaGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACCAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGCGGACTCAATATAAACCGTACCGGAAGAAATACCGAAAAACAAACCGGTCACGGAACAGACACAGTAAACAAACAgaacacgcctcaagcg                                                       >  1:395073/1‑224 (MQ=255)                                                                                                                                                                                       gCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑ggg                                                  >  2:596341/1‑239 (MQ=255)                                                                                                                                                                                        cACGACAACGTAAATGAGCCCC‑CCAGC‑AG‑C‑AGGCTAGCC‑‑GGCACCCACCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGa                                                 <  2:12028/239‑1 (MQ=255)                                                                                                                                                                                             cAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACCGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCACTCTAAGCCGTAACGGAAGAAGTACCGCTAAACGAACCGGTCACGGACGAGATACTGTCGCCATACAACGCCTGCTTCATCCCC‑GGGAAttt                                             >  2:96015/1‑238 (MQ=255)                                                                                                                                                                                                     aTCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGCTAAACGAACCGGTCACGGAAGAGATACTGTCGCCATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTtctc                                    >  2:885110/1‑239 (MQ=255)                                                                                                                                                                                                     aTCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGCCGACATACAGCGCCTCCTTCATGCGC‑GGGAATTTCCACTcct                                     >  2:689137/1‑238 (MQ=255)                                                                                                                                                                                                         gcgcA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGCTAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCCGCTTCAAGCGC‑GGGAAATTCCCCTTCCCaac                                 >  1:272930/1‑236 (MQ=255)                                                                                                                                                                                                          cgcA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTTCTCATcc                                >  2:514246/1‑238 (MQ=255)                                                                                                                                                                                                          cgcA‑CAAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCATCAGGAACAGCAGACCACCAAC‑CACTGCCGACAGACCGGTACGACCGCCAACCGCAACGCCGGCAGAGGCCCCACTATAAGCCGGACAGGAAGAACTCCCGACACACGACCCGGCCACGGAAGGAACCCTTGCGCCATACAACGCCCGCTCCCTTCGC‑GGGCATTTCCCCCTCCCAAcc                                >  1:772742/1‑238 (MQ=255)                                                                                                                                                                                                           gcA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGCAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGCCACGGCAGAGATACTGTAGACAAACAACGCCTGCTTCAT‑CGCGGGGAATTTCCCCTCCCCACCCGc                              >  1:688555/1‑239 (MQ=255)                                                                                                                                                                                                           gcA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAAAGATACTGTCGACATACAGCGCCTGCTTCATCCCC‑GGGAATTTCCCCTTCCCATccc                               >  1:48077/1‑237 (MQ=255)                                                                                                                                                                                                           gcA‑CCAG‑‑TG‑CAAGCGTAGCCT‑GGCA‑CCACCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTTCTCATCCGc                              <  2:641166/239‑1 (MQ=255)                                                                                                                                                                                                                     gg‑C‑AGCTTAGCCT‑GGCC‑CCCGCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CCCTGCCGTCAGACCGGTACGACCGCCAACCGGTACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTTCTCATCCGCCAGAc                         <  2:551062/235‑1 (MQ=255)                                                                                                                                                                                                                                gCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCgg             >  2:745884/1‑239 (MQ=255)                                                                                                                                                                                                                                   t‑GGCA‑CCATCCCCGCCAGCGGCG‑AGAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGACGTCAGGCCGGTACGACCGCCAACCGATTCGCCGGAAGAGCAATCAATGTAAGCAGCAACGAGAGAAGTACCGATAACCGAACCGGCCACGGCAGAGAAACATTCGCGACACAACGCTTACATCAGGCCC‑CGGAATATCACATCCTCCAACGCCAGCACGGCTGTATCggga           >  1:67087/1‑236 (MQ=255)                                                                                                                                                                                                                                   t‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAGAAGGACTCAATATAAGCCGTAACGGAACAAGAACCGATAAACGAAACGGTCACGAACGGGAAACGTCCGGCTAACAGCGCCAGCTTCATGCGC‑GGGAATTTCCCCTTCTCAccccccccacaagcctaaccgatcg          >  1:462279/1‑220 (MQ=255)                                                                                                                                                                                                                                             tCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGATATACTGTCGACATACAGCGCCGGCTTCATGCGC‑GGGCATTTCCCCTTCCCATCCGCCAGACCTGCTTTATCGGTCACGCCAATc  >  2:134618/1‑239 (MQ=255)                                                                                                                                                                                                                                             tCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATc  >  2:528918/1‑239 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  CACGGGCTTAGATCACGCAGACGTCCGGTACCGATTTTCATCACGCAGTAGGAGATAAAGCCCAGCGCGATACCTTCGGTAATCGAGAAGCTGAACGGCAT‑CATCACGGCGGTAATAAACGCCGGAACAGATTCAGTAAGATCCTGCCA‑GTTCACGCGTGCCAGACTGG‑AAG‑TCATCAGCACGCCAACGTAAATCAGCGCA‑CCAGC‑TG‑C‑AGCGTAGCCT‑GGCA‑CCATCCCCGCCAGCGGCG‑ACAGAAAGATAACCAGCAGGAACAGCAGACCAACAAC‑CACTGCCGTCAGACCGGTACGACCGCCAACCGATACGCCGGAAGAGGACTCAATATAAGCCGTAACGGAAGAAGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGC‑GGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATC  >  NC_000913/3895329‑3895793

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