Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,440,764 C→G 8.8% intergenic (‑463/‑276) ompC ← / → rcsD outer membrane porin protein C/phosphotransfer intermediate protein in two‑component regulatory system with RcsBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,440,7640CG 8.8% 176.2 / 6.2 57intergenic (‑463/‑276)ompC/rcsDouter membrane porin protein C/phosphotransfer intermediate protein in two‑component regulatory system with RcsBC
Reads supporting (aligned to +/- strand):  ref base C (30/22);  new base G (3/2);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

TTTCAACCGGATGCCTCGCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAA  >  minE/1440696‑1440830
                                                                    |                                                                  
tttCAACCGGATGCCTCGCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTAc                                                                    <  1:96187/69‑1 (MQ=255)
tttCAACCGGATGCCTCGCATTCGGTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACca                                                                  >  1:266635/1‑70 (MQ=255)
  tCAACCGGATGCCTAGCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACca                                                                  >  1:454012/1‑69 (MQ=255)
     aCCGGATGCCTCGCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGa                                                              >  1:856380/1‑70 (MQ=255)
     aCCGGATGCCTCGCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGa                                                              >  1:434247/1‑70 (MQ=255)
      ccGGATGCCTCGCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGaaa                                                            <  1:337713/71‑1 (MQ=255)
      ccGGATGCCTCGCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGaaa                                                            <  1:637002/71‑1 (MQ=255)
            gCCTCGCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGaaaaa                                                          >  1:3144014/1‑67 (MQ=255)
               tCGCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGaaa                                                            <  1:2264403/62‑1 (MQ=255)
                cGCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACg                                                  >  1:3134349/1‑71 (MQ=255)
                 gCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACcac                                                                 >  1:3460926/1‑55 (MQ=255)
                 gCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACtata                                                     >  1:1901099/1‑67 (MQ=255)
                 gCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACtata                                                     >  1:3781931/1‑67 (MQ=255)
                    ttCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGaaaa                                                           <  1:1630708/58‑1 (MQ=255)
                    ttCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTg                                              <  1:2373784/71‑1 (MQ=255)
                     tCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGt                                             >  1:2749523/1‑71 (MQ=255)
                     tCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGt                                             >  1:2736434/1‑71 (MQ=255)
                      cGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAAc                                                         >  1:545561/1‑58 (MQ=255)
                      cGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGt                                             >  1:302305/1‑70 (MQ=255)
                      cGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGt                                             <  1:1638843/70‑1 (MQ=255)
                      cGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGt                                             >  1:1365065/1‑70 (MQ=255)
                      cGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGt                                             >  1:2356259/1‑70 (MQ=255)
                       ggTTTTTTTTACCCTTCTGTACACACTTTTCATTATTCTGTGCTACCACAGAAAAAc                                                         >  1:2962809/1‑57 (MQ=255)
                        gTTTTTTTTACCCTTCTTTACACACTTTTCATTTTTCTGTGGTAGcaca                                                                <  1:3231489/49‑1 (MQ=37)
                        gTTTTTTTTACCCTTCTTTACACACTTTTCATTTTTCTGTGGTAGcaca                                                                <  1:206770/49‑1 (MQ=37)
                        gTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGtt                                            <  1:308099/69‑1 (MQ=255)
                        gTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTaa                                          >  1:1804927/1‑71 (MQ=255)
                           ttttttACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTa                                       >  1:3769413/1‑71 (MQ=255)
                                   ccTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGaaaaa                                                          <  1:2134269/44‑1 (MQ=255)
                                        tttACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAaataat                          >  1:2786229/1‑71 (MQ=255)
                                           acacacTTTTCATTATTCTGTGCTACCACAGAAAAAc                                                         <  1:2904622/37‑1 (MQ=255)
                                           acacacTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAATTATTTGTGAAATAGTTa                        <  1:917240/70‑1 (MQ=37)
                                            cacacTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAAc                      <  1:1093767/71‑1 (MQ=255)
                                            cacacTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAAc                      <  1:1154899/71‑1 (MQ=255)
                                                 ttttCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGtt                                            <  1:3229648/44‑1 (MQ=255)
                                                 ttttCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGtt                                            <  1:635696/44‑1 (MQ=255)
                                                 ttttCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGtt                                            <  1:2551208/44‑1 (MQ=255)
                                                 ttttCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGtt                                            <  1:3767507/44‑1 (MQ=255)
                                                 ttttCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAaat                             <  1:3234379/59‑1 (MQ=255)
                                                 ttttCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAAtt                         <  1:2204524/63‑1 (MQ=255)
                                                  tttCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCg                 >  1:1041159/1‑70 (MQ=255)
                                                   ttCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTa                        >  1:3175263/1‑62 (MQ=255)
                                                   ttCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTa                        >  1:2061771/1‑62 (MQ=255)
                                                   ttCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGc                <  1:1649211/70‑1 (MQ=255)
                                                   ttCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCa               >  1:978617/1‑71 (MQ=255)
                                                    tCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCa               >  1:2335574/1‑70 (MQ=255)
                                                    tCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCa               >  1:1328584/1‑70 (MQ=255)
                                                          ttCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTc        <  1:3429188/71‑1 (MQ=255)
                                                            cTGTGGTAGCACAGAATAACTATAACGCTTGTTAAc                                         >  1:553555/1‑36 (MQ=18)
                                                            cTGTGGTAGCACAGAATAACTATAACGCTTGTTAAc                                         >  1:934431/1‑36 (MQ=18)
                                                            cTGTGGTAGCACAGAATAACTATAACGCTTGTTAAc                                         >  1:958303/1‑36 (MQ=18)
                                                            cTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTcc       <  1:2727467/70‑1 (MQ=255)
                                                            cTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTcc       <  1:300053/70‑1 (MQ=255)
                                                             tgtgCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAg     >  1:2709049/1‑71 (MQ=255)
                                                                 cTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCaa  >  1:2853759/1‑70 (MQ=255)
                                                                  tACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCg                 >  1:1582627/1‑54 (MQ=255)
                                                                  tACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCg                 >  1:2983574/1‑54 (MQ=255)
                                                                   aCCACAGAAAAACTATAACGCTTGTTAACTATTTCACAaataat                          >  1:1451455/1‑44 (MQ=255)
                                                                    |                                                                  
TTTCAACCGGATGCCTCGCATTCGGTTTTTTTTACCCTTCTTTACACACTTTTCATTATTCTGTGCTACCACAGAAAAACTATAACGCTTGTTAACTATTTCACAAATAATTAACATCCGCATAATTTCCAGCAA  >  minE/1440696‑1440830

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