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,7600AG94.8% 217.6 / 2.0 58N431S (AAC→AGC) igaAinner membrane protein ‑ inhibits the Rcs signaling pathway
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (36/19);  total (38/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

ACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAACA  >  NC_000913/3527714‑3527804
                                              |                                            
aCCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcga                                           >  1:1211209/1‑50 (MQ=255)
 ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                          >  1:2700724/1‑50 (MQ=255)
  cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                         <  1:834486/50‑1 (MQ=255)
  cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                         >  1:466621/1‑50 (MQ=255)
  cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                         <  1:3805088/50‑1 (MQ=255)
  cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                         >  1:2969213/1‑50 (MQ=255)
  cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                         >  1:2818983/1‑50 (MQ=255)
   aaTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGc                                        >  1:279067/1‑50 (MQ=255)
   aaTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGc                                        >  1:1785021/1‑50 (MQ=255)
    aTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGcc                                       >  1:2294818/1‑50 (MQ=255)
    aTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGcc                                       <  1:1338837/50‑1 (MQ=255)
     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      <  1:198814/50‑1 (MQ=255)
     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      <  1:386947/50‑1 (MQ=255)
     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      <  1:2905587/50‑1 (MQ=255)
       cACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCTGc                                    >  1:2445757/1‑50 (MQ=255)
         ccTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                   >  1:3576344/1‑49 (MQ=255)
         ccTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTc                                  <  1:2631671/50‑1 (MQ=255)
           ttttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAt                                >  1:494989/1‑50 (MQ=255)
               ttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGc                             >  1:1609797/1‑49 (MQ=255)
                tACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCa                           >  1:4216269/1‑50 (MQ=255)
                tACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCa                           >  1:711594/1‑50 (MQ=255)
                 aCCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                          >  1:1721233/1‑50 (MQ=255)
                 aCCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                          >  1:3725584/1‑50 (MQ=255)
                 aCCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                          >  1:107487/1‑50 (MQ=255)
                  ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAtt                         <  1:2702670/50‑1 (MQ=255)
                  ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                          <  1:2668689/49‑1 (MQ=255)
                   cGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTa                        >  1:1240217/1‑50 (MQ=255)
                   cGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTa                        <  1:3727110/50‑1 (MQ=255)
                   cGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTa                        >  1:2927180/1‑50 (MQ=255)
                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:3003211/1‑50 (MQ=255)
                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:4351019/1‑50 (MQ=255)
                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       <  1:3490437/50‑1 (MQ=255)
                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:1012007/1‑50 (MQ=255)
                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:3587329/1‑50 (MQ=255)
                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       <  1:3548404/50‑1 (MQ=255)
                     tttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:1535518/1‑49 (MQ=255)
                      ttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:3443909/1‑48 (MQ=255)
                       tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:2025882/1‑47 (MQ=255)
                       tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:110251/1‑47 (MQ=255)
                       tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCa                     <  1:2946175/49‑1 (MQ=255)
                       tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCCCATTGCCATTAc                       >  1:1336335/1‑47 (MQ=255)
                       tGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTAc                       >  1:2538176/1‑47 (MQ=255)
                       tGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCa                     <  1:212009/49‑1 (MQ=255)
                        gACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAg                    >  1:2597319/1‑49 (MQ=255)
                          cTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAAt                 <  1:1712740/50‑1 (MQ=255)
                          cTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAAt                 <  1:4059197/50‑1 (MQ=255)
                           tGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATc                >  1:3125354/1‑50 (MQ=255)
                           tGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATc                >  1:1306443/1‑50 (MQ=255)
                             cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATc                <  1:3025444/48‑1 (MQ=255)
                             cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt               <  1:1599045/49‑1 (MQ=255)
                               cGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGaa            >  1:4190252/1‑50 (MQ=255)
                                gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGaa            <  1:1496325/49‑1 (MQ=255)
                                gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc           <  1:1869480/50‑1 (MQ=255)
                                gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc           >  1:108186/1‑50 (MQ=255)
                                gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc           >  1:1274407/1‑50 (MQ=255)
                                gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc           >  1:4182482/1‑50 (MQ=255)
                                 cAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt          <  1:1449345/50‑1 (MQ=255)
                                  aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg         >  1:134432/1‑50 (MQ=255)
                                  aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg         <  1:2636500/50‑1 (MQ=255)
                                  aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg         >  1:2883342/1‑50 (MQ=255)
                                   gATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTgg        >  1:2737514/1‑50 (MQ=255)
                                    atcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGt       >  1:4209398/1‑50 (MQ=255)
                                       atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa    >  1:894614/1‑50 (MQ=255)
                                          tGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaaca  >  1:1832528/1‑49 (MQ=255)
                                           ggAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaaca  >  1:2717783/1‑48 (MQ=255)
                                              |                                            
ACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAACA  >  NC_000913/3527714‑3527804

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: