Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	10,696	C→T	A221T (GCC→ACC)	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,696	0	C	T	98.2%	182.8 / NA	56	A221T (GCC→ACC)	yaaW	putative enzyme‑specific chaperone YaaW
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (29/26);  minor base G (0/1);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTG  >  NC_000913/10456‑10942                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       gCCCAGCCGTGCCGGATGAGCCAACTTATGGTTGTCCATAATTCCTCAAAAATCATCATCGTATGAATGGTGAAATAATTTCCCTGAATAACGTTAG‑GTTGTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGGCTACGATCTTTTGAGGGGAAAATGAAAACTTTCCCAGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGATTCGTc                                                                                                                                                                                                                                              <  2:517912/251‑1 (MQ=255)                    gCCAACTTAGTGTTGCCCCTAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGc                                                                                                                                                                                                                           >  1:187718/1‑251 (MQ=255)                         cTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGCGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGAGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGAcc                                                                                                                                                                                                                      >  1:370656/1‑251 (MQ=255)                               gttcgccATAATACCTCAAAAAGCATCATCGAATGAATGGGGAAATAATTTCCCGGAATAACTGTGGTGTT‑TGCAGGGCGCGGCATAAGAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCTGTTGCCGGTATCAGATCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCTGTACGGTTGGAATCGTCACGCGATAGGCGCTGTCGCTGACCGCTTTa                                                                                                                                                                                                                <  2:3665/246‑1 (MQ=255)                                  gCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAAcc                                                                                                                                                                                                             <  1:435949/251‑1 (MQ=255)                                   cccATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAAccc                                                                                                                                                                                                            >  2:306533/1‑251 (MQ=255)                                    ccATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCGCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAAcccc                                                                                                                                                                                                           <  2:445996/251‑1 (MQ=255)                                     cATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCa                                                                                                                                                                                                          <  2:139354/251‑1 (MQ=255)                                      aTAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAt                                                                                                                                                                                                         >  2:460044/1‑251 (MQ=255)                                         aTTCCTCAAAAATCATCATCGCATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAttt                                                                                                                                                                                                       >  2:390171/1‑250 (MQ=255)                                         aTTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGCATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGCTCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCAttt                                                                                                                                                                                                       >  1:544687/1‑250 (MQ=255)                                               cAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCg                                                                                                                                                                                                >  1:127710/1‑251 (MQ=255)                                               cAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCg                                                                                                                                                                                                >  1:51851/1‑251 (MQ=255)                                                 aaaaTCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGAt                                                                                                                                                                                                                                       >  2:228007/1‑210 (MQ=255)                                                 aaaaTCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGAt                                                                                                                                                                                                                                       <  1:228007/210‑1 (MQ=255)                                                         catcGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAg                                                                                                                                                                                       <  1:360590/250‑1 (MQ=255)                                                            cGAATGAATGGTGAAATAATTTCCCTGAATAACGGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCTCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCCGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGcc                                                                                                                                                                                   <  2:169829/251‑1 (MQ=255)                                                             gaatgaatGGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCt                                                                                                                                                                                  <  1:507333/251‑1 (MQ=255)                                                               atgaatgGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTGAGTGCCGCACCGACAGGGccgcc                                                                                                                                                                                <  2:579011/251‑4 (MQ=255)                                                                     ggTGAAATAATGTCCCTGAATAACTGTAGTGTT‑ATCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCGTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCTCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGcc                                                                                                                                                                           <  1:266767/250‑1 (MQ=255)                                                                           ataatTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCccgcgcc                                                                                                                                                                    <  1:390171/251‑1 (MQ=255)                                                                                           aCTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGc                                                                                                                                                    <  1:154943/251‑1 (MQ=255)                                                                                                              cgcgGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATtgct                                                                                                                                  >  1:587564/1‑251 (MQ=255)                                                                                                                 ggCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGATTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCg                                                                                                                               >  1:387901/1‑251 (MQ=255)                                                                                                                   cATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTg                                                                                                                             >  1:432540/1‑251 (MQ=255)                                                                                                                          atCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGt                                                                                                                      <  1:151728/251‑1 (MQ=255)                                                                                                                          atCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGATAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCCTTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGt                                                                                                                      >  2:563374/1‑251 (MQ=255)                                                                                                                            cAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTa                                                                                                                     >  1:321718/1‑250 (MQ=255)                                                                                                                                 aGTGGCGCAGAGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTCTCAGACCTGAGTGGAGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTAAGCTCATTGCTGCGTGGATGCGTAAAATGc                                                                                                               >  2:387923/1‑251 (MQ=255)                                                                                                                                   tGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCgg                                                                                                             >  1:11140/1‑251 (MQ=255)                                                                                                                                    ggggCAGTGTCTACGCTCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGACGATTTGCAGTACGGTTGGAATCGTCACGCGATAGCCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGCAAAATCCggg                                                                                                            >  2:457227/1‑251 (MQ=255)                                                                                                                                        cAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTc                                                                                                          >  1:588554/1‑249 (MQ=255)                                                                                                                                        cAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAg                                                                                                        <  2:321718/251‑1 (MQ=255)                                                                                                                                         aGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAAGGCGGGTCAGt                                                                                                       <  2:361829/251‑1 (MQ=255)                                                                                                                                         aGTGACTACTATCTTGTGATGGCAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCTCGAACCAGCCGGCGCAGGCAGGCGATTTTCAGTTCGGGAGGAATCGTCACGCGATAGGCGCTGTCGCTGAGCGTTTTAATCCCATTTAGTGCCTCACCTACGGGGCTTCGCAGTCCCTCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTTCTGCGTGGGTGCGTAAAATGCGGGTCAGt                                                                                                       <  1:336853/251‑1 (MQ=255)                                                                                                                                               tACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTg                                                                                                 >  1:374062/1‑251 (MQ=255)                                                                                                                                                            ggggAAAATGCAAATTTTCCCCGGTTTCCGGCATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAACTGCGGGTCCGTTGGCTGGAAAGCAAATGCg                                                                                    >  2:362773/1‑251 (MQ=255)                                                                                                                                                             gggAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCg                                                                                    >  1:418033/1‑250 (MQ=255)                                                                                                                                                             gggAAAAAGAAAAAATTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCg                                                                                    >  1:515985/1‑250 (MQ=255)                                                                                                                                                               gaaaatgaaaatTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGaca                                                                                 <  1:325538/251‑1 (MQ=255)                                                                                                                                                                       aaagtttCCCCTGTTTCCGGGATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCGCAGTGCTGCGTTGGTGCGTAAAATGCTGGTCAGTTGGCTGGAAAGCAAATGCGATACCCCtttt                                                                         <  2:496237/247‑1 (MQ=255)                                                                                                                                                                         atTTTCCCCGGTTTCCGTTATCAGACCTGAGTGGCGCTAACAATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGCAAGCAAATGCGACCCACCTTTTGc                                                                       >  2:44384/1‑251 (MQ=255)                                                                                                                                                                              ccccGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCAGGCTCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGGTGGATGGAAAGCAAATGCGACACACCTTTCGCCaat                                                                   <  2:61565/250‑1 (MQ=255)                                                                                                                                                                               cccGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGCTTTGCCGGACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAACGCAAATGCGACACACCTTTTGCCaata                                                                  >  2:369718/1‑250 (MQ=255)                                                                                                                                                                                                               aaCCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGctc                                   <  2:127710/249‑1 (MQ=255)                                                                                                                                                                                                               aaCCATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAAGCGGCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAAGGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGctc                                   <  2:51851/249‑1 (MQ=255)                                                                                                                                                                                                                  cATCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCa                              <  1:141785/251‑1 (MQ=255)                                                                                                                                                                                                                   aTCCGGCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTcc                               >  2:159953/1‑249 (MQ=255)                                                                                                                                                                                                                       ggCGCAGGCAGGCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCa                         <  1:61589/251‑1 (MQ=255)                                                                                                                                                                                                                           caggcaggCGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAGCCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCa                     <  1:256298/251‑1 (MQ=255)                                                                                                                                                                                                                             ggcaggcGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCAcc                   >  2:150174/1‑251 (MQ=255)                                                                                                                                                                                                                               caggcGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCt                 >  1:474924/1‑251 (MQ=255)                                                                                                                                                                                                                                aggcGATTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTg                >  2:214360/1‑251 (MQ=255)                                                                                                                                                                                                                                     aTTTGCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATc           >  1:337620/1‑251 (MQ=255)                                                                                                                                                                                                                                         gCAGTACGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTGGGCGGGAAAGCAAATGCGACACACCTTTTGCCAAGAATTGGTCTTTCATCATCAGCGGCAGCAGCGCTTCCAGCTCATTCACCCTGGCATCGAcc       <  2:587564/251‑1 (MQ=255)                                                                                                                                                                                                                                              aCGGTTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTg  <  2:241684/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       GCCCAGCGGTGCCGGATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTT‑TTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTG  >  NC_000913/10456‑10942

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