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.5% 263.6 / ‑3.7 66N431S (AAC→AGC) igaAinner membrane protein ‑ inhibits the Rcs signaling pathway
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (47/18);  total (47/19)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

AACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAACAAAGCG  >  NC_000913/3527713‑3527809
                                               |                                                 
aaCCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcg                                                 <  1:1045321/50‑1 (MQ=255)
 aCCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcga                                                >  1:2540136/1‑50 (MQ=255)
  ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                               <  1:3823197/50‑1 (MQ=255)
  ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                               <  1:862162/50‑1 (MQ=255)
   cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                              >  1:1189083/1‑50 (MQ=255)
   cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                              >  1:4894823/1‑50 (MQ=255)
   cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                              >  1:430755/1‑50 (MQ=255)
    aaTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGc                                             >  1:5033630/1‑50 (MQ=255)
     aTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGcc                                            >  1:5854921/1‑50 (MQ=255)
      ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                           <  1:4709297/50‑1 (MQ=255)
      ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                           <  1:5327845/50‑1 (MQ=255)
       tCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCg                                          >  1:5436361/1‑50 (MQ=255)
        cACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGc                                         >  1:5594038/1‑50 (MQ=255)
         aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                        >  1:3947703/1‑50 (MQ=255)
         aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                        <  1:3414497/50‑1 (MQ=255)
          ccTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                        >  1:148563/1‑49 (MQ=255)
           cTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCa                                      <  1:3601151/50‑1 (MQ=255)
           cTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCa                                      <  1:5212021/50‑1 (MQ=255)
           cTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCa                                      >  1:1990642/1‑50 (MQ=255)
            ttttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAt                                     >  1:101963/1‑50 (MQ=255)
            ttttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAt                                     >  1:5605319/1‑50 (MQ=255)
             tttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAtt                                    >  1:5477586/1‑50 (MQ=255)
              ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAtt                                    <  1:3027596/49‑1 (MQ=255)
              ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAt                                     <  1:5612096/48‑1 (MQ=255)
              ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTg                                   >  1:5414580/1‑50 (MQ=255)
               tttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGc                                  >  1:4332663/1‑50 (MQ=255)
               tttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGc                                  >  1:5288365/1‑50 (MQ=255)
               tttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGc                                  >  1:4675259/1‑50 (MQ=255)
                 tACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCa                                >  1:1305092/1‑50 (MQ=255)
                  aCCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCa                                <  1:693950/49‑1 (MQ=255)
                  aCCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                               >  1:3763786/1‑50 (MQ=255)
                  aCCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                               >  1:1510742/1‑50 (MQ=255)
                   ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                               <  1:4297985/49‑1 (MQ=255)
                      tttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAcc                           <  1:1954773/50‑1 (MQ=255)
                      tttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                            >  1:2210508/1‑49 (MQ=255)
                       ttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                            >  1:1364923/1‑48 (MQ=255)
                       ttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                            >  1:746663/1‑48 (MQ=255)
                        tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                            >  1:4768715/1‑47 (MQ=255)
                         gACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAg                         >  1:4447392/1‑49 (MQ=255)
                         gACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAg                         <  1:5775697/49‑1 (MQ=255)
                          aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGaa                       >  1:1365280/1‑50 (MQ=255)
                          aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGaa                       >  1:5056700/1‑50 (MQ=255)
                          aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGaa                       >  1:1476645/1‑50 (MQ=255)
                          aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGaa                       >  1:4037015/1‑50 (MQ=255)
                          aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGa                        <  1:5921396/49‑1 (MQ=255)
                             gCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt                    >  1:3344180/1‑50 (MQ=255)
                             gCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt                    >  1:1771807/1‑50 (MQ=255)
                              cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt                    <  1:4391137/49‑1 (MQ=255)
                              cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTg                   <  1:4909961/50‑1 (MQ=255)
                               tCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTg                   >  1:507041/1‑49 (MQ=255)
                               tCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTg                   >  1:1000646/1‑49 (MQ=255)
                               tCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTg                   >  1:5055820/1‑49 (MQ=255)
                                cGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGaa                 >  1:5259960/1‑50 (MQ=255)
                                 gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGa                  <  1:2561296/48‑1 (MQ=255)
                                 gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc                >  1:5765532/1‑50 (MQ=255)
                                 gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc                >  1:4528236/1‑50 (MQ=255)
                                 gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc                <  1:2974165/50‑1 (MQ=255)
                                  cAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt               >  1:1757402/1‑50 (MQ=255)
                                  cAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt               >  1:5324615/1‑50 (MQ=255)
                                   aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg              >  1:3494259/1‑50 (MQ=255)
                                   aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg              <  1:1681722/50‑1 (MQ=255)
                                   aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg              >  1:5665420/1‑50 (MQ=255)
                                   aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg              >  1:5219585/1‑50 (MQ=255)
                                   aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg              >  1:475126/1‑50 (MQ=255)
                                    gATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTgg             <  1:5194409/50‑1 (MQ=255)
                                      tcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGt            >  1:5417167/1‑49 (MQ=255)
                                      tcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGt            >  1:864541/1‑49 (MQ=255)
                                       catcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTc           >  1:1510814/1‑49 (MQ=255)
                                        atcTGGCGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTca          <  1:1742261/49‑1 (MQ=38)
                                        atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa         >  1:4240758/1‑50 (MQ=255)
                                         tcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa         >  1:3209831/1‑49 (MQ=255)
                                           tGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaaca       >  1:266151/1‑49 (MQ=255)
                                            ggAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaaca       >  1:1709299/1‑48 (MQ=255)
                                               acgacgCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAACAAAGCg  <  1:4568119/50‑1 (MQ=255)
                                               |                                                 
AACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAACAAAGCG  >  NC_000913/3527713‑3527809

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: