Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,527,760	A→G	N431S (AAC→AGC)	igaA →	inner membrane protein ‑ inhibits the Rcs signaling pathway

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,527,760	0	A	G	100.0%	72.1 / NA	19	N431S (AAC→AGC)	igaA	inner membrane protein ‑ inhibits the Rcs signaling pathway
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/7);  total (12/7)

CCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAACA  >  NC_000913/3527715‑3527804                                              |                                             ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                          >  1:957025/1‑50 (MQ=255) ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                          <  1:1785971/50‑1 (MQ=255)  cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                         >  1:2198101/1‑50 (MQ=255)     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      <  1:1184385/50‑1 (MQ=255)        aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                   <  1:2248620/50‑1 (MQ=255)          cTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCa                                 <  1:435579/50‑1 (MQ=255)           ttttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAt                                >  1:2248984/1‑50 (MQ=255)             ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAt                                <  1:2129938/48‑1 (MQ=255)             ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAt                                <  1:43162/48‑1 (MQ=255)               ttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGc                             <  1:1615357/49‑1 (MQ=255)               ttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGc                             >  1:1498456/1‑49 (MQ=255)                 aCCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                          >  1:559886/1‑50 (MQ=255)                  ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                          >  1:78069/1‑49 (MQ=255)                     tttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:1007795/1‑49 (MQ=255)                      ttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:1770742/1‑48 (MQ=255)                            gCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt               >  1:2055897/1‑50 (MQ=255)                                 cAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt          <  1:4405/50‑1 (MQ=255)                                  aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg         >  1:903717/1‑50 (MQ=255)                                      catcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTc      >  1:470963/1‑49 (MQ=255)                                       atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa    >  1:1657831/1‑50 (MQ=255)                                         cTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaac   >  1:62636/1‑49 (MQ=255)                                          tGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaaca  >  1:2136826/1‑49 (MQ=255)                                              |                                             CCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAACA  >  NC_000913/3527715‑3527804

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