Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,211,683 A→G 18.6% I465I (ATT→ATC rng ← ribonuclease G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,211,6830AG18.6% 91.8 / 17.3 43I465I (ATT→ATCrngribonuclease G
Reads supporting (aligned to +/- strand):  ref base A (20/15);  new base G (4/4);  total (24/19)
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

TCATTACGACGTCAAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCATA  >  minE/2211613‑2211751
                                                                      |                                                                    
tcatTACGACGTCAAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGaaa                                                                      <  1:404652/71‑1 (MQ=255)
tcatTACGACGTCAAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGaaa                                                                      <  1:450382/71‑1 (MQ=255)
tcatTACGACGTCAAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGaaa                                                                      <  1:1351598/71‑1 (MQ=255)
 catTACGACGTCAAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGaaaa                                                                     >  1:517866/1‑71 (MQ=255)
  atTACGACGTCAAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTGCTTTAACCTGTTTGCCAACGAAAAt                                                                    >  1:115778/1‑71 (MQ=255)
  atTACGACGTCAAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAAt                                                                    >  1:1462552/1‑71 (MQ=255)
             aaaCTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCAc                                                              >  1:197548/1‑66 (MQ=255)
             aaaCTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTcc                                                          >  1:1247890/1‑70 (MQ=255)
                cTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTc                                                                 >  1:1214750/1‑60 (MQ=255)
                      cTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGaaa                                                                      <  1:91193/49‑1 (MQ=255)
                      cTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAgt                                                >  1:1183196/1‑71 (MQ=255)
                       tGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAgtg                                               >  1:356843/1‑71 (MQ=255)
                        ggTTATAGAGCGGTTCAATTTTTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAgtgt                                              >  1:115365/1‑71 (MQ=255)
                         gTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAgtgt                                              <  1:1535818/70‑1 (MQ=255)
                         gTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAGtgtg                                             <  1:1443985/71‑1 (MQ=255)
                         gTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAgtgt                                              <  1:833834/70‑1 (MQ=255)
                          ttATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAGAtt                                                                   <  1:1027135/48‑1 (MQ=255)
                            ataGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGaaaa                                                                     >  1:1391474/1‑44 (MQ=255)
                            ataGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCa                                                      >  1:919867/1‑59 (MQ=255)
                              agagCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCAc                                                              >  1:1377906/1‑49 (MQ=255)
                               gagCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCtt                                       <  1:867178/71‑1 (MQ=255)
                                      tCAATTTGTACTTTAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAGTGTGACACTTCGCCttt                                >  1:1186360/1‑71 (MQ=255)
                                       cAATTTGTACTTTAACCTGTTTGCCAACGAAAAttt                                                                  >  1:1423090/1‑36 (MQ=255)
                                        aaTTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTc                               <  1:835020/70‑1 (MQ=255)
                                         aTTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGa                                                  >  1:648828/1‑50 (MQ=255)
                                          tttGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTcc                                                          <  1:1935384/41‑1 (MQ=255)
                                          tttGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAgtgt                                              <  1:1161618/53‑1 (MQ=255)
                                           ttGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAgg                           >  1:1329379/1‑71 (MQ=255)
                                                 tttAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGac                                                 <  1:1835361/43‑2 (MQ=255)
                                                   tAACCTGTTTGCCAACGAAGATTTCCACTTCCGCCAGCGAgtg                                               >  1:973805/1‑43 (MQ=255)
                                                   tAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAgtg                                               >  1:1557005/1‑43 (MQ=255)
                                                     aCCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGc                    <  1:2037415/68‑1 (MQ=255)
                                                     aCCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTAc                 >  1:161247/1‑71 (MQ=255)
                                                       cTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAgtgt                                              <  1:923377/40‑1 (MQ=255)
                                                       cTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAgtgt                                              <  1:2016227/40‑1 (MQ=255)
                                                         gTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGa                                            >  1:2133015/1‑40 (MQ=255)
                                                         gTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGa                                            >  1:457207/1‑40 (MQ=255)
                                                          tttGCCAACGAAAATTTCCACTTCCGCCAGCGAGtgtg                                             <  1:1994825/38‑1 (MQ=255)
                                                          tttGCCAACGAAAATTTCCACTTCCGCCAGCGAGtgtg                                             <  1:78434/38‑1 (MQ=255)
                                                          tttGCCAACGAAAATTTCCACTTCCGCCAGCGAGtgtg                                             <  1:1217101/38‑1 (MQ=255)
                                                          tttGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCg            <  1:741098/71‑1 (MQ=255)
                                                             gCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGga          <  1:1687427/70‑1 (MQ=255)
                                                              ccAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGga          <  1:2021039/69‑1 (MQ=255)
                                                                aaCGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAg      >  1:2054573/1‑71 (MQ=255)
                                                                aaCGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAg      >  1:1251650/1‑71 (MQ=255)
                                                                 aCGAAGATTTCCACTTCAGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGc     <  1:1949891/71‑1 (MQ=255)
                                                                   gAAGATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCa    >  1:1394016/1‑70 (MQ=255)
                                                                   gAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCa    >  1:123711/1‑70 (MQ=255)
                                                                   gAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCa    >  1:1678930/1‑70 (MQ=255)
                                                                    aaaaTTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCATa  <  1:681289/71‑1 (MQ=255)
                                                                      |                                                                    
TCATTACGACGTCAAACTGCTCCTGGTTATAGAGCGGTTCAATTTGTACTTTAACCTGTTTGCCAACGAAAATTTCCACTTCCGCCAGCGAGTGTGACTCTTCGCCTTTCAAGGCTTCAGCTACTGCCGGAGAAGCATA  >  minE/2211613‑2211751

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