Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	298,885	G→A	P684S (CCG→TCG)	paoC ←	aldehyde dehydrogenase: molybdenum cofactor‑binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	298,885	0	G	A	82.2%	104.8 / 11.8	45	P684S (CCG→TCG)	paoC	aldehyde dehydrogenase: molybdenum cofactor‑binding subunit
Reads supporting (aligned to +/- strand):  ref base G (5/3);  new base A (25/12);  total (30/15)
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

GAAGCGCAGAATGGTCCGGTGGTCATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCG  >  NC_000913/298637‑299128                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    gAAGCGCAGAATGGCACGGTGGTCATCGGTGAGACACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCgggg                                                                                                                                                                                                                                                   >  1:134042/1‑251 (MQ=255)      gcAGAATGGTCCGGTGGTCATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGat                                                                                                                                                                                                                                              >  1:135229/1‑251 (MQ=255)         gAATGGTCCGGTGGTCATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATg                                                                                                                                                                                                                                           <  2:105954/251‑1 (MQ=255)             ggTCCGGTGGTCAGCGGTGAGAAACGCCTGCTCTGTCTTATGCCAGTGGTCGTCTTTGTCAAAAAGCGGGTATGACAGTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACGCCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCTAGGAGGATATGGGGt                                                                                                                                                                                                                                       <  2:306435/251‑1 (MQ=255)                        cATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTatcatc                                                                                                                                                                                                                             <  2:248711/250‑1 (MQ=255)                                agaAACGCCTGCTCTGGCTTATGCCAGTGTTCGGCTTTGTCAAAAAGAGGGCATCAAATTGTTCCTCCTTAAAACACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGCGACAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATAGAGGAGGATATGGGGTCGGTATCATCCAGGAAAAt                                                                                                                                                                                                                    >  1:303387/1‑251 (MQ=255)                                   aaCGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCAc                                                                                                                                                                                                                 >  1:31218/1‑251 (MQ=255)                                         ggCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTg                                                                                                                                                                                                           <  2:425474/250‑1 (MQ=255)                                             ctGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGtt                                                                                                                                                                                                         <  2:372494/249‑1 (MQ=255)                                               ggCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTg                                                                                                                                                                                                      <  1:34019/250‑1 (MQ=255)                                                   tATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGa                                                                                                                                                                                                   <  1:223372/249‑1 (MQ=255)                                                               cGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAAc                                                                                                                                                                                      <  1:166356/250‑1 (MQ=255)                                                                        cAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTc                                                                                                                                                                             >  1:49565/1‑250 (MQ=255)                                                                         aaaaaGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCAt                                                                                                                                                                           <  1:296900/251‑1 (MQ=255)                                                                         aaaaaGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCAt                                                                                                                                                                           <  2:210477/251‑1 (MQ=255)                                                                          aaaaGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTCTCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCAt                                                                                                                                                                           >  1:235759/1‑250 (MQ=255)                                                                          aaaaGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCAt                                                                                                                                                                           >  1:445102/1‑250 (MQ=255)                                                                                   tATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCAGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGATTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCat                                                                                                                                                                 >  1:417655/1‑251 (MQ=255)                                                                                        cATTGTTCCTCCTTACACCACACCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGGGGCGTTATACAACGCGTTGGCGATAGACGCGCTCACGCCGCACAGGCCCAGCTCACCGACAGCTGTGGCCTTCATCGGGGAGAATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGTGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCgg                                                                                                                                                            >  2:284477/1‑250 (MQ=255)                                                                                                         ccACACCCGGAAGCTTATCGGCCGGCTTATCCAGAGTGATGGGAAAACCCCGTACCCGAATCCCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACCCCTGGGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGTTTTGGGATATCCGCATGACCCGGCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGccc                                                                                                                                           <  1:63982/251‑1 (MQ=255)                                                                                                            caTCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACCCCTTTGGCCTTCATCGAggagga                                                                                                                                                                                                                                               <  2:13019/148‑1 (MQ=255)                                                                                                            caTCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAggagga                                                                                                                                                                                                                                               >  1:13019/1‑148 (MQ=255)                                                                                                             aTCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc                                                                                                                                       >  1:471127/1‑251 (MQ=255)                                                                                                                        tATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCAcc                                                                                                                            >  1:489597/1‑251 (MQ=255)                                                                                                                        tATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGGCATCCAcc                                                                                                                            >  1:124893/1‑251 (MQ=255)                                                                                                                         agcGAGCAGTTTATCCAGAGTGATGGGATACTCCCGTCCCCGAATCCCGGGGGCTTTATACCCCGCGTTGGCGATCGCCGCGCTCACGCCGCACAGGCCCGGCTCACCGACCCCTTTGCCCTTCATCGAGGAGGATAGGGGGTCGTTACCACCCAGAAAAATCCCCCCCTTTTTGGGGATACCCGCATGAACCGCCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCg                                                                                                                           <  1:103117/249‑1 (MQ=255)                                                                                                                              gcagcTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGc                                                                                                                      >  1:360930/1‑251 (MQ=255)                                                                                                                                  cTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcct                                                                                                                  >  1:198548/1‑251 (MQ=255)                                                                                                                                       ccAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATc                                                                                                             >  1:413244/1‑251 (MQ=255)                                                                                                                                               aTGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgc                                                                                                     >  1:371354/1‑251 (MQ=255)                                                                                                                                                               cccGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                                                     >  2:213983/1‑251 (MQ=255)                                                                                                                                                                ccGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                                                     >  1:269007/1‑250 (MQ=255)                                                                                                                                                                    aTACCAGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAACCGGTCATCCACCTCCAGCTACTCCATCAGTGCCCCGCCCATCGCCATAGTCaattc                                                                                >  2:362355/1‑247 (MQ=255)                                                                                                                                                                      cccGGGGGCTTTATACACCGCGTTGGCGATAGCCGCGCTCCCGCCGCACAGGCCCAGCTCACCGACACTTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCa                                                                                                                              <  1:10770/202‑1 (MQ=255)                                                                                                                                                                      aCCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCa                                                                                                                              >  2:10770/1‑203 (MQ=255)                                                                                                                                                                          gTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAAt                                                                          >  1:475958/1‑251 (MQ=255)                                                                                                                                                                               gTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt                                                                     >  1:283969/1‑251 (MQ=255)                                                                                                                                                                                     cacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgc                                                               >  1:267967/1‑251 (MQ=255)                                                                                                                                                                                        cgcgTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcg                                                            >  2:10807/1‑251 (MQ=255)                                                                                                                                                                                              ggCGATAACAGCGCTCCTGCCGCACAAGCCCAGCTCCACGACCCATTTAGCCTTCATCGAGGAGGCAACGGGGTCGGTCTCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCAACTCATACCACGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATACACCGCCAGCTCCTCCATCAGTGCCGCGCCCATCCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt                                                      >  2:465416/1‑251 (MQ=255)                                                                                                                                                                                                          gcTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg                                           >  2:97298/1‑250 (MQ=255)                                                                                                                                                                                                          gcTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg                                           >  1:356053/1‑250 (MQ=255)                                                                                                                                                                                                                gccgcACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTc                                    >  1:409299/1‑251 (MQ=255)                                                                                                                                                                                                                  cccacggccccagccccccgcccccTTTGGCCTTCATCGAGGAGGATATGGGGTCGTTATCCCCCAGAAAACTCCCCTCCTTTTTTGGAATACCCGCATAAACCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt                                  <  1:10807/228‑1 (MQ=255)                                                                                                                                                                                                                                                 cTTCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGAccc   >  1:97527/1‑251 (MQ=255)                                                                                                                                                                                                                                                  ttCATCGAGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGACACCTCATACCACGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCAACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCg  >  1:1054/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    GAAGCGCAGAATGGTCCGGTGGTCATCGGTGAGAAACGCCTGCTCTGGCTTATGCCAGTGTTCGTCTTTGTCAAAAAGCGGGTATGACATTGTTCCTCCTTAAACCACATCCGGCAGCTTATCGAGCAGCTTATCCAGAGTGATGGGATAATCCCGTACCCGAATACCGGTGGCGTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACACAGCGAGCATACGCCGGACCCG  >  NC_000913/298637‑299128

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