Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,451,016	A→G	D20D (GAT→GAC)	yfcQ ←	putative fimbrial protein YfcQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,451,016	0	A	G	100.0%	102.8 / NA	30	D20D (GAT→GAC)	yfcQ	putative fimbrial protein YfcQ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/19);  total (11/19)

GACTTCAATGGTTTTACCACCGCTGATGCTACAGGTTGGCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCATCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGAAAGTCTTTCTCATTGATAGTCCACCAC  >  NC_000913/2450941‑2451089                                                                            |                                                                          gACTTCAATGGTTTTACCACCGCTGATGCTACAGGTTGGCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTcag                                                                       <  1:431811/80‑1 (MQ=255) gACTTCAATGGTTTTACCACCGCTGATGCTACAGGTTGGCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTcag                                                                       <  2:431812/80‑1 (MQ=255)      cAATGGTTTTACCACCGCTGATGCTACAGGTTGGCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTCAgcagca                                                                  <  2:431813/80‑1 (MQ=255)      cAATGGTTTTACCACCGCTGATGCTACAGGTTGGCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCAt                                                                 >  2:431857/1‑81 (MQ=255)           gTTTTACCACCGCTGATGCTACAGGTTGGCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGc                                                             <  2:431814/80‑1 (MQ=255)              ttACCACCGCTGATGCTACAGGTTGGCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTAt                                                          <  2:431815/80‑1 (MQ=255)               tACCACCGCTGATGCTACAGGTTGGCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATgg                                                        >  2:431863/1‑81 (MQ=255)               tACCACCGCTGATGCTACAGGTTGGCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTACGTCAGCAGCATGAGCTATg                                                         <  2:431816/80‑1 (MQ=255)                accaccGCTGATGCTACAGGTTGGCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGgc                                                       >  2:431851/1‑81 (MQ=255)                             cTACAGGTTGGCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACaa                                          <  1:431817/81‑1 (MQ=255)                                  ggTTGGCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAg                                     <  1:431818/81‑1 (MQ=255)                                    ttGGCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAggg                                   <  2:431819/81‑1 (MQ=255)                                     tGGCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGt                                  >  1:431859/1‑81 (MQ=255)                                      ggCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGt                                  <  2:431820/80‑1 (MQ=255)                                       gCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGtt                                 <  2:431821/80‑1 (MQ=255)                                       gCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGtt                                 <  2:431822/80‑1 (MQ=255)                                           tGACAATAATGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGa                             >  2:431865/1‑80 (MQ=255)                                                ataatGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGAAAGTCt                       >  2:431889/1‑81 (MQ=255)                                                 taatGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGAAAGTCtt                      <  2:431823/81‑1 (MQ=255)                                                  aatGTCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGAAAGTCttt                     >  1:431856/1‑81 (MQ=255)                                                      tCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGAAAGTCTTtctc                  <  2:431824/80‑1 (MQ=255)                                                      tCCCGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGAAAGTCTTTCTCa                 <  2:431825/81‑1 (MQ=255)                                                        ccGTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGAAAGTCTTTCTCAtt               >  2:431872/1‑81 (MQ=255)                                                          gTGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGAAAGTCTTTCTCATTGa             <  2:431826/81‑1 (MQ=255)                                                           tGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGAAAGTCTTTCTCATTGAt            <  1:431828/81‑1 (MQ=255)                                                           tGAAAGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGAAAGTCTTTCTCATTGAt            <  1:431827/81‑1 (MQ=255)                                                             aaaGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGAAAGTCTTTCTCATTGATAg          <  2:431829/81‑1 (MQ=255)                                                             aaaGGTGATGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGAAAGTCTTTCTCATTGATAg          >  2:431874/1‑81 (MQ=255)                                                                    aTGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGAAAGTCTTTCTCATTGATAGTccacca   >  1:431860/1‑81 (MQ=255)                                                                     tGTCTTCGTCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGAAAGTCTTTCTCATTGATAGTccaccac  >  1:431867/1‑81 (MQ=255)                                                                            |                                                                          GACTTCAATGGTTTTACCACCGCTGATGCTACAGGTTGGCGGTGACAATAATGTCCCGTGAAAGGTGATGTCTTCATCAGCAGCATGAGCTATGGCGCTACTGACACAAAGCAGGGTTAAGAAAGTCTTTCTCATTGATAGTCCACCAC  >  NC_000913/2450941‑2451089

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