Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 752,626 G→A 100% D198D (GAC→GAT mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE752,6260GA75.1% 39.2 / 9.1 36D198D (GAC→GATmfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base G (2/7);  new base A (0/27);  total (2/34)
Fisher's exact test for biased strand distribution p-value = 5.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTGCGCTGGCTGTCGACGTCAAACACCCGCAGGCTGTCGATTTCATCATCAAAGAAATCAAGACGATAAGGCAGCTCAC  >  minE/752567‑752693
                                                           |                                                                   
ggAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTGCGCTGGCTGTCGACGTCa                                                            <  1:806823/69‑1 (MQ=255)
ggAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTGCGCTGGCTGTCGACGTCa                                                            <  1:2719855/69‑1 (MQ=255)
         tGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTGCGCTGGCTGTCGACGTCAAACACCCg                                                    <  1:273683/68‑1 (MQ=255)
         tGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTGCGCTGGCTGTCGACGTCAAACACCCg                                                    <  1:563030/68‑1 (MQ=255)
         tGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTGCGCTGGCTGTCGACGTCAAACACCCg                                                    <  1:2429814/68‑1 (MQ=255)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:3286968/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:299757/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:3051088/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:3124493/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:312488/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:3201267/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:3238831/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:2972514/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:390939/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:700202/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:7070/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:864364/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:878918/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:977538/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:1927790/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:1101155/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:1558523/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:1564572/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:1576455/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:1823591/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:1070682/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:1991872/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:2009012/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:2151956/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:2528174/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:2938349/44‑1 (MQ=37)
                      gATTGATCTCTTCTACTTCCTCCAGCGTGCGCTGGCTATCGACg                                                               <  1:2939489/44‑1 (MQ=37)
                         tGATCGCTTCTACTTCCTCCAGCGTGCGCTGGCTGTCGACGTCAAACACCCGCAGGCTGTCGATTtcat                                   >  1:804686/1‑69 (MQ=255)
                         tGATCGCTTCTACTTCCTCCAGCGTGCGCTGGCTGTCGACGTCAAACACCCGCAGGCTGTCGATTtcat                                   >  1:1494052/1‑69 (MQ=255)
                                    aCTTCCTCCAGCGTGCGCTGGCTGTCGACGTCAAACACCCGCAGGCTGTCGATTTCATCATCAAAGa                          <  1:165941/67‑1 (MQ=255)
                                              gCGTGCGCTGGCTGTCGACGTCAAACACCCGCAGGCTGTCGATTTCATCATCAAAGAAATCAAGACGAt              <  1:3222298/69‑1 (MQ=255)
                                                          tgTCGACGTCAAACACCCGCAGGCTGTCGATTTCATCATCAAAGAAATCAAGACGATAAGGCAGCTCAc  <  1:2217949/69‑1 (MQ=255)
                                                          tgTCGACGTCAAACACCCGCAGGCTGTCGATTTCATCATCAAAGAAATCAAGACGATAAGGCAGCTCAc  <  1:831954/69‑1 (MQ=255)
                                                           |                                                                   
GGAAATTCGTGCGCGGGCAGCAGATTGATCGCTTCTACTTCCTCCAGCGTGCGCTGGCTGTCGACGTCAAACACCCGCAGGCTGTCGATTTCATCATCAAAGAAATCAAGACGATAAGGCAGCTCAC  >  minE/752567‑752693

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