Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,911,287	G→A	A338T (GCT→ACT)	yebQ →	putative transporter YebQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,911,287	0	G	A	96.6%	183.4 / ‑1.4	60	A338T (GCT→ACT)	yebQ	putative transporter YebQ
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (29/28);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.33e-01

CCGTACCGCTGCTACCGGTGGATTTACTGCGTATCCCGCTGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGGCTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGAACGGCTCGTCTACTGGGTCAGAGTAGCGGCGCGGCGCTGGTGGCGCTGATG  >  NC_000913/1911042‑1911532                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      ccGTACCGCTGCTACCGGTGGATTTACTGCGTATCCCGCTGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACttt                                                                                                                                                                                                                                                   >  1:294432/1‑250 (MQ=255)   gTACCGCTTCTACCGGTGGATTTACTGCGTATCCCGCTGTTTTCACTTTCTATTTGCACATCTGTGTGCTCTTACTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGACTATTTGATTGAACGCGTACATGCCGGATTGCTGGGGACTTTAgt                                                                                                                                                                                                                                                <  2:519371/251‑2 (MQ=255)    tACCGCTGCTACCGGTGGATTTACTGCGTATCCCGCTGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAggg                                                                                                                                                                                                                                               >  1:342098/1‑251 (MQ=255)                cGGTGGATTTACTGCGTATCCCGCTGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATg                                                                                                                                                                                                                                   >  1:494730/1‑251 (MQ=255)                 ggTGGATTTACTGCGTATCCCGCTGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCCGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACCTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATgg                                                                                                                                                                                                                                  <  1:304247/251‑1 (MQ=255)                         tACTGCGTATCCCGCTGTTTTCACTTTCTATTTGCACGTCTGTTTGCTCTTTCTGCGCGCAAATGCTGGCGATGGTTGCCCTGCGCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGTTTTACTTCTGACGCCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGAGTGCTGGAGTCTTTAGGGTTGTTCATCATGGCTGCggg                                                                                                                                                                                                                           <  2:128493/250‑1 (MQ=255)                          aCTGCGTATCCCGCTGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCgggg                                                                                                                                                                                                                          <  1:158733/250‑1 (MQ=255)                              cGTATCCCGCTGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCtttt                                                                                                                                                                                                                     <  1:545950/251‑1 (MQ=255)                                tATCCCGCTGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCtttttt                                                                                                                                                                                                                   >  1:390505/1‑251 (MQ=255)                                    ccGCTGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCt                                                                                                                                                                                                               >  1:375045/1‑251 (MQ=255)                                     cGCTGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAACTGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCTCTg                                                                                                                                                                                                              >  2:475194/1‑251 (MQ=255)                                     cGCTGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTg                                                                                                                                                                                                              <  2:624816/251‑1 (MQ=255)                                      gCTGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTgg                                                                                                                                                                                                             >  2:521746/1‑251 (MQ=255)                                         gTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGGCTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGt                                                                                                                                                                                                            <  1:234232/249‑1 (MQ=255)                                                                  gTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGat                                                                                                                                                                                   <  1:249998/249‑1 (MQ=255)                                                                    ttGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTTCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGTGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATc                                                                                                                                                                                <  2:345071/250‑1 (MQ=255)                                                                                     aaaTGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCgatg                                                                                                                                                              <  2:371037/251‑1 (MQ=255)                                                                                      aaTGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCgatga                                                                                                                                                             <  2:545471/251‑1 (MQ=255)                                                                                       aTGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCgatgat                                                                                                                                                            >  1:630801/1‑251 (MQ=255)                                                                                         gCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATc                                                                                                                                                           <  1:302269/250‑1 (MQ=255)                                                                                              cAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTAtgt                                                                                                                                                     >  1:601081/1‑251 (MQ=255)                                                                                               aaTNGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTAtgtg                                                                                                                                                    <  2:601081/251‑1 (MQ=255)                                                                                               aaTGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGGCTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGCTGATCTTAtgtg                                                                                                                                                    >  2:491058/1‑251 (MQ=255)                                                                                               aaTGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTAtgtg                                                                                                                                                    >  1:521784/1‑251 (MQ=255)                                                                                                    tttCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGCCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTCTTTGATTGAAAGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTCCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCt                                                                                                                                               >  1:406054/1‑251 (MQ=255)                                                                                                       cccTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTgg                                                                                                                                             >  2:243749/1‑250 (MQ=255)                                                                                                             ccTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGc                                                                                                                                      >  2:10542/1‑251 (MQ=255)                                                                                                              cTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCATGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGTTGCTGGATTTGGCt                                                                                                                                     >  1:426454/1‑251 (MQ=255)                                                                                                                               tGCTCGGGCGTAGTGAAGTCGAAACAGGATTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGATGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAAt                                                                                                                    >  1:46867/1‑251 (MQ=255)                                                                                                                                    gggCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGCTGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCAAGTCACCCAATAAcc                                                                                                                >  2:503803/1‑250 (MQ=255)                                                                                                                                        gTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACAcc                                                                                                           <  1:265892/251‑1 (MQ=255)                                                                                                                                         tAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCCCCGCTGGCAGGCTATCTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTTCCCGCGTCACCTGCGGGTATCAATATTATCTGGCCGACGATCTTATGTGGTGCTGGATGTGGATTATTCCAGTCACCCAATAACCACAcc                                                                                                           >  2:345424/1‑250 (MQ=255)                                                                                                                                            tGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCatt                                                                                                        >  1:601233/1‑250 (MQ=255)                                                                                                                                                       cAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCcgcg                                                                                            >  2:96968/1‑251 (MQ=255)                                                                                                                                                       cAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCcgcg                                                                                            >  2:576440/1‑251 (MQ=255)                                                                                                                                                              aCTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTcgcg                                                                                     >  1:527563/1‑251 (MQ=255)                                                                                                                                                                    gACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCATGCTATTTGAGTGAACTCGTACATGCAGGATTGCTGGGTACTTTACGGTTGTTCATCATGGCTGCGTGGCTTTTGTCCCTGGTTCTGCTGCCCGCGTCACCTGAGGATATCAAGATTATCGGGCCGATGATCTTATGTGGTGCTGGATTGGGCTTAGTCCAGTCGCCCAATAGCCACAGCATCAGTTGCTCCGCGCGTCGCGAACGt                                                                                <  2:54832/250‑1 (MQ=255)                                                                                                                                                                    gACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATGGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAAc                                                                                  <  2:98802/248‑1 (MQ=255)                                                                                                                                                                        ccGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCgg                                                                           <  1:85597/251‑1 (MQ=255)                                                                                                                                                                          gTGGCCGTTAGCAACGCTGGTGATGGCTCCGCTGACAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTg                                                                         >  1:215041/1‑251 (MQ=255)                                                                                                                                                                          gTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTg                                                                         <  2:342098/251‑1 (MQ=255)                                                                                                                                                                          gTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTCCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTCGCGGTg                                                                         >  2:224540/1‑251 (MQ=255)                                                                                                                                                                                        cGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATg                                                           <  2:563001/251‑1 (MQ=255)                                                                                                                                                                                               gATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGa                                                     >  1:469154/1‑250 (MQ=255)                                                                                                                                                                                                    cTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGAGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACATCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGAACGGcg                                               >  2:498533/1‑250 (MQ=255)                                                                                                                                                                                                     tCCGCTGGCAGGCTATTTGATTGAACGCGCACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGAACGGCTc                                              >  1:513475/1‑251 (MQ=255)                                                                                                                                                                                                      ccGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGAACGGCTCg                                             <  2:294432/251‑1 (MQ=255)                                                                                                                                                                                                       cGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGAACGGCTCGt                                            <  2:509415/251‑1 (MQ=255)                                                                                                                                                                                                         cTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTAt                                                                                                                                  <  2:203239/163‑1 (MQ=255)                                                                                                                                                                                                         cTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTAt                                                                                                                                  >  1:203239/1‑163 (MQ=255)                                                                                                                                                                                                             cAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCCCCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGAACGGCTCGTCTACTg                                      <  2:464911/251‑1 (MQ=255)                                                                                                                                                                                                              aGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGAACGGCTCGTCTACTgg                                     <  1:10542/251‑1 (MQ=255)                                                                                                                                                                                                                gCTATGTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGAACGGCTCGTCTACTGGGt                                   <  2:513475/251‑1 (MQ=255)                                                                                                                                                                                                                   aTTTGATTGAACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGAACGGCTCGTCTACTGGGTCa                                 <  1:96968/250‑1 (MQ=255)                                                                                                                                                                                                                     ttgattgaACGCGTACATGCAGGATTGCTGTGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCCCCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGAACGGCTCGTCTACTGGGTCaga                               <  1:405600/250‑1 (MQ=255)                                                                                                                                                                                                                       gattgaACGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGAACGGCTCGTCTACTGGGTCAGAGt                             <  2:106518/250‑1 (MQ=255)                                                                                                                                                                                                                          tgatcgcgTACATGCAGGATTGCTTGGTACTTTAAGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGTCTTTTTCCAGTCCCCTAATAACCACACCATTATTTCCTCCGCGCCTCGCGACCGTTGCGGTGGAGCCAGTTGCATGTTAGGAACGGCTCGTCTACTGGGTAAGAGTAgc                          <  1:92683/246‑1 (MQ=255)                                                                                                                                                                                                                           taaCGCGTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCCATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGCACGGCTCATCTACTGGGGCAGAGTAgcgg                        >  2:517603/2‑251 (MQ=255)                                                                                                                                                                                                                                cgTACATGCAGGATTGCTGGGGACTTTAGGGTTGTTCATAATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGAACGGCCCGTCCATTGGGGCAGAGGAacggcgagg                   >  2:646394/1‑248 (MQ=255)                                                                                                                                                                                                                                   aCATGCAGGAGTGCTGGGGACTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGTCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGAACGGCTCGTCTACTGGGTCAGAGTAGCggcgcggcgc                <  2:243814/251‑1 (MQ=255)                                                                                                                                                                                                                                                  ggggaCTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGAACGGCTCGTCTACTGGGTCAGAGTAGCGGCGCGGCGCTGGTGGCGCTGATg  <  1:317148/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      CCGTACCGCTGCTACCGGTGGATTTACTGCGTATCCCGCTGTTTTCACTTTCTATTTGCACATCTGTTTGCTCTTTCTGCGCACAAATGCTGGCAATGGTTTCCCTGCCCTTTTACCTGCAAACCGTGCTCGGGCGTAGTGAAGTCGAAACAGGTTTACTTCTGACACCGTGGCCGTTAGCAACGATGGTGATGGCTCCGCTGGCAGGCTATTTGATTGAACGCGTACATGCAGGATTGCTGGGGGCTTTAGGGTTGTTCATCATGGCTGCGGGGCTTTTTTCCCTGGTTCTGCTGCCCGCGTCACCTGCGGATATCAATATTATCTGGCCGATGATCTTATGTGGTGCTGGATTTGGCTTATTCCAGTCACCCAATAACCACACCATTATTACCTCCGCGCCTCGCGAACGTAGCGGTGGAGCCAGTGGCATGTTAGGAACGGCTCGTCTACTGGGTCAGAGTAGCGGCGCGGCGCTGGTGGCGCTGATG  >  NC_000913/1911042‑1911532

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