Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,931,779	T→C	F155S (TTT→TCT)	kup →	K(+):H(+) symporter Kup

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,931,779	0	T	C	98.0%	198.1 / ‑6.4	51	F155S (TTT→TCT)	kup	K(+):H(+) symporter Kup
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (25/25);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCCGGTGAGGGGGGGATCCTGACGTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATTTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGTGTTGTCGATTACGGGGGTCGAGGCGCTGTATGCTGATATGGGG  >  NC_000913/3931544‑3932023                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     aaCGCCGGTGAGGGGGGGATCCTGACGTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGACTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATTTATGATTCAAAAAc                                                                                                                                                                                                                                        >  2:414156/1‑250 (MQ=255)   cGCCGGTGAGGGGGGGATCCTGACGTTGATGTCGCTTGCCGGGCGTAAGACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGACCCGAAGCTGGAGACCTGGCTAGTTCCCCTCTAAATTATCGTTCTCACTTTCTTATCTATAATTCAAAAACATg                                                                                                                                                                                                                                     >  2:26353/1‑251 (MQ=255)       ggTGAGGGGGGGATCCTGACGTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAGCCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCGTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGTTCTTGAAATCGTCGCCCCGCAGCTTGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCAc                                                                                                                                                                                                                                 <  1:217883/251‑1 (MQ=255)             ggggggATCCTGACGTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTAt                                                                                                                                                                                                                           >  1:300695/1‑251 (MQ=255)                         aCGTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCt                                                                                                                                                                                                               >  1:204353/1‑251 (MQ=255)                           gTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGt                                                                                                                                                                                                             >  2:114924/1‑251 (MQ=255)                            ttGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGtt                                                                                                                                                                                                            >  1:59267/1‑251 (MQ=255)                                   cGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCg                                                                                                                                                                                                     <  1:511396/251‑1 (MQ=255)                                              cGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTgcc                                                                                                                                                                                          >  1:434884/1‑249 (MQ=255)                                                tAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCt                                                                                                                                                                                        >  1:151658/1‑251 (MQ=255)                                                 aaTACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCt                                                                                                                                                                                        <  1:91590/250‑1 (MQ=255)                                                                 ccACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGTATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCa                                                                                                                                                                       <  1:368174/251‑1 (MQ=255)                                                                           gCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGt                                                                                                                                                             <  2:59267/251‑1 (MQ=255)                                                                            cTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGtt                                                                                                                                                            >  1:253024/1‑251 (MQ=255)                                                                            cTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGtt                                                                                                                                                            >  1:395975/1‑251 (MQ=255)                                                                                 gATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACCCCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGt                                                                                                                                                        <  1:12569/250‑1 (MQ=255)                                                                                           cTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGAt                                                                                                                                                                                                   >  1:66486/1‑197 (MQ=255)                                                                                           cTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGAt                                                                                                                                                                                                   <  2:66486/197‑1 (MQ=255)                                                                                                  gcggcAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTGCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCg                                                                                                                                      <  1:471788/251‑1 (MQ=255)                                                                                                     gcAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGaa                                                                                                                                   <  2:395975/251‑1 (MQ=255)                                                                                                      cAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCCCTTTATTATCTATGATTCAAAAACATGGCACGGCTCTGGTCGGTAAGCTGTTTg                                                                                                                                                                                                          >  1:195937/1‑179 (MQ=255)                                                                                                      cAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTg                                                                                                                                                                                                          <  2:195937/179‑1 (MQ=255)                                                                                                             ttCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGtgctgc                                                                                                                            <  1:62924/250‑1 (MQ=255)                                                                                                                 aTGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATgc                                                                                                                        <  2:204353/250‑1 (MQ=255)                                                                                                                   ggTGAAGGCGACATAACACCCGCTATTTCGGTGATGGCAGCCACTGAAGTTCTGGAAATCGTCGCCCCGCATCTGGATACCTGGATAGTTCCCCTCGCAATTATCGTACTCACTTTATTATCTATGAGTCAAAAACAAGGCACTGCTAGGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTTGTTTCTGATTCTGGTATGACTGTGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCt                                                                                                                     <  1:489032/251‑1 (MQ=255)                                                                                                                      gAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGaa                                                                                                                  <  1:115441/251‑1 (MQ=255)                                                                                                                       aaGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGtgct                                                                                                                              <  2:123568/238‑1 (MQ=255)                                                                                                                       aaGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGtgct                                                                                                                              >  1:123568/1‑238 (MQ=255)                                                                                                                          tcgtcATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATcc                                                                                                               >  2:146402/1‑250 (MQ=255)                                                                                                                           cgtcATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCa                                                                                                              >  1:395809/1‑250 (MQ=255)                                                                                                                                 aaCACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTGTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTggg                                                                                                       <  1:319412/251‑1 (MQ=255)                                                                                                                                      ccGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGTAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTATCCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTg                                                                                                  <  2:285641/251‑1 (MQ=255)                                                                                                                                         cTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCAt                                                                                               >  2:420147/1‑251 (MQ=255)                                                                                                                                                   gATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTg                                                                                     >  1:260548/1‑251 (MQ=255)                                                                                                                                                          gCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCTCTGAATCCAATGTGGGCGTTGCATTGCTTCCTTTAATACaa                                                                              <  1:425144/251‑1 (MQ=255)                                                                                                                                                           ccATTGAAGGTCTTTAAATCGTCGCCCCGCAGCTGGATTCCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTTAATCCAATGTGGGCGGTGCTTTTCTTCCTTGAATACaaa                                                                             <  1:465123/251‑1 (MQ=255)                                                                                                                                                                            aTCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGc                                                            <  1:195060/251‑1 (MQ=255)                                                                                                                                                                              cgtcgCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCAGTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCAt                                                          <  1:323864/251‑1 (MQ=255)                                                                                                                                                                                  gCCCCGCAGCTGGATACCTGGATAGTCCCCCTCTCAATTCTCGTTCTCACTTTATTATCTATGCTTCACAAACATGGCACGGATATGGTCGGTACGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGCCTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGACTCCAATGTGGGCGGTGCATTTCTTCCTTGAGTACACAACGGTTTCTTTTATTGCATTAgg                                                      >  2:477714/1‑251 (MQ=255)                                                                                                                                                                                       gCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGc                                                   >  1:142769/1‑249 (MQ=255)                                                                                                                                                                                         aGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAgtg                                               >  2:66480/1‑251 (MQ=255)                                                                                                                                                                                                 aCCTGGATAGATCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTTGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTGCTTTTATTGCATTAGGGGCAGTGGTGTTGTc                                       <  2:260548/251‑1 (MQ=255)                                                                                                                                                                                                     ggATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGTCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTTGTGTTGTCGAtt                                   <  1:52813/251‑1 (MQ=255)                                                                                                                                                                                                                     aaTTATCGTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGTGTTGTCGATTACGGGGGTCGAGgcgc                   >  2:129768/1‑251 (MQ=255)                                                                                                                                                                                                                      aTTATCGTTCTCTCTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGGCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCAGTTCTTCCTTGAATACAAAACGGTTTCTTTTAGTGGATTGGGGGCAGTGGTGTTGTCGATTACGGGGGTCGAGGCGCt                  <  2:243915/251‑1 (MQ=255)                                                                                                                                                                                                                          tagTTCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGTGTTGTCGATTACGGGGGTCGAGGCGCTGTAt              >  2:450572/3‑251 (MQ=255)                                                                                                                                                                                                                             ttCTCACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGTGTTGTCGATTACGGGGGTCGAGGCGCTGTATGCt           >  1:191908/1‑251 (MQ=255)                                                                                                                                                                                                                              tctcACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGTGTTGTCGATTACGGGGGTCGAGGCGCTGTATGCTg          >  2:347459/1‑251 (MQ=255)                                                                                                                                                                                                                                 cACTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGTGTTGTCGATTACGGGGGTCGAGGCGCTGTATGCTg          <  2:445093/248‑1 (MQ=255)                                                                                                                                                                                                                                   cTTTATTATCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGTAAGCGCTGCATGCGCTGAATCCAATGTGAGCGGTGCATTTCTTCCTTGCATACAACACGGTTTCTTTTATTGCATTAGGCGCAGTGTTGTTGTCGATTACGTGGGTCGTGGCGCTGTATGCTGatatt     >  1:345138/1‑250 (MQ=255)                                                                                                                                                                                                                                      tattatCTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGTGTTGTCGATTACGGGGGTCGAGGCGCTGTATGCTGATATgggg  >  1:523787/1‑251 (MQ=255)                                                                                                                                                                                                                                        ttatctatGATTCAAAAACATGGCACGGTTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTaccc                                                                                                                                        <  2:435833/115‑4 (MQ=255)                                                                                                                                                                                                                                        ttatctatGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAAcc                                                                                                                                        >  1:435833/1‑115 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                     AACGCCGGTGAGGGGGGGATCCTGACGTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATTTATGATTCAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGTGTTGTCGATTACGGGGGTCGAGGCGCTGTATGCTGATATGGGG  >  NC_000913/3931544‑3932023

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑