Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,340,208 A→G 40.0% intergenic (‑2/‑371) bipA ← / → glnA GTP‑binding protein/glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,340,2080AG40.0% 17.0 / 31.2 30intergenic (‑2/‑371)bipA/glnAGTP‑binding protein/glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base A (5/13);  new base G (5/7);  total (10/20)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

ATGGTCTACGTGCGCGATGATGGCGATATTACGCAATTTTTCGATCACAACTTTGCCTCAGGCATTAGAAATAGCGCGTTATTGTACACGGATTAATCGCACTACAA  >  minE/2340159‑2340265
                                                 |                                                         
aTGGTCTACGTGCGCGATGATGGCGATATTACGCAATTTTTCGATCACAACTTTGCCTCAg                                                <  1:583442/61‑1 (MQ=255)
aTGGTCTACGTGCGCGATGATGGCGATATTACGCAATTTTTCGATCACAACTTTGCCTCAg                                                <  1:581521/61‑1 (MQ=255)
aTGGTCTACGTGCGCGATGATGGCGATATTACGCAATTTTTCGATCACAACTTTGCCTCAg                                                <  1:75264/61‑1 (MQ=255)
aTGGTCTACGTGCGCGATGATGGCGATATTACGCAATTTTTCGATCACAACTTTGCCTCAg                                                <  1:274348/61‑1 (MQ=255)
   gTCTACGTGCGCGATGATGGCGATATTACGCAATTTTTCGATCACAACTTTGCCTCAgg                                               <  1:437378/59‑1 (MQ=255)
          tGCGCGATGATGGCGATATTACGCAATTTTTCGATCACAACTTTGCCTCAGGCATTAGAAAt                                     <  1:19580/62‑1 (MQ=255)
          tGCGCGATGATGGCGATATTACGCAATTTTTCGATCACAACTTTGCCTCAGGCATTAGAAAt                                     <  1:263907/62‑1 (MQ=255)
              cgATGATGGCGATATTACGCAATTTTTCGATCACAACTTTGCCTCAGGCATTAGAAATAgc                                  <  1:1749/61‑1 (MQ=255)
              cgATGATGGCGATATTACGCAATTTTTCGATCACAACTTTGCCTCAGGCATTAGAAATAgc                                  <  1:233898/61‑1 (MQ=255)
              cgATGATGGCGATATTACGCAATTTTTCGATCACAACTTTGCCTCAGGCATTAGAAATAgc                                  <  1:294785/61‑1 (MQ=255)
              cgATGATGGCGATATTACGCAATTTTTCGATCACAACTTTGCCTCAGGCATTAGAAATAgc                                  <  1:304669/61‑1 (MQ=255)
               gatgatGGCGATATTACGCAATTTTTCGATCACAACTTTGCCTCAGGCATTAGAAATAgcgc                                <  1:307598/62‑1 (MQ=255)
               gatgatGGCGATATTACGCAATTTTTCGATCACAACTTTGCCTCAGGCATTAGAAATAgcgc                                <  1:74817/62‑1 (MQ=255)
                  gatgGCGATATTACGCAATTTTTCGATCACAGCTTTGCCTCAgg                                               <  1:339853/44‑1 (MQ=255)
                  gatgGCGATATTACGCAATTTTTCGATCACAGCTTTGCCTCAgg                                               <  1:113957/44‑1 (MQ=255)
                  gatgGCGATATTACGCAATTTTTCGATCACAGCTTTGCCTCAgg                                               <  1:419398/44‑1 (MQ=255)
                  gatgGCGATATTACGCAATTTTTCGATCACAGCTTTGCCTCAgg                                               <  1:542370/44‑1 (MQ=255)
                  gatgGCGATATTACGCAATTTTTCGATCACAGCTTTGCCTCAgg                                               <  1:577517/44‑1 (MQ=255)
                  gatgGCGATATTACGCAATTTTTCGATCACAGCTTTGCCTCAgg                                               <  1:182678/44‑1 (MQ=255)
                  gatgGCGATATTACGCAATTTTTCGATCACAGCTTTGCCTCAgg                                               <  1:158544/44‑1 (MQ=255)
                         atatTACGCAATTTTTCGATCACAGCTTTGCCTCAGGCATTAGa                                        >  1:214680/1‑44 (MQ=255)
                         atatTACGCAATTTTTCGATCACAGCTTTGCCTCAGGCATTAGa                                        >  1:409409/1‑44 (MQ=255)
                         atatTACGCAATTTTTCGATCACAGCTTTGCCTCAGGCATTAGa                                        >  1:194469/1‑44 (MQ=255)
                         atatTACGCAATTTTTCGATCACAGCTTTGCCTCAGGCATTAGa                                        >  1:183746/1‑44 (MQ=255)
                         atatTACGCAATTTTTCGATCACAGCTTTGCCTCAGGCATTAGa                                        >  1:518916/1‑44 (MQ=255)
                                             cacaACTTTGCCTCAGGCATTAGAAATAGCGCGTTATTGTACACGGATTAATCGCACTACaa  >  1:76602/1‑62 (MQ=255)
                                             cacaACTTTGCCTCAGGCATTAGAAATAGCGCGTTATTGTACACGGATTAATCGCACTACaa  >  1:585898/1‑62 (MQ=255)
                                             cacaACTTTGCCTCAGGCATTAGAAATAGCGCGTTATTGTACACGGATTAATCGCACTACaa  >  1:478403/1‑62 (MQ=255)
                                             cacaACTTTGCCTCAGGCATTAGAAATAGCGCGTTATTGTACACGGATTAATCGCACTACaa  >  1:354240/1‑62 (MQ=255)
                                             cacaACTTTGCCTCAGGCATTAGAAATAGCGCGTTATTGTACACGGATTAATCGCACTACa   >  1:350629/1‑61 (MQ=255)
                                                 |                                                         
ATGGTCTACGTGCGCGATGATGGCGATATTACGCAATTTTTCGATCACAACTTTGCCTCAGGCATTAGAAATAGCGCGTTATTGTACACGGATTAATCGCACTACAA  >  minE/2340159‑2340265

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: