Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,252,050 T→C 26.9% K6K (AAA→AAG yebC ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,252,0500TC26.9% 36.7 / 15.9 26K6K (AAA→AAGyebCconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (13/6);  new base C (4/3);  total (17/9)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

ATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAA  >  minE/1251985‑1252115
                                                                 |                                                                 
aTTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTAc                                                               >  1:244640/1‑70 (MQ=255)
   ttAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTAt                                                            <  1:784158/70‑1 (MQ=255)
    tAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATg                                                           <  1:2144073/70‑1 (MQ=255)
              tttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATa                                                 >  1:698812/1‑70 (MQ=255)
              tttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGAcc                                                        >  1:889206/1‑63 (MQ=255)
              tttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGAcc                                                        >  1:758892/1‑63 (MQ=255)
               ttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTAt                                                            >  1:1964868/1‑58 (MQ=255)
               ttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTTcc                                                    >  1:556770/1‑66 (MQ=255)
                  ccGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATaaaaa                                             <  1:311125/70‑1 (MQ=255)
                     cgcTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCAtt                                                                  <  1:871938/46‑1 (MQ=255)
                     cgcTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCAtt                                                                  <  1:1345735/46‑1 (MQ=255)
                     cgcTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCAtt                                                                  <  1:1027202/46‑1 (MQ=255)
                          aGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCa                                     <  1:1798935/70‑1 (MQ=255)
                            cATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGa                                                          >  1:97327/1‑47 (MQ=255)
                            cATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGa                                                          >  1:2073121/1‑47 (MQ=255)
                              tCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAg                                <  1:1984996/71‑1 (MQ=255)
                                 tGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTa                                                             >  1:1711222/1‑39 (MQ=255)
                                 tGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAatat                                          >  1:411378/1‑58 (MQ=255)
                                  gcgcAGCTTTACGATGTCTGGTGTTGGCCCATTTAc                                                               >  1:1146880/1‑36 (MQ=255)
                                   cgcAGCTTTACGATGTCTGGTGTTGGCCCATTTAc                                                               <  1:288584/35‑1 (MQ=255)
                                   cgcAGCTTTACGATGTCTGGTGTTGGCCCATTTAc                                                               <  1:1545861/35‑1 (MQ=255)
                                          ttACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCcctgcc                         >  1:1108965/1‑66 (MQ=255)
                                          ttACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCcctgcc                         >  1:930535/1‑66 (MQ=255)
                                          ttACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGtt                     >  1:1126677/1‑70 (MQ=255)
                                          ttACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGtt                     >  1:1064634/1‑70 (MQ=255)
                                                  tCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATAtctc                                       <  1:1731237/44‑1 (MQ=255)
                                                  tCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAg                                <  1:409592/51‑1 (MQ=255)
                                                        gtTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAAt       >  1:1162434/1‑70 (MQ=255)
                                                              ccATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACaa  >  1:769913/1‑69 (MQ=255)
                                                                 tttACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTc                    <  1:2185148/48‑1 (MQ=255)
                                                                 tttACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGc                <  1:978950/52‑1 (MQ=255)
                                                                 |                                                                 
ATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAA  >  minE/1251985‑1252115

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: