Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 656,044 T→C 15.7% F44L (TTC→CTC)  ycbS → predicted outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE656,0440TC15.7% 172.9 / 23.8 70F44L (TTC→CTC) ycbSpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base T (32/27);  new base C (6/5);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTGGGCTGGTGGTGTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAATAAC  >  minE/655974‑656105
                                                                      |                                                             
tggtggGCTGGTGGTGTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCttttt                                                               >  1:452138/1‑71 (MQ=255)
 ggtggGCTGGTGGTGTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCtttttt                                                              >  1:1394455/1‑71 (MQ=255)
     ggCTGGTGGTGTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCaaa                                                          <  1:2635972/71‑1 (MQ=255)
     ggCTGGTGGTGTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCaaa                                                          <  1:1232245/71‑1 (MQ=255)
      gCTGGTGGTGTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAg                                                         <  1:1010234/71‑1 (MQ=255)
          gtggtgTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTCTCAAAgatga                                                     <  1:3987055/71‑1 (MQ=255)
           tggtgTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTCTCAAAgatga                                                     <  1:3991792/70‑1 (MQ=255)
                  tgtgTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTAAAAgatga                                                     >  1:1611878/1‑63 (MQ=255)
                  tgtgTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTCTCAAAGATGATCCCTCaa                                             >  1:3732593/1‑71 (MQ=255)
                    tgtCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAgat                                                       >  1:2662511/1‑59 (MQ=255)
                           aGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCg                                    <  1:3450474/71‑1 (MQ=255)
                           aGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCg                                    <  1:1808359/71‑1 (MQ=255)
                            gCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGa                                   <  1:147080/71‑1 (MQ=255)
                            gCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGa                                   <  1:353456/71‑1 (MQ=255)
                              gTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAgatgat                                                    >  1:2642447/1‑52 (MQ=255)
                              gTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAgatgat                                                    >  1:766648/1‑52 (MQ=255)
                              gTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTAAAAg                                                         >  1:658869/1‑47 (MQ=255)
                               ttCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAgatga                                                     <  1:758641/50‑1 (MQ=255)
                               ttCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAgatga                                                     <  1:2601405/50‑1 (MQ=255)
                               ttCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAgatga                                                     <  1:2547434/50‑1 (MQ=255)
                               ttCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGAttt                                >  1:199915/1‑71 (MQ=255)
                               ttCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGAttt                                >  1:531077/1‑71 (MQ=255)
                                tCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTg                               >  1:3112677/1‑71 (MQ=255)
                                tCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTg                               >  1:1079429/1‑71 (MQ=255)
                                tCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTg                               >  1:555768/1‑71 (MQ=255)
                                 cAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCaaaca                                                        >  1:381690/1‑43 (MQ=255)
                                    tGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAAtg                                           <  1:3907903/55‑1 (MQ=255)
                                    tGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAAtg                                           <  1:2720167/55‑1 (MQ=255)
                                    tGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAAtg                                           <  1:1781869/55‑1 (MQ=255)
                                    tGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAAtg                                           <  1:3325283/55‑1 (MQ=255)
                                    tGCACAAGCTGAAACCTGGTTAGATCCTGCTTTTTTCAAAGATGATCCCTCAAtg                                           <  1:2464392/55‑1 (MQ=255)
                                     gCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCg                          >  1:1785644/1‑71 (MQ=255)
                                      cacaAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCaa                                             >  1:3853179/1‑51 (MQ=255)
                                      cacaAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGt                         >  1:4076341/1‑71 (MQ=255)
                                      cacaAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGt                         >  1:3211061/1‑71 (MQ=255)
                                      cacaAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGt                         >  1:2629271/1‑71 (MQ=255)
                                      cacaAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGt                         >  1:2492239/1‑71 (MQ=255)
                                        caAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAgatg                                                      <  1:3357854/40‑1 (MQ=255)
                                        caAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAga                                                        >  1:3248628/1‑38 (MQ=255)
                                        caAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAga                                                        >  1:2900254/1‑38 (MQ=255)
                                        caAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAga                                                        >  1:359211/1‑38 (MQ=255)
                                               aaaCCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCg                                    >  1:2618825/1‑51 (MQ=255)
                                               aaaCCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCg                                    >  1:319940/1‑51 (MQ=255)
                                               aaaCCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCg                                    >  1:3017924/1‑51 (MQ=255)
                                               aaaCCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCg                                    >  1:2096777/1‑51 (MQ=255)
                                                 aCCTGGTTCGATCCTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAg               >  1:2667469/1‑70 (MQ=255)
                                                 aCCTGGTTCGATCCTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAg               >  1:641802/1‑70 (MQ=255)
                                                      gTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACaaa         <  1:3138061/71‑1 (MQ=255)
                                                      gTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACaaa         <  1:1164064/71‑1 (MQ=255)
                                                       ttCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACaaa         >  1:1172683/1‑70 (MQ=255)
                                                       ttCGATCCTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACaa          >  1:583104/1‑69 (MQ=255)
                                                          gATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAAt     <  1:2566370/71‑1 (MQ=255)
                                                          gATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAAt     <  1:985610/71‑1 (MQ=255)
                                                          gATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAAt     <  1:3244796/71‑1 (MQ=255)
                                                          gATCCTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAAt     <  1:43502/71‑1 (MQ=255)
                                                          gATCCTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAAt     <  1:86235/71‑1 (MQ=255)
                                                           aTCCTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAAt     <  1:3369900/70‑1 (MQ=255)
                                                            tCCTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAATaa   >  1:4011646/1‑71 (MQ=255)
                                                            tCCTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAATaa   >  1:2915307/1‑71 (MQ=255)
                                                             ccTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGAtt                                 >  1:1863833/1‑40 (MQ=255)
                                                             ccTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACaaaaaa      >  1:2358198/1‑67 (MQ=255)
                                                             ccTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACaaaaaa      >  1:1633171/1‑67 (MQ=255)
                                                              cTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGAttt                                <  1:3643278/40‑1 (MQ=255)
                                                              cTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGAttt                                <  1:841986/40‑1 (MQ=255)
                                                              cTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGAttt                                <  1:834533/40‑1 (MQ=255)
                                                              cTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGAttt                                <  1:629587/40‑1 (MQ=255)
                                                              cTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAATAAc  <  1:3247018/70‑1 (MQ=255)
                                                              cTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAATAAc  <  1:2997128/70‑1 (MQ=255)
                                                                 cTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACaaaaa       >  1:1178994/1‑62 (MQ=255)
                                                                 cTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACaaaaa       >  1:3436785/1‑62 (MQ=255)
                                                                 cTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACaaaaa       >  1:521555/1‑62 (MQ=255)
                                                                 cTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACaaaaa       >  1:3869410/1‑62 (MQ=255)
                                                                      |                                                             
TGGTGGGCTGGTGGTGTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAATAAC  >  minE/655974‑656105

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