Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,603,289 T→C 67.4% V105V (GTA→GTG yfeX ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,603,2890TC67.4% 47.0 / 27.5 58V105V (GTA→GTGyfeXconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (5/14);  new base C (28/11);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 1.67e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

ACCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGAACTGGGTCGTCGGCGCAAGGCCTT  >  minE/1603250‑1603353
                                       |                                                                
aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGagaaagaa                                            <  1:3390710/62‑1 (MQ=255)
aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGagaaagaa                                            <  1:2532979/62‑1 (MQ=255)
aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGagaaagaa                                            <  1:590872/62‑1 (MQ=255)
aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTgg                                       >  1:223445/1‑67 (MQ=255)
aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTgg                                       >  1:1775274/1‑67 (MQ=255)
aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGa                                      >  1:63077/1‑68 (MQ=255)
aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGa                                      >  1:3151881/1‑68 (MQ=255)
aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACATCGTGACGCAGAGAAAGAATGTGGa                                      >  1:600134/1‑68 (MQ=255)
  cAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATc                                    <  1:15764/68‑1 (MQ=255)
    aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGACTGTGGa                                      <  1:2816/64‑1 (MQ=255)
    aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa                                      <  1:932031/64‑1 (MQ=255)
    aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa                                      <  1:864872/64‑1 (MQ=255)
    aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa                                      <  1:2999388/64‑1 (MQ=255)
    aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa                                      <  1:3059031/64‑1 (MQ=255)
    aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa                                      <  1:3331374/64‑1 (MQ=255)
    aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa                                      <  1:484490/64‑1 (MQ=255)
    aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa                                      <  1:155709/64‑1 (MQ=255)
    aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa                                      <  1:1271022/64‑1 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:1035698/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:3291628/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:3287824/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:3269829/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:417791/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:304938/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:878064/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:2980380/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:2722180/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:618212/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:2268555/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:177270/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:1447878/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:2158382/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:1444097/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:2212989/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:1142863/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:2224051/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:2278924/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:2313071/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:111971/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:2329453/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:1114636/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:1007172/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:249508/1‑57 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaaga                                             >  1:562614/1‑56 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa                                      <  1:2618757/63‑1 (MQ=255)
     aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGAAGCagag                                                  >  1:1714011/1‑51 (MQ=255)
     aGGCTTACATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa                                            >  1:425651/1‑57 (MQ=255)
             aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa                          <  1:1788709/67‑1 (MQ=255)
             aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa                          <  1:2039998/67‑1 (MQ=255)
             aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa                          <  1:2098250/67‑1 (MQ=255)
             aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa                          <  1:215028/67‑1 (MQ=255)
             aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa                          <  1:3223093/67‑1 (MQ=255)
             aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa                          <  1:1236697/67‑1 (MQ=255)
             aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa                          <  1:3068108/67‑1 (MQ=255)
             aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa                          <  1:2675735/67‑1 (MQ=255)
                 ccgccTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa                                      <  1:2904597/51‑1 (MQ=255)
                                agaAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGAACTGGGTCGTCGGCGCAAGGc     <  1:674330/69‑1 (MQ=255)
                                   aGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGAACTGGGTCGTCGGCGCAAGGCCtt  <  1:2427565/69‑1 (MQ=255)
                                       |                                                                
ACCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGAACTGGGTCGTCGGCGCAAGGCCTT  >  minE/1603250‑1603353

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: