Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	10,747	G→A	L204L (CTA→TTA)	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,747	0	G	A	95.7%	126.0 / ‑2.0	48	L204L (CTA→TTA)	yaaW	putative enzyme‑specific chaperone YaaW
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (29/16);  total (30/18)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

CAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTA‑GTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTAC  >  NC_000913/10502‑10989                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   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)                                                                                                                                                                                                                      ggcacTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCCCAGCGCCTCCCAGCCCCGCGCCGAGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCACTTGGCTGGAAAGCAACTGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTAtactccc                          >  2:679355/5‑245 (MQ=255)                                                                                                                                                                                                                          cTGCCGCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACAATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTACTGCTACTGTTCCTCGTc                      >  2:447729/1‑251 (MQ=255)                                                                                                                                                                                                                            gCCGCTGACCGCCTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCa                    <  1:17093/251‑1 (MQ=255)                                                                                                                                                                                                                             ccGCTGACCGCTTTAACGCTAATTA‑ATGCGGCGCCCTCAGGGCCTCCCAGACCCGCGCCTCGCAGCAAACCATGCCCAAGCACGCTCATTGCTTCGTGGGTGCGTAAAATGCGGGTCAATGGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCAt                   <  2:489802/251‑1 (MQ=255)                                                                                                                                                                                                                              cGCTGAACGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCACGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTTCTCGTCCAtt                  >  2:662177/1‑251 (MQ=255)                                                                                                                                                                                                                               gCTGACCGCTTTAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGGAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCAtat                 <  1:575848/251‑3 (MQ=255)                                                                                                                                                                                                                                        tttAACCCCATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAgg        >  1:278852/1‑251 (MQ=255)                                                                                                                                                                                                                                             ccccATTTA‑ATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTa   <  1:594424/251‑1 (MQ=255)                                                                                                                                                                                                                                              cccATTTA‑ATGCTGCGCCCACAGGGCATCCCAGACCCGCGCCGCGCAGCAAACCACGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCCTCGCCCGCGCGCAGACACTCCAGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTAc  >  2:366859/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   CAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTA‑GTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTAC  >  NC_000913/10502‑10989

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: ‑