Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	10,973	T→A	K128N (AAA→AAT)  I48I (ATT→ATA)	yaaW ← mbiA →	putative enzyme‑specific chaperone YaaW uncharacterized protein MbiA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	10,973	0	T	A	95.3%	123.4 / ‑0.8	43	K128N (AAA→AAT)  I48I (ATT→ATA)	yaaW mbiA	putative enzyme‑specific chaperone YaaW uncharacterized protein MbiA
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (20/21);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.42e-01

AGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATT  >  NC_000913/10746‑11208                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            aatgctgcgcccaCAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATAttcttc                                                                                                                                                                                                                                       <  2:277363/222‑1 (MQ=255) aatgctgcgcccaCAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATAttcttc                                                                                                                                                                                                                                       >  1:277363/13‑234 (MQ=255)              cAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCAGGCCCAAGCACGCTCAGTGCTGCGGGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAAc                                                                                                                                                                                                         <  2:259088/251‑1 (MQ=255)               aGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACt                                                                                                                                                                                                        >  2:148810/1‑251 (MQ=255)                 ggCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCCCCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTAAGTTCCTCGTCCATATTCTTCCAGTTCTTACGCAGAAATTGTTCCAGTAACTGt                                                                                                                                                                                                      >  2:337739/1‑251 (MQ=255)                   ccTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTg                                                                                                                                                                                                    <  1:117232/251‑1 (MQ=255)                                ggtcgcgcAGCACACCATGCCCAAGCACGCTCATTGCTGCGTGGGTTAGTAACAGGCGGGTCAATTGGCTGGAAAGCAATTGCGACACGCCTTTTGCCAATAATGTGTCTTACATCAGCAGCTGCAGCAGCTCTTCCAGCTCATGCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAAc                                                                                                                                                                                       <  2:396767/248‑1 (MQ=255)                                   cgcgcAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTACCTGTTGCTCAATCTCAAACGt                                                                                                                                                                                     >  1:383937/1‑250 (MQ=255)                                   cgcgcAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGt                                                                                                                                                                                     >  2:74714/1‑250 (MQ=255)                                   cgcgcAGCAAACCAAGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTTCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACCCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTa                                                                                                                                                                                    <  1:433007/251‑1 (MQ=255)                                    gcgcAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAg                                                                                                                                                                                   <  1:34738/251‑1 (MQ=255)                                      gcagcaAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGa                                                                                                                                                                                  >  2:385802/1‑250 (MQ=255)                                                   cccAAGCACGCTCAGTGCTGCGTGGGGGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTACCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCgg                                                                                                                                                                    <  2:163184/251‑1 (MQ=255)                                                       agcaCGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCttt                                                                                                                                                                >  1:398940/4‑251 (MQ=255)                                                           cGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCATCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGCCTTTCAGc                                                                                                                                                            >  2:106425/1‑251 (MQ=255)                                                              tCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTc                                                                                                                                                         >  1:406614/1‑251 (MQ=255)                                                                    ctgcGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGtt                                                                                                                                                                                                                 <  2:52209/189‑1 (MQ=255)                                                                    ctgcGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGtt                                                                                                                                                                                                                 >  1:52209/1‑189 (MQ=255)                                                                          gggTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGa                                                                                                                                              >  2:416554/1‑250 (MQ=255)                                                                               cGTAAAATGGGAGTCAATCGGTTGGAAAGCAAATGCGGCACGCCTTTTGCCAATACTTTGTCTTTAATAAGCAGCGGCAGCGGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCCGAAACTCCAGCTTCTGTTCCTCGTCCATACTCTTCCAGGTGTTACGCAGCCATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTATTTGTCGGCTTTCAGCTTCAATAGTTTTGAAACATc                                                                                                                                        >  2:207401/1‑251 (MQ=255)                                                                               cGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTGCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATc                                                                                                                                        <  2:377396/251‑1 (MQ=255)                                                                                gTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCACTCTCAAACGTAGACATCTCTTTGTCGGCTTTCa                                                                                                                                                              >  2:314799/1‑228 (MQ=255)                                                                                gTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCa                                                                                                                                                              <  1:314799/228‑1 (MQ=255)                                                                                     aTGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCaa                                                                                                                                  <  2:502995/251‑1 (MQ=255)                                                                                     aTGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTACAGTAACTGTTGCTCAATCTCAACCGTAGCCATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGCGCa                                                                                                                                   >  2:200077/1‑250 (MQ=255)                                                                                                                      acacCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGc                                                                                                                                                                   <  2:181467/185‑1 (MQ=255)                                                                                                                      acacCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGc                                                                                                                                                                   >  1:181467/1‑185 (MQ=255)                                                                                                                         ccTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGttt                                                                                              >  1:338026/1‑251 (MQ=255)                                                                                                                            tttGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTttg                                                                                             <  1:148810/249‑1 (MQ=255)                                                                                                                                 cAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCg                                                                                       >  1:327282/1‑250 (MQ=255)                                                                                                                                  aataatTTGGCTTTCATCAGCAGCGGCAGCAGCTCTTCGAGCTCATTCCCCCTGGCATCGCCCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCCGTTTGTTGGCGac                                                                                     <  2:444096/251‑2 (MQ=255)                                                                                                                                          gccatgcctcctcagcGGCTGCAGCTCTTCCAGCTCAGTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCGCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTGCAATCGTTTTGAAACAGCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGGTTGTTGGCGATACTATCg                                                                              <  2:77600/238‑1 (MQ=255)                                                                                                                                               tcatcaGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGccacc                                                                         <  1:385802/250‑1 (MQ=255)                                                                                                                                                catcaGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTACATATTCTTCCAGGTGTTACGCAGAAATTATTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGATTTCAGCTTCACTCGTTTTGAAACATCGAGCAATATGGCCCGATACAATTTACCGTGTCCGCGCTGTTTGTTGGCGATACTATCGccacca                                                                        >  2:108081/1‑250 (MQ=255)                                                                                                                                                                       ccAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATc                                                <  1:181078/251‑1 (MQ=255)                                                                                                                                                                          gCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCTTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTTAAACGTAGACATCTCTTTGTCGGCTGTCAGCTTCAATCGTTTTGAAACAGCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGc                                             <  2:398940/251‑1 (MQ=255)                                                                                                                                                                          gCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGc                                             <  1:50612/251‑1 (MQ=255)                                                                                                                                                                           cTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGACCATATTCATCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCt                                            >  1:351129/1‑251 (MQ=255)                                                                                                                                                                                   cccTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTg                                    <  1:30414/251‑1 (MQ=255)                                                                                                                                                                                                 gcgTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCGCCAGCAATCAGCTGCCAGTTTCGGCGATGTTGCTc                      <  1:410359/251‑1 (MQ=255)                                                                                                                                                                                                      cAGAAACTACGGCTTATGTTCCTCGTCCATATTCTTCCAGGTGCTACTCAGAAATTGTTTCAGTAACTGTTGCTCAATATCAAACGTAGACATCGCGTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCTATACTATCGCCACCAAAATGCTGGAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCGCGGGAt                 <  2:229798/251‑1 (MQ=255)                                                                                                                                                                                                                     aTGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAACATGCTGTAATTCTCCAGCACTCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGAtt  >  2:401600/1‑251 (MQ=255)                                                                                                                                                                                                                     aTGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGAtt  >  1:334339/1‑251 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            AGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATT  >  NC_000913/10746‑11208

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