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%	140.7 / NA	35	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 (22/13);  total (22/13)

CCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAA  >  NC_000913/3527715‑3527802                                              |                                           ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcga                                         >  1:1520955/1‑49 (MQ=255)   aaTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                       >  1:1297051/1‑49 (MQ=255)       cACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                    <  1:791330/48‑1 (MQ=255)       cACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCg                                   >  1:2007726/1‑49 (MQ=255)        aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCg                                   <  1:2059266/48‑1 (MQ=255)        aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGc                                  >  1:595617/1‑49 (MQ=255)          cTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                 >  1:1116499/1‑48 (MQ=255)           ttttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCa                               >  1:547968/1‑49 (MQ=255)           ttttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCa                               >  1:1258107/1‑49 (MQ=255)            tttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAt                              >  1:1305922/1‑49 (MQ=255)              tttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTg                            >  1:1777074/1‑49 (MQ=255)                  ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                        >  1:810385/1‑49 (MQ=255)                  ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                        >  1:1989280/1‑49 (MQ=255)                   cGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAtt                       <  1:1574865/49‑1 (MQ=255)                   cGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAtt                       <  1:1963207/49‑1 (MQ=255)                   cGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAtt                       <  1:572753/49‑1 (MQ=255)                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTa                      >  1:564426/1‑49 (MQ=255)                       tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAcc                    <  1:355099/48‑1 (MQ=255)                       tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                     >  1:698708/1‑47 (MQ=255)                       tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                     >  1:1164254/1‑47 (MQ=255)                        gACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                     >  1:141187/1‑46 (MQ=255)                         aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAg                  <  1:1860864/48‑1 (MQ=255)                         aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGa                 <  1:1436047/49‑1 (MQ=255)                          cTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGaa                >  1:1021601/1‑49 (MQ=255)                          cTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGaa                <  1:1134282/49‑1 (MQ=255)                            gCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATc              >  1:522755/1‑49 (MQ=255)                            gCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATc              >  1:1434133/1‑49 (MQ=255)                             cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATc              <  1:1172435/48‑1 (MQ=255)                              tCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt             >  1:1304719/1‑48 (MQ=255)                                gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGaa          >  1:578522/1‑49 (MQ=255)                                 cAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc         >  1:654337/1‑49 (MQ=255)                                  aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt        <  1:1585515/49‑1 (MQ=255)                                    atcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTgg      <  1:30041/49‑1 (MQ=255)                                     tcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGt     >  1:1221490/1‑49 (MQ=255)                                       atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTca   <  1:661490/49‑1 (MQ=255)                                       atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTc    >  1:623853/1‑48 (MQ=255)                                         cTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa  >  1:150123/1‑48 (MQ=255)                                              |                                           CCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAA  >  NC_000913/3527715‑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: ‑