Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	11,021	T→C	E112E (GAA→GAG)  I64I (ATT→ATC)	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	11,021	0	T	C	95.7%	129.5 / ‑2.1	48	E112E (GAA→GAG)  I64I (ATT→ATC)	yaaW mbiA	putative enzyme‑specific chaperone YaaW uncharacterized protein MbiA
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (25/20);  total (26/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 = 1.44e-01

CGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGG  >  NC_000913/10777‑11267                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       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)                                           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                                                                                                                                                                                                                              >  1:181467/1‑185 (MQ=255)                                                                                       acacCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGc                                                                                                                                                                                                                              <  2:181467/185‑1 (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)                                                                                                                                                                                                            catgtgttACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTGCAGCTTCAATCTTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCAGCGATTTAAACAGTTCGTTGCGTATCat                                       <  1:273054/245‑2 (MQ=255)                                                                                                                                                                                                            cAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGGCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAg                                       <  1:273060/251‑1 (MQ=255)                                                                                                                                                                                                                  ttACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGATCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCt                                 >  1:220465/1‑251 (MQ=255)                                                                                                                                                                                                                    aCGCAGAAATTTTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATTCAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTgg                               <  2:382725/251‑1 (MQ=255)                                                                                                                                                                                                                        aGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGagag                           <  1:223137/251‑1 (MQ=255)                                                                                                                                                                                                                          aaaTTGTTCCAGTAACTGTTGCTCACTCTCAAACGAAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCCGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATACCCTTCCATCGATTTACACAGTTAGTTGAGTATCAGTACGCTGGAGAGGc                         >  2:453553/1‑251 (MQ=255)                                                                                                                                                                                                                                 tCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGGGtttt                  >  2:126596/1‑251 (MQ=255)                                                                                                                                                                                                                                    aGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGcc               >  2:1825/1‑251 (MQ=255)                                                                                                                                                                                                                                    aGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCAAGCAAAGTGGCCCGATACAATTTAACGTGTCCGCGCAGTATGTTGGCGATACTATCGCCACCAAAATGCTGAAATTATCCAGCACTCAGCTGCCAGTTGAGGCGAGGTCGCTCGGGATGCCCTTCAATCGAATTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGCGTTTTGcc               >  1:249489/1‑251 (MQ=255)                                                                                                                                                                                                                                          tGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGATCAAAATGTCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTACATCGATTTAAACCGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTTCCtttttt         >  2:406806/1‑250 (MQ=255)                                                                                                                                                                                                                                          tGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTACATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGCGAGGCGCGTTTTGCCTTTTTc         >  2:387390/1‑251 (MQ=255)                                                                                                                                                                                                                                            ttGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCAACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGAATCAGTACGCTGGAGAGGCGAGTTTTGACTTTTTCAt       >  1:366807/1‑251 (MQ=255)                                                                                                                                                                                                                                                 cAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATgg  >  2:156886/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       CGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGG  >  NC_000913/10777‑11267

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