Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	46,976	A→G	L390L (CTA→CTG)	yaaU →	putative transporter YaaU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	46,976	0	A	G	97.7%	133.7 / NA	44	L390L (CTA→CTG)	yaaU	putative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (26/17);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

CAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGTCGTATTGGCACCATTGTTTCGACCTGGGCACTACCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAGACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAA  >  NC_000913/46748‑47217                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  cAGCTTGGCCATGATGCCGCTGGAGCTGGCGGTTTTGGTGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCGTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGAACTGGGCACTGCCGATCTTTATCAATAATTACg                                                                                                                                                                                                                             <  2:164780/251‑1 (MQ=255)                           gtcgggTTTGGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATgggg                                                                                                                                                                                                   <  1:276034/247‑1 (MQ=255)                               ggTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGCACTCTTCCCGACCGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTCAGCCGTATTGGCACCCTTGTTTCGCCCTGGGCACTGCCGCTCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATgggg                                                                                                                                                                                                   >  2:276034/1‑247 (MQ=255)                                             tCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTAtctc                                                                                                                                                                                        <  1:619924/244‑1 (MQ=255)                                             tCCCGGATATGGGGATCTGGATGGTCGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGCTCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTAtctc                                                                                                                                                                                        >  2:619924/1‑244 (MQ=255)                                                cGGATGTGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTTTGCCTTTTTCTCTGGCTGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCGGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGTACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGc                                                                                                                                                                              <  1:629325/251‑1 (MQ=255)                                                      tGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCttgttg                                                                                                                                                                        >  1:489144/1‑251 (MQ=255)                                                               ggctggTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGt                                                                                                                                                                <  2:130141/250‑1 (MQ=255)                                                                   ggTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGt                                                                                                                                                           >  1:208595/1‑251 (MQ=255)                                                                     tAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTACTATGCAGTGGCTCTATCCTCATGCACTCTTCCCGAAAGATATCCGCGCCTCTGCAGTGGGCGTGATTCTGTCCTTAAGCCGTATTGGCACCATCGTTTCAACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCgatt                                                                                                                                                         >  1:603478/1‑248 (MQ=255)                                                                                      gtgtATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAgg                                                                                                                                        >  2:223585/1‑251 (MQ=255)                                                                                          aTGCCTTGTTCTCTGGCGGGCCGGGTAAGTTGCAGTGGCTCTATCCTAATGAACCCTTCCCGACAGTTACCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCTTATTGGTACCATTGTATCGACCGGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCGAATGGGGGCGGGTAGCTCGCTGTGGGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATg                                                                                                                                    <  1:268296/251‑1 (MQ=255)                                                                                            gCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGCGACTCGAGGGATGtt                                                                                                                                  >  2:16945/1‑250 (MQ=255)                                                                                                           gggCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAAcc                                                                                                                     >  1:424084/1‑249 (MQ=255)                                                                                                              ccGGGCAATTTGCAGTGGCTCTATCCTAATGACCTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGGGACTCGAGGGATGTCACTGGCGCAAACCAGCa                                                                                                                 >  2:625085/1‑250 (MQ=255)                                                                                                                gggTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAata                                                                                                              <  1:134995/251‑1 (MQ=255)                                                                                                                   tAATTTGCAGTGGCTCTGTCATAATGACATCTTCCCGACAGATATCCGCGCCTATGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCCCTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATg                                                                                                            >  2:319053/1‑250 (MQ=255)                                                                                                                         gCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGAGATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTTTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAGCCAGCAATATGACGATcc                                                                                                     <  1:301201/251‑1 (MQ=255)                                                                                                                               gCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGc                                                                                               <  2:424084/251‑1 (MQ=255)                                                                                                                                           tGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATCTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGt                                                                                   >  2:241300/1‑251 (MQ=255)                                                                                                                                              aCTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGGGCCTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGGGAATgggggaa                                                                                 >  2:123136/1‑247 (MQ=255)                                                                                                                                              aCTCTTCCCGACAGATAGCCGCGCCTCGGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTTCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCATAGAATGGGGTaaa                                                                                <  2:315783/251‑1 (MQ=255)                                                                                                                                                    cccGAAAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTc                                                                          >  2:607837/1‑251 (MQ=255)                                                                                                                                                      cGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTGAAGCCGTATGGGCACCATTGTTTCGACCTGGGCACTGCAGATGTTTATCAATAATTGCGGTATCAGGAACACGATGCTAATGGGGGGGGGTATCTCGCTGTTTGGCTTGTTGATGTCCGTAGCGTTTGCCCCGGAGACGCTATGGATGGCACTGGCGCAAACCAGCAATATGACGATCCGCTGGCAGAGAATGGGGTAAATTGTTCAg                                                                        <  1:366680/251‑1 (MQ=255)                                                                                                                                                          aGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGTGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGGAGCGTTTGCCCCGGAGACTCGAGGTATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTTTTCAGAttt                                                                    <  2:49796/251‑1 (MQ=255)                                                                                                                                                           gATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGACGCTAATGGGGGCGGGTAtct                                                                                                                                                                                         <  2:156233/133‑1 (MQ=255)                                                                                                                                                           gATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGACGCTAATGGGGGCGGGTAtct                                                                                                                                                                                         >  1:156233/1‑133 (MQ=255)                                                                                                                                                                   cgcCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCtttt                                                           >  1:151563/1‑251 (MQ=255)                                                                                                                                                                       tctGCCGTGGGCGTGATTATGTACTTAAGCCGTCTTGGCAGCATTGCTTCGACCTGGGCACTGCAGATCTTTATCAGTAATTACGGTATCAGTACCACGATGCTACTGGGGGCGGGTATCTCGCTGTTTGGCCTGTTGATTCCCGTAGCGTTTGCCACGGAGACTGGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGAAACCTGTTCAGATTTCACTCTTTTCTg                                                        >  2:299733/1‑250 (MQ=255)                                                                                                                                                                          gCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGa                                                       <  1:9264/248‑1 (MQ=255)                                                                                                                                                                             gTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAat                                                 <  1:16945/251‑1 (MQ=255)                                                                                                                                                                                    tGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCCAACCAGCAATATGACGCTCAGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTg                                          >  2:464690/1‑251 (MQ=255)                                                                                                                                                                                     gATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGCTCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTg                                          >  1:49930/1‑250 (MQ=255)                                                                                                                                                                                         aTGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGTGGGATGTCACTGGCGCAAACAAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTTTCTTTTCTGAATCAATATTATTGACtat                                     >  2:11533/1‑251 (MQ=255)                                                                                                                                                                                         aTGTCCTTAAGCCGTATTGGCACCATTGTTTCGAACTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTTTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGAAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACtat                                     <  2:397168/251‑1 (MQ=255)                                                                                                                                                                                                tAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGCCGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCcgc                              >  2:83614/1‑251 (MQ=255)                                                                                                                                                                                                 aagccGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTACTGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGCTCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTATCTTTTCTGAATCAATATTATTGACTATAAGCcgcg                             >  2:57603/5‑251 (MQ=255)                                                                                                                                                                                                     cGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGa                          >  1:352475/1‑250 (MQ=255)                                                                                                                                                                                                           ggCACCATTGTTTCGCCCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGCCGATCCGCGGGCCGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATg                   >  2:526516/1‑251 (MQ=255)                                                                                                                                                                                                               ccATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACt                >  2:235359/1‑250 (MQ=255)                                                                                                                                                                                                                cATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTac              >  1:122068/1‑251 (MQ=255)                                                                                                                                                                                                                    ggttCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGTGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACGGGCGCAATCCAGCAAGATGACGATCCGCGGGCAGATAATGGGGTAAATTGTTCAGATTGCTCTCTGTACTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACtt          <  2:90491/249‑1 (MQ=255)                                                                                                                                                                                                                            ctgggccatACCGATCTTTATCCATAATTCCGGTATCAGTAACACGATGCTAATGGGGGCGGGTATATCGCTGTTTGGCTTGATGATTTCCGTCGCGTTTGCCCCGGAGACTCGCGGGATGTCACTGGCGAAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGGAAATATATGACTACACTATGTGGGa   >  2:96859/9‑250 (MQ=255)                                                                                                                                                                                                                             tGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGaa  >  1:510918/1‑250 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  CAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGTCGTATTGGCACCATTGTTTCGACCTGGGCACTACCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAGACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAATATTATTGACTATAAGCCGCGTGAATATATGACTACACTTTGTGGGAA  >  NC_000913/46748‑47217

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