Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,736,586 G→T 7.4% A55E (GCG→GAG)  rnc ← RNase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,736,5860GT 7.4% 228.1 / 3.5 67A55E (GCG→GAG) rncRNase III
Reads supporting (aligned to +/- strand):  ref base G (32/30);  new base T (3/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCCACGGACCAGCGTGGCGCGCATCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGC  >  minE/1736516‑1736656
                                                                      |                                                                      
gCCACGGACCAGCGTGGCGCGCATCCGGCTCATATCGCCTTCATCCACACGAGCGAAACGGTGATAAAgcg                                                                        >  1:790360/1‑71 (MQ=255)
 ccACGGACCAGCGTGGCGCGCATCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAgcgc                                                                       <  1:1033213/71‑1 (MQ=255)
   aCGGACCAGCGTGGCGCGCATCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCAt                                                                     >  1:1221253/1‑71 (MQ=255)
   aCGGACCAGCGTGGCGCGCATCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCAt                                                                     >  1:19545/1‑71 (MQ=255)
   aCGGACCAGCGTGGCGCGCATCCGGCTCATATCGCATTCATCCACACGAGGGAAACGGTGATAAAGCGCAt                                                                     >  1:1544319/1‑71 (MQ=255)
     ggACCAGCGTGGCGCGCATCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTg                                                                   <  1:1062133/71‑1 (MQ=255)
     ggACCAGCGTGGCGCGCATCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTg                                                                   >  1:208588/1‑71 (MQ=255)
     ggACCAGCGTGGCGCGCATCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTg                                                                   <  1:1369670/71‑1 (MQ=255)
     ggACCAGCGTGGCGCGCATCCGGCTCATATAGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTg                                                                   >  1:910798/1‑71 (MQ=255)
      gACCAGCGTGGCGCGCATCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTgg                                                                  >  1:885264/1‑71 (MQ=255)
       aCCAGCGTGGCGCGCATCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGc                                                                 <  1:68693/71‑1 (MQ=255)
       aCCAGCGTGGCGCGCATCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGc                                                                 <  1:494076/71‑1 (MQ=255)
       aCCAGCGTGGCGCGCATCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGc                                                                 <  1:1639107/71‑1 (MQ=255)
              tGGCGCGCATCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACg                                                          <  1:95491/71‑1 (MQ=255)
                  gcgcATCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACg                                                          >  1:146382/1‑67 (MQ=255)
                      aTCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACg                                                          <  1:2373027/63‑1 (MQ=255)
                      aTCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACg                                                          <  1:1238474/63‑1 (MQ=255)
                      aTCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAg                                                  <  1:154441/71‑1 (MQ=255)
                       tCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCg                                                                <  1:2379161/56‑1 (MQ=255)
                         cGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACg                                                          >  1:180440/1‑60 (MQ=255)
                         cGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACg                                                          >  1:271853/1‑60 (MQ=255)
                         cGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGaa                                                <  1:537613/70‑1 (MQ=255)
                         cGGCTCATATCGCCTTCATCCACACGAGGGAAACGGAGATAAAGCGCATTGGCGATaa                                                            >  1:321135/1‑58 (MQ=255)
                         cGGCTCATATCGCCTTAATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACg                                                          >  1:2552428/1‑60 (MQ=255)
                          ggCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAg                                                       >  1:310659/1‑62 (MQ=255)
                          ggCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAg                                                       >  1:2748833/1‑62 (MQ=255)
                               atatCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCTCATTGGCGATAACGTAGCTCAGaa                                                >  1:3227617/1‑64 (MQ=255)
                               atatCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATagag                                           <  1:775034/69‑1 (MQ=255)
                                tatCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCTCATTGGCGATAACg                                                          <  1:2907160/53‑1 (MQ=255)
                                 atCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATag                                             <  1:3232428/65‑1 (MQ=255)
                                  tctccTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGcc                                      >  1:764369/4‑71 (MQ=255)
                                  tCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCTCATTGGCGATAACGTAGCTCAGAATAGAGTCGc                                       >  1:1762110/1‑70 (MQ=255)
                                  tCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATa                                              <  1:2421551/63‑1 (MQ=255)
                                  tCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATa                                              <  1:2717594/63‑1 (MQ=255)
                                  tCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGcc                                      >  1:2303604/1‑71 (MQ=255)
                                      cTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTaaa                                  <  1:200772/71‑1 (MQ=255)
                                         cATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATaa                                                            >  1:2097859/1‑42 (MQ=255)
                                         cATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACg                                                          >  1:261607/1‑44 (MQ=255)
                                         cATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACg                                                          >  1:88246/1‑44 (MQ=255)
                                         cATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTaaa                                  <  1:3032948/68‑1 (MQ=255)
                                          aTCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAAt                               <  1:3145344/70‑1 (MQ=255)
                                           tCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATaga                                            <  1:2698181/56‑1 (MQ=255)
                                             cacaCGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAAt                                               >  1:2770968/1‑51 (MQ=255)
                                             cacaCGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTc                                         >  1:826446/1‑57 (MQ=255)
                                               cacGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTc                             >  1:2609329/1‑67 (MQ=255)
                                                 cGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATa                                              >  1:347714/1‑48 (MQ=255)
                                                 cGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATa                                              >  1:883178/1‑48 (MQ=255)
                                                 cGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACg                       <  1:493909/71‑1 (MQ=255)
                                                    gggAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTc                    <  1:519912/71‑1 (MQ=255)
                                                      gataCGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGt                  <  1:3200786/68‑1 (MQ=255)
                                                      gAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGt                  <  1:2870161/71‑1 (MQ=255)
                                                      gAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGt                  <  1:2403772/71‑1 (MQ=255)
                                                      gAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGt                  <  1:182821/71‑1 (MQ=255)
                                                       aaaCGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATa                                              >  1:2740081/1‑42 (MQ=255)
                                                        aaCGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTa                <  1:1600234/71‑1 (MQ=255)
                                                            gTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGtt            <  1:1009971/71‑1 (MQ=255)
                                                            gTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGtt            <  1:1079629/71‑1 (MQ=255)
                                                               aTAAAGCTCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGttt           <  1:3365777/69‑1 (MQ=255)
                                                               aTAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTAc         >  1:761650/1‑71 (MQ=255)
                                                               aTAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTAc         >  1:2079369/1‑71 (MQ=255)
                                                                tAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTc                    <  1:436218/59‑1 (MQ=255)
                                                                 aaaGCTCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTAc         >  1:1997124/1‑69 (MQ=255)
                                                                 aaaGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTaaa                         >  1:1430272/1‑53 (MQ=255)
                                                                 aaaGCGAATTGGCGATAACGTAGCTCAGAATAGAGTCg                                        >  1:1131898/1‑38 (MQ=255)
                                                                  aaGCGCATTGGCGATAACGTAGCTCAGAATAGAGTc                                         >  1:780922/1‑36 (MQ=255)
                                                                  aaGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCg                                        >  1:1242191/1‑37 (MQ=255)
                                                                  aaGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCg                                        >  1:1749556/1‑37 (MQ=255)
                                                                  aaGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCg                                        >  1:2187188/1‑37 (MQ=255)
                                                                   aGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACg                       <  1:602386/53‑1 (MQ=255)
                                                                     cgcATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTgg   <  1:1825211/71‑1 (MQ=255)
                                                                     cgcATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTgg   <  1:1536930/71‑1 (MQ=255)
                                                                      gcATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCt                            <  1:3136866/45‑1 (MQ=255)
                                                                      gcATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGc  <  1:3149470/71‑1 (MQ=255)
                                                                      gcATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGc  <  1:1292637/71‑1 (MQ=255)
                                                                      gcATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGc  <  1:117790/71‑1 (MQ=255)
                                                                      |                                                                      
GCCACGGACCAGCGTGGCGCGCATCCGGCTCATATCGCCTTCATCCACACGAGGGAAACGGTGATAAAGCGCATTGGCGATAACGTAGCTCAGAATAGAGTCGCCTAAAAATTCTAAACGCTCGTTATGTTTACTGCTGGC  >  minE/1736516‑1736656

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: