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,7600AG85.8% 530.5 / 67.2 183N431S (AAC→AGC) igaAinner membrane protein ‑ inhibits the Rcs signaling pathway
Reads supporting (aligned to +/- strand):  ref base A (10/16);  new base G (99/58);  total (109/74)
Fisher's exact test for biased strand distribution p-value = 2.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 203 total aligned reads displayed.

AACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAACA  >  NC_000913/3527713‑3527804
                                               |                                            
 aCCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcga                                           >  2:2531315/1‑50 (MQ=255)
 aCCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAacga                                           <  2:5692212/50‑1 (MQ=255)
  ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                          >  2:5123977/1‑50 (MQ=255)
  ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                          >  1:1622374/1‑50 (MQ=255)
  ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                          >  2:2657960/1‑50 (MQ=255)
  ccAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgac                                          <  1:2553420/50‑1 (MQ=255)
   cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                         >  2:4733891/1‑50 (MQ=255)
   cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                         >  2:3780760/1‑50 (MQ=255)
   cAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAgcgacg                                         <  2:2137644/50‑1 (MQ=255)
    aaTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGc                                        <  1:5344366/50‑1 (MQ=255)
     aTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGcc                                       >  1:2217847/1‑50 (MQ=255)
     aTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGcc                                       >  1:1988413/1‑50 (MQ=255)
     aTTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGcc                                       <  2:4111535/50‑1 (MQ=255)
      ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      <  2:2333709/50‑1 (MQ=255)
      ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      <  1:2496213/50‑1 (MQ=255)
      ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      <  1:1464894/50‑1 (MQ=255)
      ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      <  2:5768942/50‑1 (MQ=255)
      ttCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGccc                                      <  2:1791032/50‑1 (MQ=255)
      gtCACCTTTTTTTCCGTTTGACTGCTCGCAGATCATCTGGAGCGACGccc                                      <  1:5755310/49‑1 (MQ=255)
         aCCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCt                                   >  1:2031652/1‑50 (MQ=255)
          ccTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTc                                  <  1:533273/50‑1 (MQ=255)
          ccTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTc                                  <  2:1398004/50‑1 (MQ=255)
          ccTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTc                                  <  2:3962576/50‑1 (MQ=255)
           cTTTTTTGCCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCa                                 <  2:1748926/50‑1 (MQ=255)
           cTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCa                                 >  2:5623990/1‑50 (MQ=255)
            ttttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAt                                >  1:3325857/1‑50 (MQ=255)
            ttttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAt                                >  1:5601222/1‑50 (MQ=255)
             tttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAtt                               >  1:1074575/1‑50 (MQ=255)
             tttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAtt                               >  2:4501709/1‑50 (MQ=255)
             tttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAtt                               >  2:2995833/1‑50 (MQ=255)
              ttttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCAtt                               <  1:3808889/49‑1 (MQ=255)
              ttttACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTg                              >  1:5324135/1‑50 (MQ=255)
               tttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTg                              <  2:416160/49‑1 (MQ=255)
               tttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGc                             >  2:2059587/1‑50 (MQ=255)
               tttACCGTTTGACTGCTCGCAGATCATCTGGAACTACGCCCGCTCATTg                              <  1:2276363/49‑1 (MQ=255)
                ttACCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGc                             >  2:1643791/1‑49 (MQ=255)
                  aCCGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                          >  2:4352553/1‑50 (MQ=255)
                   ccGTTTGACTGCTCGCAGATCCTCTGGAGCGACGCCCGCTCATTGCCAt                          >  2:5079094/1‑49 (MQ=38)
                   ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                          >  1:4290267/1‑49 (MQ=255)
                   ccGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCAt                          <  1:3185490/49‑1 (MQ=255)
                   ccGTTTGACTGCTCGCAGATCATCTGGAGAGACGCCCGCTCATTGCCAt                          >  2:1461203/1‑49 (MQ=38)
                    cGTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTa                        >  1:2616140/1‑50 (MQ=255)
                     gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       <  2:4367529/50‑1 (MQ=255)
                     gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       <  2:211927/50‑1 (MQ=255)
                     gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:4906559/1‑50 (MQ=255)
                     gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  2:10299/1‑50 (MQ=255)
                     gTTTGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATCGCCATTAc                       <  2:61922/50‑1 (MQ=255)
                       ttGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:3483907/1‑48 (MQ=255)
                        tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:1888491/1‑47 (MQ=255)
                        tGACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTAc                       >  1:4589554/1‑47 (MQ=255)
                        tGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTAc                       >  2:1551511/1‑47 (MQ=255)
                        tGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCa                     <  2:4135642/49‑1 (MQ=255)
                         gACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAg                    >  1:1253820/1‑49 (MQ=255)
                         gACTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAg                    >  1:5384673/1‑49 (MQ=255)
                         gACTGCTCGCAGATCATCTGGAGCGACGCCCGCCCATTGCCATTACCAg                    >  2:3545467/1‑49 (MQ=255)
                          aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGaa                  >  1:4707814/1‑50 (MQ=255)
                          aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGa                   <  2:5554878/49‑1 (MQ=255)
                          aCTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGa                   <  1:1939640/49‑1 (MQ=255)
                          aCTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGa                   <  2:3189096/49‑1 (MQ=255)
                           cTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAAt                 <  1:4345724/50‑1 (MQ=255)
                           cTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAAt                 <  2:2394858/50‑1 (MQ=255)
                           cTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAAt                 <  1:288182/50‑1 (MQ=255)
                           cTGCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAAt                 >  2:2819881/1‑50 (MQ=255)
                           cTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAAt                 >  1:3008637/1‑50 (MQ=255)
                            tacTCGCAGTACATCTGGAGCGACGCCCTCTCATTGCCATTACCAGAATc                >  2:4881476/3‑50 (MQ=255)
                             gCTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt               >  1:4190575/1‑50 (MQ=255)
                              cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATc                <  2:710706/48‑1 (MQ=255)
                              cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt               <  2:1211553/49‑1 (MQ=255)
                              cTCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCt               <  2:2726628/49‑1 (MQ=255)
                               tCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTg              >  1:2749042/1‑49 (MQ=255)
                               tCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTg              >  1:561848/1‑49 (MQ=255)
                               tCGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTg              >  1:5810277/1‑49 (MQ=255)
                                cGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGaa            >  2:2202874/1‑50 (MQ=255)
                                cGCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGaa            >  2:3300695/1‑50 (MQ=255)
                                 gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc           <  2:5473340/50‑1 (MQ=255)
                                 gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc           >  1:335430/1‑50 (MQ=255)
                                 gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc           >  1:3721291/1‑50 (MQ=255)
                                 gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc           <  1:4779062/50‑1 (MQ=255)
                                 gCAGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAAc           <  1:5460266/50‑1 (MQ=255)
                                 gCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGa             <  2:1867790/48‑1 (MQ=255)
                                  cAGATCATCTGGAGTGACGCCCGCTCATTGCCATTACCAGAATCTGAACt          >  2:5567847/1‑50 (MQ=39)
                                   aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg         >  1:4976978/1‑50 (MQ=255)
                                   aGATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTg         >  1:5167533/1‑50 (MQ=255)
                                    gATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTgg        >  1:4000093/1‑50 (MQ=255)
                                    gATCATCTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTgg        <  2:5075304/50‑1 (MQ=255)
                                    gATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTgg        <  2:800072/50‑1 (MQ=255)
                                     atcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGt       <  1:451136/50‑1 (MQ=255)
                                     atcatcTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGt       <  1:1571904/50‑1 (MQ=255)
                                      tcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTc      <  2:1078943/50‑1 (MQ=255)
                                      tcatcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTc      <  1:3219506/50‑1 (MQ=255)
                                       catcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTc      >  2:2859843/1‑49 (MQ=255)
                                       catcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTc      >  2:5232091/1‑49 (MQ=255)
                                        atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa    <  1:4182998/50‑1 (MQ=255)
                                        atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTca     <  2:4037182/49‑1 (MQ=255)
                                        atcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTca     <  1:2451303/49‑1 (MQ=255)
                                         tcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa    >  2:43969/1‑49 (MQ=255)
                                         tcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa    >  2:4276734/1‑49 (MQ=255)
                                         tcTGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaa    >  1:5236899/1‑49 (MQ=255)
                                           tGGAGCGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaaca  >  1:1863758/1‑49 (MQ=255)
                                             gAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTcaaca  <  1:913950/47‑1 (MQ=255)
                                               |                                            
AACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATTACCAGAATCTGAACTGGTCAACA  >  NC_000913/3527713‑3527804

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