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_0009133,527,7600AG100.0% 163.9 / NA 43N431S (AAC→AGC) igaAinner membrane protein ‑ inhibits the Rcs signaling pathway
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/18);  total (25/18)

CCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAACA  >  NC_000913/3527715‑3527804
                                             |                                            
ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                          >  1:1908405/1‑50 (MQ=255)
   aTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGcc                                       >  1:1069317/1‑50 (MQ=255)
    ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      <  1:1058091/50‑1 (MQ=255)
     tCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCg                                     >  1:2566880/1‑50 (MQ=255)
       aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGc                                    <  1:1120487/49‑1 (MQ=255)
       aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGc                                    <  1:1440092/49‑1 (MQ=255)
       aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                   >  1:531896/1‑50 (MQ=255)
       aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                   >  1:1116065/1‑50 (MQ=255)
       aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                   <  1:2054535/50‑1 (MQ=255)
       aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                   <  1:2503368/50‑1 (MQ=255)
         cTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCa                                 >  1:1472374/1‑50 (MQ=255)
         cTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCa                                 >  1:353139/1‑50 (MQ=255)
             tttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGc                             >  1:1311417/1‑50 (MQ=255)
               tACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGcc                            <  1:2531553/49‑1 (MQ=255)
                aCCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCa                           <  1:2673290/49‑1 (MQ=255)
                 ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAtt                         <  1:1406130/50‑1 (MQ=255)
                  cGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTa                        <  1:349092/50‑1 (MQ=255)
                  cGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTa                        >  1:1106305/1‑50 (MQ=255)
                  cGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTa                        <  1:1813922/50‑1 (MQ=255)
                   gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:1891214/1‑50 (MQ=255)
                     ttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:1081313/1‑48 (MQ=255)
                      tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCa                     <  1:2692737/49‑1 (MQ=255)
                        aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGaa                  >  1:702876/1‑50 (MQ=255)
                         cTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAAt                 >  1:971194/1‑50 (MQ=255)
                           gCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt               >  1:1721063/1‑50 (MQ=255)
                            cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt               <  1:622075/49‑1 (MQ=255)
                            cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTg              <  1:1862071/50‑1 (MQ=255)
                            cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTg              <  1:856648/50‑1 (MQ=255)
                             tCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTg              >  1:2520235/1‑49 (MQ=255)
                               gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGa             <  1:1691540/48‑1 (MQ=255)
                               gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGa             <  1:1574261/48‑1 (MQ=255)
                               gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc           >  1:1811271/1‑50 (MQ=255)
                               gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc           >  1:714665/1‑50 (MQ=255)
                                cAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt          >  1:1027090/1‑50 (MQ=255)
                                 aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg         <  1:1889859/50‑1 (MQ=255)
                                 aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg         >  1:2120667/1‑50 (MQ=255)
                                 aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg         >  1:193127/1‑50 (MQ=255)
                                  gATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTgg        >  1:1636720/1‑50 (MQ=255)
                                    tcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGt       >  1:156955/1‑49 (MQ=255)
                                      atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa    >  1:1715774/1‑50 (MQ=255)
                                      atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTca     <  1:2375860/49‑1 (MQ=255)
                                        cTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaac   >  1:170331/1‑49 (MQ=255)
                                         tGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaaca  >  1:2719926/1‑49 (MQ=255)
                                          ggAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaaca  >  1:670714/1‑48 (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 < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: