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.4% 502.7 / ‑2.0 126N431S (AAC→AGC) igaAinner membrane protein ‑ inhibits the Rcs signaling pathway
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (64/60);  total (66/60)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-02

Only 100 of 141 total aligned reads displayed.

ACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAAC  >  NC_000913/3527714‑3527803
                                              |                                           
aCCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcga                                          >  2:2166600/1‑50 (MQ=255)
 ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                         <  2:4260593/50‑1 (MQ=255)
 ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                         >  2:4061231/1‑50 (MQ=255)
 ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                         >  2:1596157/1‑50 (MQ=255)
 ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                         <  1:2579846/50‑1 (MQ=255)
 ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                         <  1:1952720/50‑1 (MQ=255)
  cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                        >  1:2219897/1‑50 (MQ=255)
  cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                        >  1:3322047/1‑50 (MQ=255)
   aaTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGc                                       <  2:5080927/50‑1 (MQ=255)
   aaTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGc                                       <  2:5074389/50‑1 (MQ=255)
    aTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGcc                                      >  2:4828883/1‑50 (MQ=255)
    aTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGcc                                      >  1:5484756/1‑50 (MQ=255)
    aTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGcc                                      <  1:1480060/50‑1 (MQ=255)
    aTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGcc                                      <  1:1108957/50‑1 (MQ=255)
    aTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGcc                                      >  1:4974265/1‑50 (MQ=255)
     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                     >  1:2307802/1‑50 (MQ=255)
     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                     <  1:4133916/50‑1 (MQ=255)
     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                     <  2:5569854/50‑1 (MQ=255)
     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                     <  2:159837/50‑1 (MQ=255)
     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                     <  2:4190934/50‑1 (MQ=255)
     ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                     <  2:4330284/50‑1 (MQ=255)
       cACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCg                                    <  2:1104296/49‑1 (MQ=255)
       cACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCg                                    <  2:4667349/49‑1 (MQ=255)
       cACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGc                                   >  2:1976562/1‑50 (MQ=255)
        aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGc                                   <  1:2598654/49‑1 (MQ=255)
        aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                  >  2:1947958/1‑50 (MQ=255)
        aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                  <  2:2011569/50‑1 (MQ=255)
         ccTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                  >  1:3638953/1‑49 (MQ=255)
         ccTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTc                                 <  1:4458797/50‑1 (MQ=255)
          cTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCa                                >  2:444629/1‑50 (MQ=255)
          cTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCa                                >  2:2537769/1‑50 (MQ=255)
          cTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCa                                <  2:3076132/50‑1 (MQ=255)
           ttttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAt                               >  2:2444581/1‑50 (MQ=255)
            tttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAtt                              >  1:1641390/1‑50 (MQ=255)
            tttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAtt                              >  2:2421895/1‑50 (MQ=255)
             ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAtt                              <  1:676399/49‑1 (MQ=255)
             ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAtt                              <  2:169780/49‑1 (MQ=255)
             ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAtt                              <  2:4512843/49‑1 (MQ=255)
             ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAtt                              <  2:2283211/49‑1 (MQ=255)
             ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAt                               <  2:3754336/48‑1 (MQ=255)
             ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTg                             >  1:1070173/1‑50 (MQ=255)
             ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTg                             >  2:401799/1‑50 (MQ=255)
              tttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTg                             <  2:3512842/49‑1 (MQ=255)
              tttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTg                             <  1:3176683/49‑1 (MQ=255)
              tttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGc                            >  2:4428013/1‑50 (MQ=255)
               ttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGc                            <  1:5473075/49‑1 (MQ=255)
                tACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGcc                           <  2:3336593/49‑1 (MQ=255)
                tACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGcc                           <  1:3625035/49‑1 (MQ=255)
                tACCGTTTGACGGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGcc                           <  1:4912546/49‑1 (MQ=39)
                 aCCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCa                          <  1:23537/49‑1 (MQ=255)
                  ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAtt                        <  2:5556917/50‑1 (MQ=255)
                  ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                         >  2:2884909/1‑49 (MQ=255)
                  ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                         >  1:2792230/1‑49 (MQ=255)
                   cGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTa                       >  1:128517/1‑50 (MQ=255)
                   cGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTa                       >  2:3889393/1‑50 (MQ=255)
                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                      <  2:2479974/50‑1 (MQ=255)
                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                      >  2:733289/1‑50 (MQ=255)
                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                      <  1:3320/50‑1 (MQ=255)
                    gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                      <  1:1278096/50‑1 (MQ=255)
                     tttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAcc                     <  1:606949/50‑1 (MQ=255)
                     tttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAcc                     <  1:3119088/50‑1 (MQ=255)
                     tttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAcc                     <  1:2730027/50‑1 (MQ=255)
                     tttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                      >  1:5524386/1‑49 (MQ=255)
                     tttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                      >  1:2662153/1‑49 (MQ=255)
                       tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                      >  2:2001829/1‑47 (MQ=255)
                       tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCa                    <  1:976360/49‑1 (MQ=255)
                       tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCa                    <  1:4067510/49‑1 (MQ=255)
                        gACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAg                   <  2:5427422/49‑1 (MQ=255)
                         aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGaa                 <  2:5658953/50‑1 (MQ=255)
                         aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGaa                 >  2:4940484/1‑50 (MQ=255)
                         aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGaa                 >  1:3171315/1‑50 (MQ=255)
                         aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGaa                 >  2:628490/1‑50 (MQ=255)
                           tGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATc               >  1:155798/1‑50 (MQ=255)
                            gCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt              >  1:4482255/1‑50 (MQ=255)
                             cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt              >  1:4664946/1‑49 (MQ=255)
                             cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt              >  1:3007246/1‑49 (MQ=255)
                             cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt              >  1:1964981/1‑49 (MQ=255)
                             cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTg             <  2:2863351/50‑1 (MQ=255)
                               cGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCATAATCTGaa           >  1:4695113/1‑50 (MQ=39)
                               cGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGaa           >  2:3154285/1‑50 (MQ=255)
                                gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGa            <  2:5435730/48‑1 (MQ=255)
                                gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc          <  1:4511668/50‑1 (MQ=255)
                                gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc          >  2:3392981/1‑50 (MQ=255)
                                gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc          <  1:4725270/50‑1 (MQ=255)
                                gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc          >  1:529667/1‑50 (MQ=255)
                                gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc          >  1:5450230/1‑50 (MQ=255)
                                 cAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt         >  2:349181/1‑50 (MQ=255)
                                 cAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt         >  1:2866201/1‑50 (MQ=255)
                                 cAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt         <  1:2169697/50‑1 (MQ=255)
                                 cAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACt         >  1:2112035/1‑50 (MQ=255)
                                  aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg        >  2:3434562/1‑50 (MQ=255)
                                  aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg        >  1:936499/1‑50 (MQ=255)
                                    atcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGt      >  1:3836275/1‑50 (MQ=255)
                                     tcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGt      >  1:4791848/1‑49 (MQ=255)
                                     tcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTc     <  2:3979507/50‑1 (MQ=255)
                                       atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa   >  2:4910114/1‑50 (MQ=255)
                                       atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTca    <  2:2891223/49‑1 (MQ=255)
                                       atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTca    <  2:2209564/49‑1 (MQ=255)
                                        tcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa   >  1:1121450/1‑49 (MQ=255)
                                         cTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaac  >  1:4580150/1‑49 (MQ=255)
                                              |                                           
ACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAAC  >  NC_000913/3527714‑3527803

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: