Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 656,044 T→C 26.5% 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,0440TC26.5% 77.9 / 33.0 49F44L (TTC→CTC) ycbSpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base T (20/16);  new base C (7/6);  total (27/22)
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

GTGGGCTGGTGGTGTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAATAACG  >  minE/655976‑656106
                                                                    |                                                              
gtggGCTGGTGGTGTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCtttttt                                                               <  1:1051245/70‑1 (MQ=255)
 tggGCTGGTGGTGTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCa                                                             <  1:2096009/71‑1 (MQ=255)
 tggGCTGGTGGTGTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCa                                                             <  1:2114898/71‑1 (MQ=255)
   ggCTGGTGGTGTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCaaa                                                           <  1:962655/71‑1 (MQ=255)
   ggCTGGTGGTGTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCaa                                                            <  1:528101/70‑1 (MQ=255)
     cTGGTGGTGTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAg                                                          <  1:2029191/70‑1 (MQ=255)
        gtggtgTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAgatga                                                      <  1:911895/71‑1 (MQ=255)
           gtgtTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATc                                                    >  1:555091/1‑70 (MQ=255)
              tttGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTc                                                              <  1:1193308/57‑1 (MQ=255)
              tttGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTc                                                <  1:138402/71‑1 (MQ=255)
              tttGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTc                                                <  1:741583/71‑1 (MQ=255)
                tgtgTCGGGAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCaa                                              >  1:1265127/1‑71 (MQ=255)
                            gTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAAtg                                            >  1:1200930/1‑61 (MQ=255)
                            gTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAAtg                                            >  1:1408396/1‑61 (MQ=255)
                            gTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCACAAt                                             >  1:575675/1‑60 (MQ=255)
                             ttCAATGCATAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAg                                                          >  1:1059754/1‑46 (MQ=255)
                             ttCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCaaa                                                           <  1:1149640/45‑1 (MQ=255)
                             ttCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCaaa                                                           <  1:195678/45‑1 (MQ=255)
                             ttCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAg                                                          >  1:1858299/1‑46 (MQ=255)
                             ttCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGcc                                      <  1:1936620/66‑1 (MQ=255)
                               cAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTg                                >  1:203272/1‑70 (MQ=255)
                                  tGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAAt                                             <  1:630401/54‑1 (MQ=255)
                                    cacaAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAgatgat                                                     >  1:1710497/1‑44 (MQ=255)
                                    cacaAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAgatgat                                                     >  1:1693029/1‑44 (MQ=255)
                                    cacaAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAgatga                                                      >  1:232718/1‑43 (MQ=255)
                                    cacaAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGAttt                                 >  1:83702/1‑64 (MQ=255)
                                    cacaAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGAttt                                 >  1:831984/1‑64 (MQ=255)
                                    cacaAGCTGAAACCTGGTTCGATCCTGCTTTTCTCaaa                                                           >  1:1782643/1‑38 (MQ=255)
                                    cacaAGCTGAAACCTGGTTCGATCCTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGt                          >  1:1018730/1‑71 (MQ=255)
                                      caAGCTGAAACCTGGTTCGATCCTGCTTTTCTCAAAgat                                                        <  1:839206/39‑1 (MQ=255)
                                       aaGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTc                       <  1:1484130/71‑1 (MQ=255)
                                         gCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAgatg                                                       >  1:1602541/1‑37 (MQ=255)
                                          cTGAAACCTGGTTCGATCCTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGaa                    <  1:1577307/71‑1 (MQ=255)
                                            gAAACCTGGTTCGATCCTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGaa                    <  1:987035/69‑1 (MQ=255)
                                             aaaCCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGaaa                   <  1:533672/69‑1 (MQ=255)
                                             aaaCCTGGTTCGATCCTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGaaaa                  <  1:1099017/70‑1 (MQ=255)
                                             aaaCCTGGTTCGATCCTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTACGaaaa                  <  1:284618/70‑1 (MQ=255)
                                              aaCCTGGTTCGATCCTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAg                >  1:1052816/1‑71 (MQ=255)
                                                 cTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGAc             <  1:199145/71‑1 (MQ=255)
                                                   ggTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACa            >  1:1952688/1‑70 (MQ=255)
                                                     ttCGATCCTGCTTTTTTCAAAGATGATCCCTCAAtgg                                           >  1:75081/1‑37 (MQ=255)
                                                     ttCGATCCTGCTTTTTTCAAAGATGATCCCTCAAtgg                                           >  1:1234603/1‑37 (MQ=255)
                                                        gATCCTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAAt      <  1:272917/71‑1 (MQ=255)
                                                           ccTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAATAAc   >  1:1162064/1‑71 (MQ=255)
                                                           ccTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAATAAc   >  1:1474754/1‑71 (MQ=255)
                                                            cTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACaaaaaa       >  1:983626/1‑66 (MQ=255)
                                                            cTGCTTTTCTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAATAACg  >  1:938523/1‑71 (MQ=255)
                                                            cTGCTTTTCTCAAAGATGATCACTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAATAACg  >  1:595448/1‑71 (MQ=255)
                                                              gCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAATaa    >  1:1020671/1‑67 (MQ=255)
                                                               cTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAATAACg  >  1:2086327/1‑68 (MQ=255)
                                                                    |                                                              
GTGGGCTGGTGGTGTTTGTGTCGGCAGCGTTCAATGCACAAGCTGAAACCTGGTTCGATCCTGCTTTTTTCAAAGATGATCCCTCAATGGTGGCCGATTTGTCTCGTTTCGAAAAAGGACAAAAAATAACG  >  minE/655976‑656106

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