Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,459,326	G→A	V330I (GTT→ATT)	gspE →	Type II secretion system protein GspE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,459,326	0	G	A	100.0%	168.5 / NA	52	V330I (GTT→ATT)	gspE	Type II secretion system protein GspE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/26);  total (26/26)

GGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGATGTTTAAATATCATCAGCTCGCGGTGACAACAATTG  >  NC_000913/3459083‑3459566                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 gggCCGACAGGCTCCGGTAAAAGCACCACGCTATACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACACGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGc                                                                                                                                                                                                                                           >  1:223234/1‑251 (MQ=255)   gCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCCGGCCt                                                                                                                                                                                                                                         >  2:126500/1‑251 (MQ=255)    ccGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTc                                                                                                                                                                                                                                        <  1:186107/251‑1 (MQ=255)      gACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGTTGTCGTCATGGCGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGc                                                                                                                                                                                                                                      >  1:416321/1‑251 (MQ=255)       aCAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCt                                                                                                                                                                                                                                     <  1:160578/251‑1 (MQ=255)                 gTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGTCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCAt                                                                                                                                                                                                                           <  1:101037/251‑1 (MQ=255)                      aGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGt                                                                                                                                                                                                                      <  2:336885/251‑1 (MQ=255)                           cacGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGGGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGACCTCGCTCACCGGCCATCTGGTActct                                                                                                                                                                                                                 >  1:87234/1‑251 (MQ=255)                               ctctACGCCATACTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAAACGCGAAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAac                                                                                                                                                                                                             >  1:294635/1‑251 (MQ=255)                                        aTCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCacac                                                                                                                                                                                                    <  1:27152/251‑1 (MQ=255)                                           cTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTa                                                                                                                                                                                                  >  1:409998/1‑250 (MQ=255)                                             ttCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAAc                                                                                                                                                                                                >  1:95171/1‑250 (MQ=255)                                                    cTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATc                                                                                                                                                                                        >  1:202774/1‑251 (MQ=255)                                                          aCTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGgcgc                                                                                                                                                                                  >  1:341895/1‑251 (MQ=255)                                                                    gCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCgg                                                                                                                                                                        <  1:94762/251‑1 (MQ=255)                                                                        tatTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTcc                                                                                                                                                                    >  1:290891/1‑251 (MQ=255)                                                                         atTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCcg                                                                                                                                                                   <  1:173864/251‑1 (MQ=255)                                                                              gACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACa                                                                                                                                                              <  2:221671/251‑1 (MQ=255)                                                                              gACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGAACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGAc                                                                                                                                                               >  1:224208/1‑250 (MQ=255)                                                                                  gTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCATATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATggg                                                                                                                                                          <  1:204112/251‑1 (MQ=255)                                                                                   tAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGGCTGCACGTGTGGATAGGTGTTTAGCTCGCGGCCTGCGCGCCATACTTCGCCAGGAACCGGATGTAGACATGGGGGGGTGAATTCGTGCAAAAGAAACCGCGCAGCTTGCGATTCAGGCTTAGCTCAACGGGCATGTGGTACTCTCAACACTCCACGCTAACGGTGTATCAGGAGCAgataccctgtgccgccacaggggc                                                                                                                                                         >  2:89846/1‑228 (MQ=255)                                                                                      aagaTCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACGGTTGTGGATATGTCGTTCGCTCGCGGCCTGCGCGCAGTTGTTCGACAGAACCCGGATGTCGTCGTGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCGCGCTCACCGGCGGGGTGGTACTTTCAAGCCTGCACACTAACAGTGCATAAGGCGCGGTGAACTGGCTCCGTGACGTGGACGTc                                                                                                                                                      <  2:246013/251‑1 (MQ=255)                                                                                      aagaTCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAGCCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACGCCACGCTAACAGTGCAGCAGGCGCAGTGACCCTGCTGCGCGACATGGGCGTc                                                                                                                                                      <  2:230823/251‑1 (MQ=255)                                                                                           cccGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATc                                                                                                                                                 <  1:34207/251‑1 (MQ=255)                                                                                              gTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCAtt                                                                                                                                              >  1:263255/1‑251 (MQ=255)                                                                                                ggAATATGAGCTGGAAGGCATTGGTCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAATTTCGTGATACAGAAACCGCGCAGATTGCGATGCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGTCTCCGCGACATGGGCGTCGAATCATTc                                                                                                                                             <  1:415893/250‑1 (MQ=255)                                                                                                  aaTATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTg                                                                                                                                          <  2:416321/251‑1 (MQ=255)                                                                                                     atGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCtt                                                                                                                                       >  1:414960/1‑251 (MQ=255)                                                                                                           tGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCt                                                                                                                                 >  2:64691/1‑251 (MQ=255)                                                                                                            ggAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCt                                                                                                                                 <  2:95171/250‑1 (MQ=255)                                                                                                             gAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTc                                                                                                                               <  2:137969/251‑1 (MQ=255)                                                                                                             gAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCACTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCCCACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTc                                                                                                                               >  1:276936/1‑251 (MQ=255)                                                                                                                 gCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTg                                                                                                                           <  2:414960/251‑1 (MQ=255)                                                                                                                      gggCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATTCTTCGCCAGGACCCGGATGTCGTCATTGTGGGTGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCTCACTAACAGTGCATCAGGCGCAGTGACCGGTCTCCTCGACATGGGCGTCGAATCATTCCTGCTGTCGTCTTCCCTGGCAg                                                                                                                       <  2:196357/250‑1 (MQ=255)                                                                                                                              gaaGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTGGCCAAGACCCGGATGTCGTCATGGTGGGGGAAAATCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCGTCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCCGTGACCCGGCTCCGCGCCAGTGTCGTCGAATCATTCCTGCTTTCGTCTTCCCGAGCAGGGATTATcg                                                                                                              >  2:299200/3‑251 (MQ=255)                                                                                                                              gCAGGTGAATACCTGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATcg                                                                                                              <  1:145996/251‑1 (MQ=255)                                                                                                                                  gTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCa                                                                                                          <  1:114872/251‑1 (MQ=255)                                                                                                                                          ccGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTgg                                                                                                  <  1:52565/251‑1 (MQ=255)                                                                                                                                                      gTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCtgt                                                                                      >  2:232159/1‑251 (MQ=255)                                                                                                                                                              ctcgcGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCgta                                                                               >  2:109378/1‑248 (MQ=255)                                                                                                                                                                    gCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGa                                                                        <  2:95180/251‑1 (MQ=255)                                                                                                                                                                     ccTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGAc                                                                       >  1:405940/1‑251 (MQ=255)                                                                                                                                                                            gcCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGTTACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCAc                                                                >  1:131813/1‑251 (MQ=255)                                                                                                                                                                                tACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGccc                                                            >  1:366597/1‑251 (MQ=255)                                                                                                                                                                                                       gtcATGGCGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAgcgc                                         >  2:367125/1‑247 (MQ=255)                                                                                                                                                                                                         aaTGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGAGTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCTTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCCCCCCAACAAGCTCAGATGt                                   <  2:401341/250‑1 (MQ=255)                                                                                                                                                                                                          aTGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGATGtt                                  >  1:339047/1‑251 (MQ=255)                                                                                                                                                                                                           tggtggGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTAATCTCAACACTCCACACTAACAGTGCATCAGGAGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGAAACGTCTGGTTCGTCGCATGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCCACAAGCGCAGATGttt                                 >  1:364816/1‑251 (MQ=255)                                                                                                                                                                                                                      aTTCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGATGTTTAAATAtcatca                      <  2:411717/251‑1 (MQ=255)                                                                                                                                                                                                                        tCGTGATACAGAAACCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCGCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGATGTTTAAATATCATCAg                     <  2:391983/250‑1 (MQ=255)                                                                                                                                                                                                                                   aaaCCGCGCAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGATGTTTAAATATCATCAGCTCGCGGTGac          >  2:163951/1‑250 (MQ=255)                                                                                                                                                                                                                                          gcAGATTGCGATTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACCCTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGATGTTTAAATATCATCAGCTCGTGGTGACAACAATTg  <  2:276936/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 GGGCCGACAGGCTCCGGTAAAAGCACCACGCTCTACGCCATCCTTTCGGCGCTGAATACTCCCGGCCGCAATATTCTGACGGTAGAAGATCCCGTGGAATATGAGCTGGAAGGCATTGGGCAAACGCAGGTGAATACCCGTGTGGATATGTCTTTCGCTCGCGGCCTGCGCGCCATACTTCGCCAGGACCCGGATGTCGTCATGGTGGGGGAAATTCGTGATACAGAAACCGCGCAGATTGCGGTTCAGGCCTCGCTCACCGGCCATCTGGTACTCTCAACACTCCACACTAACAGTGCATCAGGCGCAGTGACCCGGCTCCGCGACATGGGCGTCGAATCATTCCTGCTTTCGTCTTCCCTGGCAGGGATTATCGCGCAACGTCTGGTTCGTCGCCTGTGTCCGCAATGCCGACAATTCACGCCCGTATCACCCCAACAAGCGCAGATGTTTAAATATCATCAGCTCGCGGTGACAACAATTG  >  NC_000913/3459083‑3459566

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