Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,747,793 T→C 14.7% I103T (ATC→ACC)  rplO → 50S ribosomal subunit protein L15

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,747,7930TC14.7% 180.4 / 21.5 75I103T (ATC→ACC) rplO50S ribosomal subunit protein L15
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base C (2/9);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 9.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

CCTGGCTAAAGTAGAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGTAACTGTTCGTGGCCTGCG  >  minE/2747728‑2747862
                                                                 |                                                                     
ccTGGCTAAAGTAGAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATc                                                                      >  1:239432/1‑67 (MQ=255)
 cTGGCTAAAGTAGAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATc                                                                      >  1:486673/1‑66 (MQ=255)
   ggCTAAAGTAGAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATc                                                                >  1:2136845/1‑70 (MQ=255)
   ggCTAAAGTAGAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATc                                                                >  1:1670129/1‑70 (MQ=255)
    gCTAAAGTAGAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGa                                                              <  1:1206183/71‑1 (MQ=255)
        aaGTAGAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGt                                                            <  1:904558/69‑1 (MQ=255)
        aaGTAGAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGt                                                            <  1:312703/69‑1 (MQ=255)
        aaGTAGAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTc                                                          >  1:118427/1‑71 (MQ=255)
        aaGTAGAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTc                                                          >  1:623214/1‑71 (MQ=255)
        aaGTAGAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTACCCAGATCGAGt                                                            <  1:303240/69‑1 (MQ=255)
        aaGTAGAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTACCCAGATCGAGt                                                            <  1:3286505/69‑1 (MQ=255)
           tAGAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGt                                                            <  1:3325084/66‑1 (MQ=255)
           tAGAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTcgcg                                                       <  1:492644/71‑1 (MQ=255)
             gAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGaa                                                     <  1:623841/71‑1 (MQ=255)
             gAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGa                                                      <  1:2154339/70‑1 (MQ=255)
              aaGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGaaa                                                    <  1:722342/71‑1 (MQ=255)
                ggCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGt                                                  >  1:1479299/1‑71 (MQ=255)
                ggCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGt                                                  <  1:133117/71‑1 (MQ=255)
                ggCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTACCCa                                                                    <  1:42025/53‑1 (MQ=255)
                ggCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTACCCa                                                                    <  1:1967081/53‑1 (MQ=255)
                 tcGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATc                                                                <  1:294420/55‑1 (MQ=255)
                 gCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTACCCAGATCGAGTTCGCGAAAGt                                                  >  1:298118/1‑70 (MQ=255)
                   ggTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGAt                                               <  1:3933748/71‑1 (MQ=255)
                     tgtAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGa                                                      <  1:3119482/62‑1 (MQ=255)
                     tgtAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTACCCAGATCGa                                                              <  1:2935188/54‑1 (MQ=255)
                       tagtagACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGaaa                                                    >  1:2253276/1‑62 (MQ=255)
                       tagtagACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCctg                                           >  1:374336/1‑71 (MQ=255)
                         gtagACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCctggc                                         <  1:2826717/71‑1 (MQ=255)
                         gtagACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCctggc                                         <  1:215569/71‑1 (MQ=255)
                          tagACCTGAACACGCTGAAAGCGGCTAACATTATCGGTAt                                                                       <  1:2869646/40‑1 (MQ=255)
                          tagACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCGAGATCGAGTTcg                                                         <  1:221649/54‑1 (MQ=255)
                          tagACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCctggc                                         <  1:2762356/70‑1 (MQ=255)
                              ccTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCctggctgg                                      >  1:877494/1‑69 (MQ=255)
                              ccTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCtggctggc                                     >  1:615620/1‑70 (MQ=255)
                              ccTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCtggctggc                                     >  1:3792270/1‑70 (MQ=255)
                              ccTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCtggctggc                                     >  1:3055467/1‑70 (MQ=255)
                              ccTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCg                                    <  1:3470962/71‑1 (MQ=255)
                              ccTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCg                                    <  1:1714685/71‑1 (MQ=255)
                              ccTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCg                                    <  1:1869568/71‑1 (MQ=255)
                               cTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGAt                                                                 >  1:3217629/1‑41 (MQ=255)
                                tGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTcgcg                                                       <  1:1333043/50‑1 (MQ=255)
                                 gAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGa                                                      <  1:1160736/50‑1 (MQ=255)
                                 gAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAg                                 >  1:1371193/1‑71 (MQ=255)
                                  aaCACGCTGAAAGCGGCTAACATTATCGGTACCCAGATCGAGTTCGCGaaa                                                    >  1:3726816/1‑51 (MQ=255)
                                   acacGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTa                               <  1:1648609/71‑1 (MQ=255)
                                    cacGCTGAAAGCGGCTAACATTATCGGTATCCAGATc                                                                >  1:2987342/1‑37 (MQ=255)
                                    cacGCTGAAAGCGGCTAACATTATCGGTATCCAGATc                                                                >  1:3715212/1‑37 (MQ=255)
                                     acGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAAc                             <  1:1796750/71‑1 (MQ=255)
                                     acGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAAc                             <  1:1463025/71‑1 (MQ=255)
                                       gCTGAAAGCGGCTAACATTATCGGTACCCAGATCGAGtt                                                           <  1:1092495/39‑1 (MQ=255)
                                        cTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGaaa                                                    >  1:1458812/1‑45 (MQ=255)
                                        cTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGaaa                                                    >  1:3407222/1‑45 (MQ=255)
                                        cTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGaaa                                                    >  1:2078544/1‑45 (MQ=255)
                                        cTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGaa                                  >  1:2860689/1‑63 (MQ=255)
                                          gAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGAc                          >  1:2976228/1‑69 (MQ=255)
                                           aaaGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGt                                <  1:1603565/62‑1 (MQ=255)
                                            aaGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCg                      <  1:1201713/71‑1 (MQ=255)
                                             aGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTaa                              >  1:1229657/1‑62 (MQ=255)
                                              gCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGt                    <  1:2949329/71‑1 (MQ=255)
                                              gCGGCTAACATTATCGGTACCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGt                    <  1:3013451/71‑1 (MQ=255)
                                                ggCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGTaa                  >  1:3465994/1‑71 (MQ=255)
                                                ggCTAACATTATCGGTACCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGTaa                  <  1:519194/71‑1 (MQ=255)
                                                 gCTAACATTATCGGTACCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGTaa                  <  1:1777567/70‑1 (MQ=255)
                                                   tAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGTaa                  >  1:2669473/1‑68 (MQ=255)
                                                    aaCATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGTAACTGt              >  1:1096887/1‑71 (MQ=255)
                                                    aaCATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGTAACTGt              >  1:178238/1‑71 (MQ=255)
                                                      cATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATc                                              <  1:3926036/37‑1 (MQ=255)
                                                      cATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGTAACTGTTc            <  1:2406771/71‑1 (MQ=255)
                                                       aTTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCgg                     >  1:1115998/1‑61 (MQ=255)
                                                        ttATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGTAACTGTTCg           <  1:2723071/70‑1 (MQ=255)
                                                        ttATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGTAACTGTTCGt          <  1:2668757/71‑1 (MQ=255)
                                                         tATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAAc                             >  1:570453/1‑51 (MQ=255)
                                                              gTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGTAACTGt              <  1:2542835/61‑1 (MQ=255)
                                                                aTCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGTAACTGTTCGTGGCCTGc   >  1:1935109/1‑70 (MQ=255)
                                                                aTCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGTAACTGTTCGTGGCCTGCg  >  1:3309554/1‑71 (MQ=255)
                                                                aTCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGTAACTGTTCGTGGCCTGCg  >  1:1772614/1‑71 (MQ=255)
                                                                 |                                                                     
CCTGGCTAAAGTAGAAGGCGGTGTAGTAGACCTGAACACGCTGAAAGCGGCTAACATTATCGGTATCCAGATCGAGTTCGCGAAAGTGATCCTGGCTGGCGAAGTAACGACTCCGGTAACTGTTCGTGGCCTGCG  >  minE/2747728‑2747862

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: