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,7600AG98.0% 198.9 / ‑3.5 51N431S (AAC→AGC) igaAinner membrane protein ‑ inhibits the Rcs signaling pathway
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (35/15);  total (35/16)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

ACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAACA  >  NC_000913/3527714‑3527804
                                              |                                            
aCCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcga                                           <  1:2443525/50‑1 (MQ=255)
 ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                          >  1:406956/1‑50 (MQ=255)
  cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                         >  1:1154072/1‑50 (MQ=255)
  cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                         >  1:1237320/1‑50 (MQ=255)
    aTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGcc                                       >  1:2737812/1‑50 (MQ=255)
     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      >  1:1320988/1‑50 (MQ=255)
     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      >  1:512931/1‑50 (MQ=255)
     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      >  1:1468425/1‑50 (MQ=255)
     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      <  1:1920886/50‑1 (MQ=255)
        aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                   >  1:866091/1‑50 (MQ=255)
         ccTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                   >  1:3481433/1‑49 (MQ=255)
         ccTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTc                                  <  1:2971300/50‑1 (MQ=255)
         ccTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTc                                  <  1:690603/50‑1 (MQ=255)
           ttttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAt                                >  1:2795384/1‑50 (MQ=255)
             ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTg                              >  1:527310/1‑50 (MQ=255)
             ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTg                              >  1:2527070/1‑50 (MQ=255)
             ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTg                              >  1:3531204/1‑50 (MQ=255)
                tACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCa                           >  1:1441639/1‑50 (MQ=255)
                  ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAtt                         <  1:2614987/50‑1 (MQ=255)
                  ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                          >  1:2900938/1‑49 (MQ=255)
                  ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                          >  1:2313805/1‑49 (MQ=255)
                  ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                          <  1:705664/49‑1 (MQ=255)
                   cGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTa                        <  1:1403626/50‑1 (MQ=255)
                   cGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTa                        <  1:2350928/50‑1 (MQ=255)
                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:3895817/1‑50 (MQ=255)
                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       <  1:1480144/50‑1 (MQ=255)
                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:776607/1‑50 (MQ=255)
                     tttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:611676/1‑49 (MQ=255)
                     tttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:275333/1‑49 (MQ=255)
                      ttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:3317032/1‑48 (MQ=255)
                       tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:1914027/1‑47 (MQ=255)
                       tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:654250/1‑47 (MQ=255)
                        gACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAg                    >  1:2875497/1‑49 (MQ=255)
                        gACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAg                    >  1:3691823/1‑49 (MQ=255)
                         aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGaa                  <  1:683125/50‑1 (MQ=255)
                          cTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAAt                 <  1:2939610/50‑1 (MQ=255)
                          cTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAAt                 <  1:1729472/50‑1 (MQ=255)
                          cTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAAt                 >  1:2241127/1‑50 (MQ=255)
                            gCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt               >  1:1924192/1‑50 (MQ=255)
                             cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt               <  1:1219623/49‑1 (MQ=255)
                              tCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTg              >  1:1634892/1‑49 (MQ=255)
                                gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGaa            <  1:1297694/49‑1 (MQ=255)
                                gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc           >  1:3078456/1‑50 (MQ=255)
                                 cAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt          <  1:2368333/50‑1 (MQ=255)
                                 cAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt          <  1:1049587/50‑1 (MQ=255)
                                  aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg         >  1:297723/1‑50 (MQ=255)
                                  aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg         >  1:765138/1‑50 (MQ=255)
                                   gATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTgg        >  1:2337279/1‑50 (MQ=255)
                                       atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa    >  1:656342/1‑50 (MQ=255)
                                       atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa    >  1:3688354/1‑50 (MQ=255)
                                       atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTca     <  1:3332084/49‑1 (MQ=255)
                                         cTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaac   >  1:1324772/1‑49 (MQ=255)
                                           ggAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaaca  >  1:3443935/1‑48 (MQ=255)
                                           ggAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaaca  >  1:3326134/1‑48 (MQ=255)
                                           ggAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaaca  >  1:873528/1‑48 (MQ=255)
                                              |                                            
ACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAACA  >  NC_000913/3527714‑3527804

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