Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,809,085	G→C	A182A (GCG→GCC)	waaA →	KDO transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,809,085	0	G	C	98.2%	164.3 / ‑4.5	57	coding (546/1278 nt)	waaA	KDO transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (29/27);  minor base . (1/0);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

AGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATT  >  NC_000913/3808844‑3809323                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                aGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATgaa                                                                                                                                                                                                                                       <  1:124840/251‑1 (MQ=255)      gTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATgaaga                                                                                                                                                                                                                                     >  1:429876/1‑248 (MQ=255)         tATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATg                                                                                                                                                                                                                                <  1:89128/250‑1 (MQ=255)              gCCGTATGACCTGCCCGAGGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCGAACCTGATTGCGGCGTTACATAAACGTAAAATTCCGCTGGTGAGCGCTAACGCGCGACGCTCTGCCCGCTCGGCCGCAGGTGATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTg                                                                                                                                                                                                                             <  2:209947/248‑1 (MQ=255)               ccGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AAt                                                                                                                                                                                                                                          >  2:251609/1‑234 (MQ=255)               ccGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AAt                                                                                                                                                                                                                                          <  1:251609/234‑1 (MQ=255)               ccGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACg                                                                                                                                                                                                                         >  2:408724/1‑251 (MQ=255)                  tATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGttt                                                                                                                                                                                                                      >  2:311668/1‑251 (MQ=255)                   aTGACCTGCCCGATGCGCTCAACCGTCTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGGTATGCCAAACTGGGAAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGtttt                                                                                                                                                                                                                     >  2:177895/1‑251 (MQ=255)                             cGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTgg                                                                                                                                                                                                           <  1:371931/251‑1 (MQ=255)                             cGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTgg                                                                                                                                                                                                           <  1:341556/251‑1 (MQ=255)                                      cAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGGCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCaaa                                                                                                                                                                                                    <  1:408724/249‑1 (MQ=255)                                                    aTAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGATCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAAcc                                                                                                                                                                                    <  2:230792/251‑1 (MQ=255)                                                       aaGTCGACCCTAAACGGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGtt                                                                                                                                                                                 <  2:373452/251‑1 (MQ=255)                                                                  aaaCTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCt                                                                                                                                                                      >  1:371564/1‑251 (MQ=255)                                                                  aaaCTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGGCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCTCGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGcc                                                                                                                                                                       <  2:32600/250‑1 (MQ=255)                                                                   aactgggtttGATTGTGGAAACCGAACTAGGGCGTAACCTGATATCGGCGAGACAGAACCGTCACATGCCTCTGGTGACCGCTATCGCGCGAATCTCTGCCCGCTCGGCCGCAGGTTATTCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTTATTGCCGCCCCA‑AATGAAGAAGATGGTGCACGGTTTGTGGCGCTGGGCTCAAAAAAGAACAAGGTCACCTTTTCCGGTAGCCTg                                                                                                                                                                     <  2:269578/243‑1 (MQ=255)                                                                   aaCTGGTGTTGATTATGGAAACCGAAATATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTg                                                                                                                                                                     >  1:7568/1‑251 (MQ=255)                                                                             gATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGat                                                                                                                                                            >  2:305821/1‑250 (MQ=255)                                                                                 aTGGAAACCGAACTATGGCCTAACCTGATTGCGGAGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATCTTTc                                                                                                                                                       >  2:169485/1‑251 (MQ=255)                                                                                  tGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCt                                                                                                                                                      >  2:130821/1‑251 (MQ=255)                                                                                        ccGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATAGTGCTGTAAcg                                                                                                                                                <  1:141057/251‑1 (MQ=255)                                                                                          gAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAAcgcc                                                                                                                                              <  1:311668/251‑1 (MQ=255)                                                                                                    cTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGgct                                                                                                                                    >  1:51035/1‑251 (MQ=255)                                                                                                      aaCCTGATTGCGGCGCTAAATAAACGCAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGACAAACTGGGCAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCACAAAATAACCAGGTAACCGTTCCCGGTAGCCTGAATTTCGCTATTTCTGTAACGCCGCAGTTGgctgc                                                                                                                                  >  2:289548/1‑251 (MQ=255)                                                                                                       aCCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGgctgct                                                                                                                                 >  1:163872/1‑251 (MQ=255)                                                                                                            aTTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGc                                                                                                                            <  2:279318/251‑1 (MQ=255)                                                                                                                 ggcgCTACCTAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGCCTCTCTGCCCGCTCGGCAGCAGGTTATGACAAACTGGGTAAATTCGTCCGTAGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGATACCGGTAGCCTGAAATTCGCTATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTg                                                                                                                        >  1:182848/1‑250 (MQ=255)                                                                                                                       aCATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCt                                                                                                                   <  2:7568/249‑1 (MQ=255)                                                                                                                                  aaaTTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTaa                                                                                                                               >  2:153065/1‑226 (MQ=255)                                                                                                                                  aaaTTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTaa                                                                                                                               <  1:153065/226‑1 (MQ=255)                                                                                                                                  aaaTTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTgg                                                                                                      >  1:190323/1‑251 (MQ=255)                                                                                                                                    aTTCCGCTGGAGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑ACTGAAGCAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTggg                                                                                                     >  2:161967/1‑250 (MQ=255)                                                                                                                                                  cGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCgg                                                                                      >  1:40555/1‑251 (MQ=255)                                                                                                                                                   gCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGt                                                                                     >  2:122541/1‑251 (MQ=255)                                                                                                                                                      aaCGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATg                                                                                  >  2:122616/1‑251 (MQ=255)                                                                                                                                                        cgcgcgACTCTCTGCCCGCTCGGCCGCCGGTTCTGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGAGACGGATTACGCTGATTGCCGCCCAA‑AATGAAGAAGCTGGTGCACGTTTTGTGGCGCTGGGCGCAAAAACTACCCAGGTAACCGTTACCGGTATCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCAGCCAGTGGGCCCCAACCCGCCCGGTATgg                                                                                 >  2:430922/1‑250 (MQ=255)                                                                                                                                                             gACTCTCTGCACGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGC‑CAACAATGAAGAAGCTGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGcc                                                                           >  2:167191/1‑251 (MQ=255)                                                                                                                                                                             ggCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGGATGGATTGCCACCAGCACTCACGatg                                                           >  1:152495/1‑249 (MQ=255)                                                                                                                                                                               ccGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGc                                                         >  1:148201/1‑251 (MQ=255)                                                                                                                                                                               ccGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCGGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGc                                                         <  2:371564/251‑1 (MQ=255)                                                                                                                                                                                        aTGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAgt                                                <  1:122541/251‑1 (MQ=255)                                                                                                                                                                                           ccAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCCGCACTCCCGAAGGCGAagggggggtg                                             >  2:416733/1‑243 (MQ=255)                                                                                                                                                                                                  gggTAAAGTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTTGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTTACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATc                                       <  1:394745/250‑1 (MQ=255)                                                                                                                                                                                                       aaTTCGTCCGTCGCTTGCTGCGTCGTATTACGCTTATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCACTGGGCGCAAAAAATAACCAGGTACCCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTTGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAATGCGAAGAGAGTGTTGTGATCGCCGca                                 <  2:273994/251‑1 (MQ=255)                                                                                                                                                                                                       aaTTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcgccgc                                  >  1:430504/1‑250 (MQ=255)                                                                                                                                                                                                                  cGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAA‑AATGACGAAGATGGTGCACGTTTTGTGGCGCTAGGCGCAAGAAATACCCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCCGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGCAGAGAGTGTGGTGATCGCCACACATCAAGCAtt                      >  1:404114/1‑251 (MQ=255)                                                                                                                                                                                                                       gctgcGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTa                  <  1:223410/250‑1 (MQ=255)                                                                                                                                                                                                                          gcGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGCAGGCGAAGAGAGTGTGGTGATCGCCGCCCATCAGGCATTGTTACAGc              >  2:352488/1‑251 (MQ=255)                                                                                                                                                                                                                          gcGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACa                <  1:182821/249‑1 (MQ=255)                                                                                                                                                                                                                          gcGTCGTATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTGTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTCCCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAGGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGc              <  2:173414/251‑1 (MQ=255)                                                                                                                                                                                                                            gtcgtATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCCCATCAGGCATTGTTACAGCa             >  2:17185/1‑250 (MQ=255)                                                                                                                                                                                                                             tcgtATTCCGCTGCTTGCCGCCCAA‑AATGACGCAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCATGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAAt           >  2:420499/1‑251 (MQ=255)                                                                                                                                                                                                                              cgtATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAAtt          <  2:124556/251‑1 (MQ=255)                                                                                                                                                                                                                              cgtATTACGCTGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGGAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAAtt          <  2:121099/251‑1 (MQ=255)                                                                                                                                                                                                                                 aTTACGCGGATTGCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTccc       <  1:17185/251‑1 (MQ=255)                                                                                                                                                                                                                                      gCTGATTTCCGCCCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAACCAGGTAACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCTCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAAtt  <  1:262609/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                AGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAA‑AATGAAGAAGATGGTGCACGTTTTGTGGCGCTGGGCGCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATT  >  NC_000913/3808844‑3809323

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