Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,097,909 T→C 37.5% intergenic (+32/‑397) ydiK → / → ydiL predicted inner membrane protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,097,9090TC37.5% 15.9 / 22.7 24intergenic (+32/‑397)ydiK/ydiLpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (9/6);  new base C (4/5);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGT  >  minE/1097839‑1097979
                                                                      |                                                                      
ctcGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGt                                                                        >  1:1992355/1‑71 (MQ=255)
            ggCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGCTAACCACCGCt                                                             >  1:646773/1‑70 (MQ=255)
                  aTCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGat                                                      <  1:845856/71‑1 (MQ=255)
                   tCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAac                                                                    <  1:1924486/56‑1 (MQ=255)
                            cGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGata                                                     <  1:1915246/62‑1 (MQ=255)
                            cGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGCTAACCACCGCt                                                             <  1:1006657/54‑1 (MQ=255)
                              aaTAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGCTAACCACCGATaa                                                           >  1:947467/1‑54 (MQ=255)
                               aTAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGctct                                         >  1:1222076/1‑71 (MQ=255)
                                tAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCtt                                        <  1:1444061/71‑1 (MQ=255)
                                        ttCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCatat                                >  1:2156794/1‑71 (MQ=255)
                                             atcaGGCGGCTCTGCCGCCTGATTGCTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAttt                              <  1:397356/68‑1 (MQ=255)
                                                  gcggcTCTGCCGCCTGATTGCTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCtt                      <  1:1371706/71‑1 (MQ=255)
                                                   cggcTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTa                     >  1:1876153/1‑71 (MQ=255)
                                                       tctGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGc                 >  1:350302/1‑71 (MQ=255)
                                                         tGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCa               >  1:1044463/1‑71 (MQ=255)
                                                            cgccTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTcc                           <  1:1099315/56‑1 (MQ=255)
                                                              ccTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAg          >  1:253851/1‑71 (MQ=255)
                                                              ccTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAg          >  1:2055154/1‑71 (MQ=255)
                                                                tGATTGCTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAttt                              <  1:844493/49‑1 (MQ=255)
                                                                tGATTGCTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATAttt                              <  1:880120/49‑1 (MQ=255)
                                                                 gATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAg          <  1:348267/68‑1 (MQ=255)
                                                                  aTTGCTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTgg      >  1:687346/1‑71 (MQ=255)
                                                                  aTTGCTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTgg      >  1:1648764/1‑71 (MQ=255)
                                                                     gTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTgg      >  1:1849663/1‑68 (MQ=255)
                                                                     gTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtag   >  1:2005487/1‑71 (MQ=255)
                                                                      ttAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGgtagt  <  1:2186606/71‑1 (MQ=255)
                                                                      |                                                                      
CTCGAAGAACTTGGCGAAATCGAGAAACCGAATAAGTAATTTCTCATCAGGCGGCTCTGCCGCCTGATTGTTAACCACCGCTAATTGATATTTGCTTGCTCTTCCCCATATTTTCCTGCTTACCGCCATTCAGCTGGTAGT  >  minE/1097839‑1097979

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: