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%	153.9 / NA	40	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 (25/15);  total (25/15)

CCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAACA  >  NC_000913/3527715‑3527804                                              |                                             ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcga                                           >  1:1844116/1‑49 (MQ=255)    aTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGc                                        >  1:2557077/1‑49 (MQ=255)     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGcc                                       <  1:2469545/49‑1 (MQ=255)     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGcc                                       >  1:2435564/1‑49 (MQ=255)      tCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      >  1:1759962/1‑49 (MQ=255)       cACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      <  1:385629/48‑1 (MQ=255)       cACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      <  1:2273758/48‑1 (MQ=255)        aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGc                                    >  1:667604/1‑49 (MQ=255)         ccTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                   <  1:758571/49‑1 (MQ=255)         ccTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                   >  1:2122236/1‑49 (MQ=255)          cTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTc                                  <  1:2595452/49‑1 (MQ=255)           ttttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCa                                 >  1:923279/1‑49 (MQ=255)            tttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAt                                >  1:2516519/1‑49 (MQ=255)             ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCa                                 <  1:1955756/47‑1 (MQ=255)              tttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTg                              >  1:1964437/1‑49 (MQ=255)              tttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTg                              >  1:2854271/1‑49 (MQ=255)               ttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTg                              <  1:756193/48‑1 (MQ=255)                 aCCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGcc                            <  1:89692/48‑1 (MQ=255)                   cGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAtt                         <  1:887558/49‑1 (MQ=255)                   cGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAtt                         <  1:2819826/49‑1 (MQ=255)                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTa                        >  1:1667686/1‑49 (MQ=255)                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTa                        >  1:2744978/1‑49 (MQ=255)                     tttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       <  1:1481443/49‑1 (MQ=255)                        gACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:1694662/1‑46 (MQ=255)                        gACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCa                     <  1:729730/48‑1 (MQ=255)                         aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAg                    >  1:2716009/1‑48 (MQ=255)                         aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAg                    >  1:2528105/1‑48 (MQ=255)                         aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGa                   <  1:1592445/49‑1 (MQ=255)                          cTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGaa                  >  1:95849/1‑49 (MQ=255)                           tGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAAt                 >  1:288069/1‑49 (MQ=255)                                gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGaa            <  1:1713846/49‑1 (MQ=255)                                 cAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc           >  1:563990/1‑49 (MQ=255)                                 cAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc           >  1:102668/1‑49 (MQ=255)                                  aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt          >  1:2190646/1‑49 (MQ=255)                                    atcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTgg        >  1:1940059/1‑49 (MQ=255)                                     tcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGt       <  1:96135/49‑1 (MQ=255)                                      catcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGt       >  1:2549334/1‑48 (MQ=255)                                        tcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa    >  1:544154/1‑49 (MQ=255)                                         cTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa    >  1:2796949/1‑48 (MQ=255)                                         cTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa    >  1:2374874/1‑48 (MQ=255)                                           ggAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaaca  >  1:2394527/1‑48 (MQ=255)                                              |                                             CCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAACA  >  NC_000913/3527715‑3527804

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