Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	10,718	G→A	S213S (AGC→AGT)	yaaW ←	putative enzyme‑specific chaperone YaaW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	10,718	0	G	A	95.2%	108.3 / ‑1.9	43	S213S (AGC→AGT)	yaaW	putative enzyme‑specific chaperone YaaW
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (24/16);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

CAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTA‑GTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATG  >  NC_000913/10477‑10960                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    cAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTg                                                                                                                                                                                                                                            <  1:522076/251‑1 (MQ=255)      tAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACTGCTGTAATTGTCACGCGATAGGCACTGCCGCTGACCg                                                                                                                                                                                                                                        <  1:211911/250‑1 (MQ=255)      tAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCg                                                                                                                                                                                                                                        <  1:509807/250‑1 (MQ=255)                    aTTCCGCAAAAATCATCATCGAAGGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTa                                                                                                                                                                                                                         <  2:342067/251‑1 (MQ=255)                          cAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTA‑GTGc                                                                                                                                                                                                                    <  2:288273/249‑1 (MQ=255)                          cAAAAAGCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGctg                                                                                                                                                                                                                  <  1:238047/251‑3 (MQ=255)                            aaaaTCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGctgcg                                                                                                                                                                                                                >  1:611597/1‑247 (MQ=255)                                 catcatcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTA‑GTGCCGCACCTAc                                                                                                                                                                                                           <  1:307868/251‑1 (MQ=255)                                  atcatcGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGctgcgcccac                                                                                                                                                                                                           >  1:400146/1‑241 (MQ=255)                                                   gAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCAATGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACccgc                                                                                                                                                                                         >  1:671098/1‑251 (MQ=255)                                                         atTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATCATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGACTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGATTTAACCCCATTTA‑ATTCTGCGCCCACAGGGCCTCACAGACccgcgccgcg                                                                                                                                                                                   >  2:594424/1‑251 (MQ=255)                                                                     aacagTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCTCGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCTCGCAGCAAACCAt                                                                                                                                                                        <  1:544018/246‑1 (MQ=255)                                                                         gTAGTGTTTTCAGGGCGCGGCCTAATAATCAGCCAGTGGGGCAGTGTCTCCGCTCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGAACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCCTTTA‑ATGCTGCGCCCACAGAGCCTCCCAGACCCGACCCGCGCAGCAAACCATGcacc                                                                                                                                                                   >  2:573100/1‑248 (MQ=255)                                                                         gTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGccc                                                                                                                                                                    >  1:428882/1‑250 (MQ=255)                                                                         gTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGccc                                                                                                                                                                    >  1:543978/1‑250 (MQ=255)                                                                                     gggCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCCAGCACGCTCAtt                                                                                                                                                       >  1:175103/1‑251 (MQ=255)                                                                                                   aatCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCg                                                                                                                                         >  1:340001/1‑251 (MQ=255)                                                                                                    atCAGCCAGTGGGGCAGTGTCTTCGAGCTTTGGAGGGGAAAATGAAAATTTGCCCCGGGTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGTGCTTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACTCTCATTGCTGCGTGGGTGCGt                                                                                                                                        <  2:543978/251‑1 (MQ=255)                                                                                                      cAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCCCTGCCGCTGACCGCTTTAACCCCATTTA‑TTGCTGCGCACACAGGGCCTCCCAGaccc                                                                                                                                                                                           >  2:13605/1‑198 (MQ=255)                                                                                                      cAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGaccc                                                                                                                                                                                           <  1:13605/198‑1 (MQ=255)                                                                                                        gCCAGTGGGGCAGTGTCTACGATCTTGTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGTTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTaaaa                                                                                                                                    <  2:479303/251‑1 (MQ=255)                                                                                                         ccAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAACTTTTCCCCGGTTTACGGTGTCAGCCCTGAGTGGAGCTAACCATCCGGCGCAGGCAGGCGCTTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCACATTTA‑ATGCTGCGCCCACAGGGCATCCCAGACCCCCTCCGCGCAGCAAACCATGCCCAAGCCCGCTCATTGCTACGAGGGGGCGTaatat                                                                                                                                   >  2:594745/1‑248 (MQ=255)                                                                                                           aGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTTCTGCGTGGCGGCGTAAAATGc                                                                                                                                 >  2:92522/1‑251 (MQ=255)                                                                                                                       tCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGGGGGTGCGGAAAATGCGGGTCAATTGGc                                                                                                                     >  2:86445/1‑251 (MQ=255)                                                                                                                            gATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCCCGCTCATTGATGCg                                                                                                                                                 >  2:233777/1‑218 (MQ=255)                                                                                                                            gATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCg                                                                                                                                                 <  1:233777/218‑1 (MQ=255)                                                                                                                                              tgaaaatTTTCCCCGGTTTACGGTGTCCGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGAACGCTTTAACCCCATTTA‑ATGCTGCGCCCCCAGGGCCTCCCCGACCCGAGCCGCGCAGCAAACCAAGCCCCAGCACGCTCATTGCTGCGTGGGt                                                                                                                                            >  2:27616/1‑205 (MQ=255)                                                                                                                                              tgaaaatTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGt                                                                                                                                            <  1:27616/205‑1 (MQ=255)                                                                                                                                                        ccccGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCt                                                                                                                                                                                                    >  1:378857/1‑139 (MQ=255)                                                                                                                                                           cGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGACAGGCACTGCCGCTGGCCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATACCCAAGCACGCCAATTGCCGCGTGGGAGCGGCAACTGCGGGTCAATTGCCTGGAAAGCAAATGCGACACCCCtttacccatttatt                                                                                 >  2:580999/1‑240 (MQ=255)                                                                                                                                                            ggTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGCCACGCCTTTTGCCAATAAttt                                                                                >  1:477370/1‑251 (MQ=255)                                                                                                                                                                                aGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGgcc                                                            >  1:421635/1‑250 (MQ=255)                                                                                                                                                                                       cTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCTCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCAAGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCcgcagcgc                                                     >  2:303943/1‑249 (MQ=255)                                                                                                                                                                                        tAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCAGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGAGACCCGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGctct                                                    >  2:169097/1‑251 (MQ=255)                                                                                                                                                                                               ccGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCt                                             >  1:562704/1‑251 (MQ=255)                                                                                                                                                                                                 ggCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTc                                            >  1:242494/1‑250 (MQ=255)                                                                                                                                                                                                               aTTTGAAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCGGACCGCTTTAACCCCATTTA‑ATGCTTCGCCCACCGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAATCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCAt                              <  2:200148/250‑1 (MQ=255)                                                                                                                                                                                                                   caaGTACGGCTGGAAGTGTCACGCGATAGGCACGGCCGCTGACCGCTTTAACCCCATCTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCTCGCTCAGTGCTGCGTGGGTGCGTAAAATGCGTGTCAAGGGGCTGGAAAGCAAATGCGACACGCCGTTTGCCAATAATTTGTCTTTCATCAGCCGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGAATCGAc                          <  2:242494/248‑1 (MQ=255)                                                                                                                                                                                                                          ggCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCa                  <  1:93534/251‑1 (MQ=255)                                                                                                                                                                                                                                 attgTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTc           >  1:654241/4‑251 (MQ=255)                                                                                                                                                                                                                                    tctccgcgATAGGCACTACCTCTAACCGCTTTAACCACATTTACATCTCGC‑CCTCCAGGGCCGCCCATGCCCGCGCCTCACACAAAACCATGCCCAAGCACGCTCATTGCTGCGTGGTTGCTTAAAATGCGGGTCAATTGGCTTGAAAGCAAATGCGACACGCCTTTTGCCAACTATTTGTCTTCCATCAGCAGCGGCAGCAGCTCGTCCAGCGTAGTCACCCTGGCATCGACCGCGTGCAGAAAcgcca         <  2:344785/246‑5 (MQ=255)                                                                                                                                                                                                                                      cACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGc       <  1:92522/250‑1 (MQ=255)                                                                                                                                                                                                                                          cgATAGGCACTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCCGCTCTTCCAGCTCATTCACCATGGCATCGACCGCGTGCAGAAACTCCTGCTTATg  >  2:9616/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    CAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTA‑GTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATG  >  NC_000913/10477‑10960

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