Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,440,185 G→A 7.6% D39D (GAC→GAT ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,440,1850GA 7.6% 283.4 / 8.2 92D39D (GAC→GATompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (45/40);  new base A (5/2);  total (50/42)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

TGAAGCCAAGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCG  >  minE/1440118‑1440251
                                                                   |                                                                  
tGAAGCCAAGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTc                                                                  >  1:225469/1‑70 (MQ=255)
  aaGCCAAGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTa                                                                <  1:4625807/70‑1 (MQ=255)
  aaGCCAAGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTa                                                                <  1:4991474/70‑1 (MQ=255)
    gCCAAGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACtt                                                             <  1:3949881/71‑1 (MQ=255)
     ccAAGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACttt                                                            <  1:1641231/71‑1 (MQ=255)
       aaGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTa                                                           >  1:3048947/1‑70 (MQ=255)
       aaGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTAc                                                          >  1:5072382/1‑71 (MQ=255)
        aGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGt                                                                   >  1:3094189/1‑61 (MQ=255)
        aGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTc                                                                  <  1:1982002/62‑1 (MQ=255)
        aGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTc                                                                  <  1:3971222/62‑1 (MQ=255)
        aGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTc                                                                  <  1:5561920/62‑1 (MQ=255)
         gACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCg                                                        >  1:4437036/1‑71 (MQ=255)
          aCGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCg                                                        >  1:4789128/1‑70 (MQ=255)
           cGCATGTATGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTa                                                      <  1:4228482/71‑1 (MQ=255)
             cATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACa                                                    <  1:5384074/71‑1 (MQ=255)
               tGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTc                                                                  <  1:4241986/55‑1 (MQ=255)
               tGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTa                                                                >  1:4168794/1‑57 (MQ=255)
               tGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTa                                                                >  1:2733721/1‑57 (MQ=255)
               tGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTa                                                                >  1:834612/1‑57 (MQ=255)
               tGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCg                                                        <  1:4823278/65‑1 (MQ=255)
               tGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGa                                                  <  1:1611298/71‑1 (MQ=255)
                 tAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATc                                                >  1:371044/1‑71 (MQ=255)
                 tAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATc                                                >  1:757451/1‑71 (MQ=255)
                   ggTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTAc                                                          >  1:243293/1‑59 (MQ=255)
                   ggTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTAc                                                          >  1:595585/1‑59 (MQ=255)
                     tCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAAt                                            <  1:106870/71‑1 (MQ=255)
                      cTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGTCCGTCTACTTTACCGTACAGATCTAAtt                                           >  1:46530/1‑71 (MQ=255)
                         gTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCGGTCTACTTTACCGTACAGATCTAAttt                                          <  1:2837429/69‑1 (MQ=255)
                         gTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAAttt                                          <  1:5245665/69‑1 (MQ=255)
                         gTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgt                                        <  1:4054538/71‑1 (MQ=255)
                          tCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgtt                                       <  1:3972393/71‑1 (MQ=255)
                          tCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgtt                                       <  1:4415298/71‑1 (MQ=255)
                             ccATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGcc                                    >  1:3906280/1‑71 (MQ=255)
                             ccATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGcc                                    <  1:624218/71‑1 (MQ=255)
                              cATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACtt                                                             >  1:4776643/1‑45 (MQ=255)
                              cATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCg                                   >  1:689984/1‑71 (MQ=255)
                                tCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAg                                                   <  1:5521815/53‑1 (MQ=255)
                                tCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAg                                                   <  1:5394859/53‑1 (MQ=255)
                                tCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAg                                                   <  1:5167924/53‑1 (MQ=255)
                                tCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCAGTCTACTTTACCGTACAGATCTAAtt                                           <  1:3704217/61‑1 (MQ=255)
                                 cTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCt                                >  1:3663035/1‑71 (MQ=255)
                                 cTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCt                                >  1:2346539/1‑71 (MQ=255)
                                  tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACtt                                                             >  1:4413362/1‑41 (MQ=255)
                                  tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTa                                                      >  1:4845125/1‑48 (MQ=255)
                                  tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTa                                                      >  1:78174/1‑48 (MQ=255)
                                  tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                               >  1:3347306/1‑71 (MQ=255)
                                  tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                               >  1:1934871/1‑71 (MQ=255)
                                  tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                               >  1:2625383/1‑71 (MQ=255)
                                  tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                               >  1:2906427/1‑71 (MQ=255)
                                  tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                               >  1:4393193/1‑71 (MQ=255)
                                  tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                               >  1:431957/1‑71 (MQ=255)
                                    cATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCTTCTttg                             <  1:180053/71‑1 (MQ=255)
                                    cATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCt                                <  1:2752969/68‑1 (MQ=255)
                                    cATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTttg                             <  1:3816773/71‑1 (MQ=255)
                                    cATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTttg                             <  1:696502/71‑1 (MQ=255)
                                     aTCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCg                                                        <  1:4857748/43‑1 (MQ=255)
                                     aTCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCg                                                        <  1:914361/43‑1 (MQ=255)
                                      tCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAAt                                            <  1:780484/54‑1 (MQ=255)
                                      tCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgt                                        >  1:4338540/1‑58 (MQ=255)
                                       cTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTAc                                                          >  1:5208711/1‑39 (MQ=255)
                                         ttgttgTCAGAGAAATAGTGCAGGCCGTCTACttt                                                            <  1:5200954/35‑1 (MQ=255)
                                         ttgttgTCAGAGAAATAGTGCAGGCCGTCTACttt                                                            <  1:1651955/35‑1 (MQ=255)
                                          tgttgtCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTa                        <  1:4519547/70‑1 (MQ=255)
                                          tgttgtCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTa                        >  1:2978238/1‑70 (MQ=255)
                                           gttgtCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTaa                       <  1:3901941/70‑1 (MQ=255)
                                            ttatCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATaa                       >  1:4584766/4‑69 (MQ=2)
                                            ttatCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAAc                     >  1:5075339/4‑71 (MQ=9)
                                             tgtCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTttgtt                           >  1:5530487/1‑64 (MQ=255)
                                              atCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATttgt                                        >  1:1072543/2‑50 (MQ=2)
                                               tCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACtt                   <  1:2968300/70‑1 (MQ=255)
                                               tCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACtt                   >  1:323894/1‑70 (MQ=255)
                                                 agagAAATAGTGCAGGCCGTCTACTTTACCGTACa                                                    >  1:3062996/1‑35 (MQ=255)
                                                 agagAAATAGTGCAGGCCGTCTACTTTACCGTACa                                                    >  1:3978740/1‑35 (MQ=255)
                                                 agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAAt                                            <  1:4981965/43‑1 (MQ=255)
                                                 agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgt                                        >  1:63663/1‑47 (MQ=255)
                                                 agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgt                                        >  1:1637962/1‑47 (MQ=255)
                                                 agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTcag                >  1:2106970/1‑71 (MQ=255)
                                                 agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTcag                >  1:4428153/1‑71 (MQ=255)
                                                   aaaaaTAATGCAGGCCATCAACTTTGCCATACAGATCTAAttt                                          <  1:3358368/41‑1 (MQ=2)
                                                    gaAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTcagcag             >  1:1661661/1‑71 (MQ=255)
                                                    aaaaTAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGt                                  >  1:5335882/2‑50 (MQ=14)
                                                    aaaaTAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTcggctg             >  1:1430497/2‑69 (MQ=9)
                                                     aaaTAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTc                  >  1:3737092/1‑65 (MQ=255)
                                                         aGTGCAGGCCGTCTACTTTTCCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGttt        >  1:1521748/1‑71 (MQ=255)
                                                         aGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTca                 >  1:4284700/1‑62 (MQ=255)
                                                          gTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTaaa                      >  1:2672177/1‑56 (MQ=255)
                                                          gTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtg       <  1:3563717/71‑1 (MQ=255)
                                                          gTGCAGGCCGTCTACTTTAACGTACAGATCTAATTTGTTGCCg                                   >  1:1951569/1‑43 (MQ=255)
                                                             cAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTca                 <  1:4235278/58‑1 (MQ=255)
                                                              aGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgctgc   <  1:482148/71‑1 (MQ=255)
                                                               ggCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTgcg  <  1:767531/71‑1 (MQ=255)
                                                               ggCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGtgcagca  <  1:2463377/71‑2 (MQ=12)
                                                                   gTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                               <  1:1041773/38‑1 (MQ=255)
                                                                   |                                                                  
TGAAGCCAAGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCG  >  minE/1440118‑1440251

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: