Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,999,664	C→T	S561L (TCG→TTG)	uvrD →	DNA helicase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,999,664	0	C	T	93.5%	118.8 / 1.0	47	S561L (TCG→TTG)	uvrD	DNA helicase II
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (20/23);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

ACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACGGCAACGCGCCAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTCGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCCGCGCGATGCAGAAACTGACGCTGACCTACGCGGAAACCCGCCGTCTGTATGGTAAAGAGGTTTACCATCGCCCGTCGCGCTTTATCGGCGAGCTGCCGGAAGAGTG  >  NC_000913/3999423‑3999901                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                              aCTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACGGCAACGCGCCAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAgg                                                                                                                                                                                                                                      <  1:100594/251‑1 (MQ=255)        gggTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACGGCAACGCGCCAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGAGGTGAAGTTGATGACGCTACACTTGGCGAAAGGCCTgg                                                                                                                                                                                                                                 >  1:193303/1‑249 (MQ=255)                  aGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACGGCAACGCGCCAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAgg                                                                                                                                                                                                                     <  1:326951/251‑1 (MQ=255)                  aGACTCCGGCATGCGTACCATGTATGAGCAGGAGAAGGACGAAACAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACGGCAACGCGCCAGTTCCGCTACAACGGAGAAGGCGAAGATTTAATGCCGCTGCAGGCGTTCATCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGGGCTGTTGATGACGCTACACTTGGCGAAAGGGCTGGGGTTCCCGCa                                                                                                                                                                                                                       >  2:209982/1‑249 (MQ=255)                   gACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGTCGCGTATCGCAAACTTAGAGGAACTGGTGACGGCAACGCGCCAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGgt                                                                                                                                                                                                                    >  1:228333/1‑251 (MQ=255)                       ccGGACTGCGTACCATGTATGAGGAGTAGAAGGGCGAAAAAGGTCAGACGAGTAGCGAAAACTTAGAGGAACAGGGGACGGCAACGCGCCAGTGCTGCTTCAACGAAGAAGACGAAGATTTAATGCAGCTGCAGGCGTTCCTAGCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTTCACTTGGCGAACGTCCTGGAGTTCCCGCAGgtgt                                                                                                                                                                                                                  <  1:348955/249‑1 (MQ=255)                         ggacTGCGTACCAAGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGGGACGGCAACGCGCCAGTTCAGCTACAACGAAGAAGACGAAGATGTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTAt                                                                                                                                                                                                              <  2:22008/248‑1 (MQ=255)                          gCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACGGCAACGCGCCAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTGTATc                                                                                                                                                                                                             <  1:257889/251‑1 (MQ=255)                              gCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACGGCAACGCGCCAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATCCCTGGCAGGATGCGGTGCAGTGGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGt                                                                                                                                                                                                           <  1:123707/249‑1 (MQ=255)                                gTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACGGCAACGCGCCAGTCCAGCTACAACGAAGAAGACGAAGATTTAAGGCCGCGGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTGGATGACGCTACACTTGGCGAAAGGCCTTGAGTTCCCGCAGGTGTTTATCGTTGGt                                                                                                                                                                                                       <  2:316580/251‑1 (MQ=255)                                           gAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACGGCAACGCGCCAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAgg                                                                                                                                                                                             <  1:209982/250‑1 (MQ=255)                                                  agaAGGGCGAAAAAGGTCAGACGCGTATAGAAAACTTAGAGGAACTGGTGACGGCAACGCGCCAGTTCAGCTACAAAGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTACTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGCAGCGGTGCAGTTGATGAAGCTACAATTGGAGAAAGGCCTGGAGTTCCCGCAGGAGATTATCGTTGGTATGGAAGAGGACATggtc                                                                                                                                                                                     >  1:191517/1‑248 (MQ=255)                                                          gAAAAAGGTCAGAAGTGTATAGAAAACTTAGAGGAACTGGTGACGGCAACGTGCCAGGTCAGCTACAACGAAGAAGACGAAGATGTAAGGCCGCTGCATGCGTTCCTCTCCCAGGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCGGGCAGGATGCGGTGCAGTGGATGACGTTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATTGAAGAGGGCATGTGCTCAAGccg                                                                                                                                                                             <  2:311378/251‑2 (MQ=255)                                                                         cgtacagAAAACTTAGAGGAACTGGTTACGGCAACGCGCCAGTTCAGCTACAACGAAGAAGACGAATATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTAGCGTGGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGa                                                                                                                                                              <  2:239675/245‑1 (MQ=255)                                                                              cGAAAACTTAGAGGAACTGGTGACGGCAACGCGCCAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAggcg                                                                                                                                                         <  1:320463/251‑1 (MQ=255)                                                                                aaaaCTTAGAGGAACTGGTGACGGCAACGCGCCAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAggcggg                                                                                                                                                       >  1:218514/1‑251 (MQ=255)                                                                                  aaCTTGGAGGAACCGGTGACGGCAACGCGCCAGTTCATCTACAACGAAGAAGACGAATATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAggcgggcg                                                                                                                                                     <  1:244002/251‑1 (MQ=255)                                                                                    cTTAGAGGAACTGGTGACGGCAACGAGCCAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGCGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTc                                                                                                                                                   >  1:111486/1‑251 (MQ=255)                                                                                            aactgtggtctgacactcgcCAGTTCAGTTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGGTTCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCGGTCTGATACCTGGCAGGATGCGGTTCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTAGCGTTTGTATGGAAGAGGGCATGTTCCAAAGCCAGATGTCGCTGGATGAAGGCGTGCGTCTGgaagaa                                                                                                                                           <  2:283513/234‑1 (MQ=255)                                                                                              cTGGTGACGGCAACGCGCCAGTTCAGCTACAACGAAGACGAAGAAGATATAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGCCAGGATGCGGTGCAGTTGATGCAGCAACACTTGGCGAAAGGACTGGTGTTCCAGCGGGTGTATATCGATGGCATGGAAGAGGGCATCTTCCCAAGCCAGATGTTGCTGGACAAAggcgagggttgggaagtaa                                                                                                                                          >  2:208382/1‑235 (MQ=255)                                                                                                   gACGGCAACGCGCCAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGcc                                                                                                                                    <  1:344704/251‑1 (MQ=255)                                                                                                           cgcgCCAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTCTCGTTGGTATGGTAGAGGGCATGTTCCCAAGCCAGCTGTCGCTGGATGAAGGTGGGCGTCTGGAAGAAGAAGCCCGTCTGGcc                                                                                                                            >  2:316531/1‑251 (MQ=255)                                                                                                           cgcgCCAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGcc                                                                                                                            >  1:14134/1‑251 (MQ=255)                                                                                                                cAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGTTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTCGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGGCGTCTGGCCTACGt                                                                                                                       >  1:77107/1‑251 (MQ=255)                                                                                                                              cGAAGAAGACGAAGATTTAATGCCGCTGCCGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCcgcg                                                                                                         >  1:268907/1‑251 (MQ=255)                                                                                                                                aagaagACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCcgcgcg                                                                                                       >  1:5202/1‑251 (MQ=255)                                                                                                                                                  tGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGAGCGGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCCGGTGTTTATCGTTGGTATGGAAGAGGGAATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCCGCGCGATGCAGAAACTGACGCTg                                                                                     >  1:147397/1‑251 (MQ=255)                                                                                                                                                   gCCGCTGCAGGCGTGCCGCTCCCACGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACTGGGCAGGAGGCGGTGCAGATGATGACGCTACACTGGGCGAAATGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTGCCCAAGCCATATGTCGCTTGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTGGTGCCTTCGTTGGCGTAACCCGCTCGATGGAGAAACTGACTCTg                                                                                     <  2:54832/250‑1 (MQ=255)                                                                                                                                                       cTGCAGGCGGTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTAGCGTTGGTAGGGAAGAGGGCCTGTTCCCAAGCCAGATGTCGCTGGGTTAAGGCGGGCGGCTGGGAGTAGAACGCTGTCTGGCCAACGTTGGCGTAACCCGCCGGAAGCAGAAACTTGCGCTGACCt                                                                                 >  2:19181/1‑250 (MQ=255)                                                                                                                                                             gCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCCGCGCGATGCAGAAACTGACGCTGACCTACGCGGa                                                                          <  1:197869/251‑1 (MQ=255)                                                                                                                                                               gTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCCGCGCGATGCAGAAACTGACGCTGACCTACGCGGaaa                                                                        <  1:244551/251‑1 (MQ=255)                                                                                                                                                                 tcctccCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGTCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCaa                                                                                                                                                                                 <  2:156596/144‑1 (MQ=255)                                                                                                                                                                 tCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCaa                                                                                                                                                                                 >  1:156596/1‑144 (MQ=255)                                                                                                                                                                  ccTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCCGCGCGATGCAGAAACTGACGCTGACCTACGCGGAAAccc                                                                     <  1:78906/251‑1 (MQ=255)                                                                                                                                                                  ccTCCCCCATGCGGCACTGGAAGTAGGTGAAGGGCAGGCGGATTCCTGGCCGGAGTCGGTTGAGGTGATGCCGCTACACTTGGCGAAAGGTCTGGAGTTCCCGCAGGTGTTTATCGGTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCCGCTCGATGCAGAAACTTCCGCTGACCTACGCGGAAAcc                                                                      <  2:193303/250‑1 (MQ=255)                                                                                                                                                                     ctcCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCCGCGCGATGCAGAAACTGACGCTGACCTACGCGGAAACccgcc                                                                  >  1:147483/1‑251 (MQ=255)                                                                                                                                                                         cATGCGGCCCTGGAAGCAGTTGAAGGGCAGGCGGATACCTGGGAGGATGCTGTGCAGTTGATGACGTTACACGCGGAGAAAGGCCTGGAGTGCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGACAGATGTTGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGGCGTCTGGCCTACGTTGGCGTAACCCGCGCGATGCAGAAACTGACGCTGACCTACGCGGAAACCCGCCGTCt                                                              <  2:77107/251‑1 (MQ=255)                                                                                                                                                                                   tGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCCGCGCGATGCAGAAACTGACGCTGACCTACGCGGAAACCCGCCGTCTGTATGGTaa                                                     >  1:331595/1‑250 (MQ=255)                                                                                                                                                                                    ggAAGCAGGTGAATGGCAGGCGGATACCGGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGGCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGAGGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCCGCGCGATGCAGAAACTGACGCTGACCTACGCGGAAACCCGCCGTCTGTATGGTAAag                                                   <  2:129700/251‑1 (MQ=255)                                                                                                                                                                                        gCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTCGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCCGCGCGATGCAGAAACTGACGCTGACCTACGCAGAAACCCGCCGTCTGTATGGTAAAGAGGt                                               >  1:169504/1‑251 (MQ=255)                                                                                                                                                                                                          gATACCTGGCAAGATGCGACGCAGTTGATGACGCTACACATGGCGAAATGCCGGGAGTTCCCGCAGGCGGTTATAGTTGGTATGGAAGACGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAATAACGCCGTCTGGCCTACGTTGGCGTAACCCGGGAGAAGAAGACACTGAAGCTGACCTAAGCGCGAACCCGCCGGCCGTATGGAAAAGGGGTTTCACACCGCCAGTcgcg                             >  2:200971/1‑251 (MQ=255)                                                                                                                                                                                                             aCCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCCGCGCGATGCAGAAACTGACGCTGACCTACGCGGAAACCCGCCGTCTGTATGGTAAAGAGGTTTACCATCGCCCGTCGCGCtt                          <  1:315112/251‑1 (MQ=255)                                                                                                                                                                                                                   cAGGATGCGGTGCATTTGATGACGCTTGACTTGGCGAAATGCCTGGAGTTCCCTCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCGGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCCTCGCGATGCAGAAACTGACGCTTACCTACGCGGAAACCCGCCGTCTGTATTGTAAAGAGGTTTACCATTGACCGTCGCGCTTTATCgg                    <  2:111486/251‑1 (MQ=255)                                                                                                                                                                                                                                aGTTGATGACGCTACACTTGGCGAAAGGCCTGGAGTTCCAGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGCAGAACGCCGTCTGGCCTACGTTGGCGTAACCCGCGCGATGCAGAAACTGACGCTGACCTACGCGGAAACCCGCCGTCTGTATGGTAAAGAGGTTTACCATCGCCCGTCGCGCTTTATCGGCGAGCTGCCGGa        >  2:197300/1‑250 (MQ=255)                                                                                                                                                                                                                                    gatgaCGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCCGCGCGATGCAGAAACTGACGCTGACCTACGCGGAAACCCGCCGTCTGTATGGTAAAGAGGTTTACCATCGCCCGTCGCGCTTTATCGGCGAGCTGCCGGAagag    <  1:256979/250‑1 (MQ=255)                                                                                                                                                                                                                                    gatgaCGATACACTTGGCAAAAGGCCTGGAGTTCACGCAGTTGTTTATCGTTGGTATGGAAGCGGGCATGTTCACAAGCCCGATGTCGCTGGCTGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTACCCCGCGCGATGCAGAAACTGACGATGACCTACGCGGAATCACGCCGTCTGTATAGAAAAGCGGTTTACAATCGCACGACGCGATTTATCGGCAAGCTACCgggagggt   >  2:211583/1‑245 (MQ=255)                                                                                                                                                                                                                                     atgaCGCTACACTTGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCCGCGCGATGCAGAAACTGACGCTGACCTACGCGGAAACCCGCCGTCTGTATGGGAAAGAGGTTTACCATCGCCCGTCGCGCTTTATCGGCGAGCTGCCGGAAGAgtg  >  1:207154/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                              ACTGACCGGGTAATTAAAGACTCCGGCCTGCGTACCATGTATGAGCAGGAGAAGGGCGAAAAAGGTCAGACGCGTATCGAAAACTTAGAGGAACTGGTGACGGCAACGCGCCAGTTCAGCTACAACGAAGAAGACGAAGATTTAATGCCGCTGCAGGCGTTCCTCTCCCATGCGGCACTGGAAGCAGGTGAAGGGCAGGCGGATACCTGGCAGGATGCGGTGCAGTTGATGACGCTACACTCGGCGAAAGGCCTGGAGTTCCCGCAGGTGTTTATCGTTGGTATGGAAGAGGGCATGTTCCCAAGCCAGATGTCGCTGGATGAAGGCGGGCGTCTGGAAGAAGAACGCCGTCTGGCCTACGTTGGCGTAACCCGCGCGATGCAGAAACTGACGCTGACCTACGCGGAAACCCGCCGTCTGTATGGTAAAGAGGTTTACCATCGCCCGTCGCGCTTTATCGGCGAGCTGCCGGAAGAGTG  >  NC_000913/3999423‑3999901

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