Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,834,068 G→A 22.2% L56L (CTG→CTA ygbM → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,834,0680GA22.2% 66.9 / 15.8 36L56L (CTG→CTAygbMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/6);  new base A (5/3);  total (27/9)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGTTTCCCTATAACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGAC  >  minE/1834002‑1834135
                                                                  |                                                                   
tCTGTTTCCCTATAACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTAgcg                                                                  <  1:362209/70‑1 (MQ=255)
    tttCCCTATAACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGtt                                                             <  1:2493620/71‑1 (MQ=255)
       cccTATAACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCaa                                                          >  1:3361306/1‑71 (MQ=255)
        ccTATAACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTAGCGCTGTTCAAt                                                         <  1:76116/71‑1 (MQ=255)
          tataACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTg                                                                     >  1:2962447/1‑57 (MQ=255)
          tataACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATAc                                                       >  1:3040470/1‑71 (MQ=255)
          tataACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTAGCGCTGTTCaa                                                          <  1:3357274/68‑1 (MQ=255)
            taACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTAGCGCTGTTCAATa                                                        >  1:135315/1‑68 (MQ=255)
                tactCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTcc                                                 >  1:3701500/1‑71 (MQ=255)
                tactCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTcc                                                 >  1:3544837/1‑71 (MQ=255)
                   tCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTc                                                            <  1:2264087/57‑1 (MQ=255)
                     cACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCaa                                                          >  1:101734/1‑57 (MQ=255)
                     cACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCg                                               >  1:3272279/1‑68 (MQ=255)
                         cTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGatat                                        >  1:1187577/1‑71 (MQ=255)
                            cAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGgc                                                                   >  1:877730/1‑41 (MQ=255)
                             aaaTCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAAt                                    >  1:2414328/1‑71 (MQ=255)
                             aaaTCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAAt                                    >  1:3314554/1‑71 (MQ=255)
                               aTCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGc                                  >  1:1974493/1‑71 (MQ=255)
                                 ccAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTa                                      >  1:968820/1‑65 (MQ=255)
                                 ccAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTa                                      >  1:1993020/1‑65 (MQ=255)
                                     aaGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTaa                                     <  1:466075/62‑1 (MQ=255)
                                     aaGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTaa                                     <  1:3028301/62‑1 (MQ=255)
                                     aaGCAACTGGAGCAAAACCATCTGACACTAGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCgggg                              >  1:3360936/1‑69 (MQ=255)
                                      aGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGaa                           >  1:3848287/1‑71 (MQ=255)
                                              gAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGAt                    >  1:504514/1‑70 (MQ=255)
                                                gCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATAtt                                       <  1:423225/49‑1 (MQ=255)
                                                 cAAAACCATCTGACACTAGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATc                >  1:1122930/1‑71 (MQ=255)
                                                  aaaaCCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATcc               <  1:3130588/71‑1 (MQ=255)
                                                   aaaCCATCTGACACTAGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCg              >  1:288238/1‑71 (MQ=255)
                                                   aaaCCATCTGACACTAGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCg              >  1:2863498/1‑71 (MQ=255)
                                                     aCCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCg                                 >  1:2885087/1‑50 (MQ=255)
                                                        aTCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCtt         >  1:2183903/1‑71 (MQ=255)
                                                          cTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTcc       >  1:3901482/1‑71 (MQ=255)
                                                          cTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTcc       >  1:989089/1‑71 (MQ=255)
                                                              cacTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGa   >  1:2211168/1‑71 (MQ=255)
                                                               acTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGAc  >  1:3270243/1‑71 (MQ=255)
                                                                cTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGc             >  1:299797/1‑59 (MQ=255)
                                                                  |                                                                   
TCTGTTTCCCTATAACTACTCCACCCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGAC  >  minE/1834002‑1834135

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: