Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,663,990 G→A 17.6% G325E (GGG→GAG)  rtcB → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,663,9900GA17.6% 114.0 / 17.4 51G325E (GGG→GAG) rtcBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (15/27);  new base A (6/3);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

TGTCTGCGCGTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCGAACTAA  >  minE/2663921‑2664056
                                                                     |                                                                  
tgtCTGCGCGTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTggg                                                                   >  1:3065008/1‑71 (MQ=255)
 gtCTGCGCGTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGc                                                                  >  1:195513/1‑71 (MQ=255)
  tCTGCGCGTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGc                                                                  <  1:983888/70‑1 (MQ=255)
      cgcgTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTg                                                                <  1:1476652/68‑1 (MQ=255)
        cgTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGaaa                                                           <  1:2232928/71‑1 (MQ=255)
          tGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATg                                                         <  1:3204735/71‑1 (MQ=255)
             tGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGAGCTGGGAAATg                                                         >  1:437851/1‑68 (MQ=255)
              ggTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGaa                                                       >  1:1871514/1‑69 (MQ=255)
               gTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGAGCTGGGAAATGAag                                                      <  1:2736640/69‑1 (MQ=255)
                 cAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTc                                                  >  1:3940286/1‑71 (MQ=255)
                  aaTATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTTGGCTGGGAAATGAAGAGTCg                                                 <  1:1027827/71‑1 (MQ=255)
                  aaTATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCg                                                 >  1:40308/1‑71 (MQ=255)
                    tatGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGtt                                               <  1:2214856/71‑1 (MQ=255)
                      tGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGtt                                               >  1:2191142/1‑69 (MQ=255)
                      tGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCt                                             <  1:3360339/71‑1 (MQ=255)
                      tGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCt                                             <  1:3007106/71‑1 (MQ=255)
                        gAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGt                                           <  1:2792630/71‑1 (MQ=255)
                         aaTTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGt                                           >  1:1388129/1‑70 (MQ=255)
                           ttattCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGAGCTGGGAAATGAag                                                      >  1:1940886/1‑57 (MQ=255)
                              ttCCCGCTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGt                                           <  1:2603657/65‑1 (MQ=255)
                                cccGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGAGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAg                                   >  1:1426356/1‑71 (MQ=255)
                                 ccGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGAGCt                                                                 >  1:3880088/1‑40 (MQ=255)
                                   ggTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGAGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGcc                                 >  1:2480941/1‑70 (MQ=255)
                                     ttCGATGGGAGCAAAAAGCTTTATCGTCCGTGAGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCa                                <  1:1269009/69‑1 (MQ=255)
                                          tGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGt                                           <  1:444547/53‑1 (MQ=255)
                                          tGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGt                                           <  1:823017/53‑1 (MQ=255)
                                          tGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGcc                         <  1:1083902/71‑1 (MQ=255)
                                           gggAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTc                                              >  1:3474021/1‑49 (MQ=255)
                                           gggAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGt                            <  1:1812938/67‑1 (MQ=255)
                                            ggAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTc                                         >  1:2393166/1‑53 (MQ=255)
                                            ggAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCa                                    >  1:2992674/1‑58 (MQ=255)
                                            ggAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcg                        <  1:2500493/70‑1 (MQ=255)
                                            ggAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcg                        <  1:2982745/70‑1 (MQ=255)
                                              aGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgggc                     <  1:802019/71‑1 (MQ=255)
                                               gCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGc                                     <  1:265659/54‑1 (MQ=255)
                                               gCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGc                                     <  1:1382742/54‑1 (MQ=255)
                                               gCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgggcg                    <  1:2935270/71‑1 (MQ=255)
                                                   aaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGt                                       >  1:1913938/1‑48 (MQ=255)
                                                    aaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTaa               <  1:2249705/71‑1 (MQ=255)
                                                     aGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAAt              >  1:2731358/1‑71 (MQ=255)
                                                      gCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATg             <  1:2819381/71‑1 (MQ=255)
                                                        tttATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTa                <  1:2857186/66‑1 (MQ=255)
                                                         ttATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGc          >  1:2402911/1‑71 (MQ=255)
                                                          tATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAg           >  1:3238509/1‑69 (MQ=255)
                                                          tATCGTCCGTGAGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGcc         >  1:2921123/1‑71 (MQ=255)
                                                            tCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCGa       <  1:1596499/71‑1 (MQ=255)
                                                             cGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCg        <  1:1195964/69‑1 (MQ=255)
                                                              gTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCGAAc     >  1:354753/1‑71 (MQ=255)
                                                               tCCGTTGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCGAACt    <  1:440708/71‑1 (MQ=255)
                                                                 cGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCGAACTaa  <  1:2702686/71‑1 (MQ=255)
                                                                 cGTGAGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCGAACTa   <  1:1203654/70‑1 (MQ=255)
                                                                  gTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCGAACTaa  >  1:3208119/1‑70 (MQ=255)
                                                                    gggCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCGAACTa   <  1:3174323/67‑1 (MQ=255)
                                                                     |                                                                  
TGTCTGCGCGTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCGAACTAA  >  minE/2663921‑2664056

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: