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%	114.0 / NA	29	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 (18/11);  total (18/11)

ACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAA  >  NC_000913/3527714‑3527802                                               |                                           aCCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcg                                          <  1:498632/49‑1 (MQ=255)    aaTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                       >  1:893527/1‑49 (MQ=255)      ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGcc                                     >  1:1303687/1‑49 (MQ=255)      ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGcc                                     <  1:433908/49‑1 (MQ=255)        cACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                    <  1:1046625/48‑1 (MQ=255)         aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGc                                  >  1:658676/1‑49 (MQ=255)          ccTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                 <  1:170957/49‑1 (MQ=255)           cTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                 >  1:1315628/1‑48 (MQ=255)               tttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAtt                             <  1:2064966/48‑1 (MQ=255)                ttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGc                           >  1:978542/1‑49 (MQ=255)                        tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                     >  1:1005453/1‑47 (MQ=255)                        tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                     >  1:1174907/1‑47 (MQ=255)                          aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAg                  >  1:3447/1‑48 (MQ=255)                            tGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAAt               >  1:1191664/1‑49 (MQ=255)                             gCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATc              >  1:860397/1‑49 (MQ=255)                              cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATc              <  1:590061/48‑1 (MQ=255)                                 gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTg            <  1:373844/47‑1 (MQ=255)                                 gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGaa          >  1:1737757/1‑49 (MQ=255)                                 gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGaa          >  1:1834200/1‑49 (MQ=255)                                 gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGaa          >  1:422072/1‑49 (MQ=255)                                 gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGaa          >  1:724125/1‑49 (MQ=255)                                   aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt        >  1:1642862/1‑49 (MQ=255)                                   aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt        <  1:1142683/49‑1 (MQ=255)                                   aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt        >  1:703142/1‑49 (MQ=255)                                    gATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg       <  1:1738212/49‑1 (MQ=255)                                      tcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGt     >  1:776830/1‑49 (MQ=255)                                      tcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGt     <  1:1513751/49‑1 (MQ=255)                                       catcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGt     >  1:472433/1‑48 (MQ=255)                                        atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTca   <  1:349749/49‑1 (MQ=255)                                        atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTc    <  1:1844865/48‑1 (MQ=255)                                         tcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa  >  1:855790/1‑49 (MQ=255)                                               |                                           ACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAA  >  NC_000913/3527714‑3527802

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