Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 432,031:1 +T 6.7% coding (397/1407 nt) ybfM → predicted outer membrane porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE432,0311.T 6.7% 225.5 / 13.9 60coding (397/1407 nt)ybfMpredicted outer membrane porin
Reads supporting (aligned to +/- strand):  ref base . (22/34);  new base T (2/2);  total (24/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

CGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCC  >  minE/431961‑432095
                                                                       |                                                                
cgaaATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGcc                                                                   >  1:2332484/1‑71 (MQ=255)
 gaaATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCt                                                                 <  1:2324288/71‑1 (MQ=255)
   aaTCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTg                                                                >  1:2737545/1‑70 (MQ=255)
       gcgTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATaa                                                           <  1:1235869/71‑1 (MQ=255)
          ttttCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTg                                                                <  1:2104448/63‑1 (MQ=255)
          ttttCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTg                                                                <  1:3360568/63‑1 (MQ=255)
           tttCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCt                                                       >  1:2731017/1‑71 (MQ=255)
           tttCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCt                                                       >  1:2944055/1‑71 (MQ=255)
            ttCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATaa                                                           >  1:876417/1‑66 (MQ=255)
            ttCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCt                                                       <  1:1909196/70‑1 (MQ=255)
             tCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTTTATAAAGCTg                                                      <  1:2414328/70‑1 (MQ=255)
             tCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTg                                                      <  1:3701358/70‑1 (MQ=255)
             tCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTg                                                      <  1:3277109/70‑1 (MQ=255)
             tCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTg                                                      <  1:2351553/70‑1 (MQ=255)
             tCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTg                                                      <  1:874656/70‑1 (MQ=255)
             tCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTg                                                      <  1:690346/70‑1 (MQ=255)
             tCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGc                                                     >  1:3119130/1‑71 (MQ=255)
                 aaaGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCgg                                                   >  1:214251/1‑69 (MQ=255)
                  aaGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAg                                                         >  1:2259662/1‑62 (MQ=255)
                   aGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCaa                                               <  1:438937/71‑1 (MQ=255)
                        taaAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTa                                          <  1:3587885/71‑1 (MQ=255)
                         aaaGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAg                                                         >  1:2919235/1‑55 (MQ=255)
                                aTGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTaaa                                        <  1:2424114/65‑1 (MQ=255)
                                 tGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGc                                                  <  1:504535/54‑1 (MQ=255)
                                 tGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGc                                                  <  1:402322/54‑1 (MQ=255)
                                 tGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGc                                                  <  1:3535114/54‑1 (MQ=255)
                                 tGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGc                                                  <  1:1156006/54‑1 (MQ=255)
                                 tGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGc                                                  <  1:1852606/54‑1 (MQ=255)
                                  gACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCt                                                                 >  1:574879/1‑38 (MQ=255)
                                       aGACTGGTCCGGCGACAAAAGCGGTATAAGCCTTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGtt                            <  1:335641/71‑1 (MQ=255)
                                       aGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGttt                           <  1:3059203/71‑1 (MQ=255)
                                       aGACTAGTCCGGCGACAAAAGCGGTATAAGCCTTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGtt                            <  1:1328061/71‑1 (MQ=255)
                                           tGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAAt                                             >  1:1972271/1‑49 (MQ=255)
                                           tGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGc                       >  1:3239222/1‑71 (MQ=255)
                                             gTCCGGCGACAAAAGCGGTATAAGCCTTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGc                       >  1:3295594/1‑70 (MQ=255)
                                             gTCCGGCGACAAAAGCGGTATAAGCCTTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGc                       >  1:3291091/1‑70 (MQ=255)
                                             gTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGc                       >  1:3344938/1‑69 (MQ=255)
                                             gTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGc                       >  1:3349202/1‑69 (MQ=255)
                                             gTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCg                      <  1:2533132/70‑1 (MQ=255)
                                               ccGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAgg                   >  1:2324528/1‑71 (MQ=255)
                                                cGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAggg                  <  1:2998378/71‑1 (MQ=255)
                                                 ggCGACAAAAGCGGTATAATCCTGTATAAAGCTGCGGCCAAAtt                                            <  1:2913294/44‑1 (MQ=255)
                                                 ggCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAAtt                                            <  1:3444406/44‑1 (MQ=255)
                                                 ggCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAAtt                                            <  1:1019666/44‑1 (MQ=255)
                                                 ggCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAAtt                                            <  1:735620/44‑1 (MQ=255)
                                                       aaaaGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGt                             >  1:1950797/1‑53 (MQ=255)
                                                       aaaaGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGt                             >  1:2577400/1‑53 (MQ=255)
                                                       aaaaGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGt                             >  1:2086112/1‑53 (MQ=255)
                                                       aaaaGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTa           <  1:2990425/71‑1 (MQ=255)
                                                       aaaaGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTa           >  1:3805328/1‑71 (MQ=255)
                                                         aaGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTaaa                                        <  1:1505799/40‑1 (MQ=255)
                                                         aaGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTaaa                                        <  1:2070583/40‑1 (MQ=255)
                                                          aGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACAt        >  1:297889/1‑71 (MQ=255)
                                                          aGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACAt        >  1:1613362/1‑71 (MQ=255)
                                                           gCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGTTTTGGGCGAGGGCAGGTTACAtt       <  1:2004268/71‑1 (MQ=255)
                                                           gCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACAtt       <  1:3536073/71‑1 (MQ=255)
                                                            cGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTc      <  1:1213000/71‑1 (MQ=255)
                                                                 taATCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGATTTGGCGAGGGCAGGTAACATTCAgag  >  1:1983285/1‑68 (MQ=255)
                                                                 taAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTAc          <  1:1988752/62‑1 (MQ=255)
                                                                 taAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCa     >  1:3647158/1‑67 (MQ=255)
                                                                 taAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGcc  >  1:3494661/1‑70 (MQ=255)
                                                                 taAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGcc  <  1:2122608/70‑1 (MQ=255)
                                                                 taAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGcc  <  1:2713633/70‑1 (MQ=255)
                                                                      cTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTc      >  1:2645334/1‑61 (MQ=255)
                                                                       |                                                                
CGAAATCGCGTTTTCAAAAAGTAATAAAGCCTATGACGAAGACTGGTCCGGCGACAAAAGCGGTATAAGCCTGTATAAAGCTGCGGCCAAATTTAAATACGGTCCGGTTTGGGCGAGGGCAGGTTACATTCAGCC  >  minE/431961‑432095

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: