Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,582,181 G→A 30.5% L352L (CTG→TTG)  ligA ← DNA ligase, NAD(+)‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,582,1810GA30.5% 74.4 / 34.6 79L352L (CTG→TTG) ligADNA ligase, NAD(+)‑dependent
Reads supporting (aligned to +/- strand):  ref base G (13/41);  new base A (22/2);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 2.74e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTA  >  minE/1582113‑1582236
                                                                    |                                                       
tttATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:2438728/69‑1 (MQ=255)
tttATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:2085606/69‑1 (MQ=255)
tttATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:2073525/69‑1 (MQ=255)
tttATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:618139/69‑1 (MQ=255)
tttATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:792652/69‑1 (MQ=255)
tttATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:1791184/69‑1 (MQ=255)
      gCCAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGcc                                                   <  1:2775234/69‑1 (MQ=255)
      gCCAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGcc                                                   <  1:2828017/69‑1 (MQ=255)
      gCCAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGcc                                                   <  1:726617/69‑1 (MQ=255)
      gCCAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGcc                                                   <  1:153746/69‑1 (MQ=255)
      gCCAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGcc                                                   <  1:628998/69‑1 (MQ=255)
       ccAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGcc                                                   <  1:1035302/68‑1 (MQ=255)
        cAATGCGTAAACCAAGACGTTCGATTTCATCGGCATTGTGTAAGGTTGCGTTACTCACCAGCACGcg                                                   <  1:603487/67‑2 (MQ=255)
             cGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:3157884/56‑1 (MQ=255)
             cGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:352885/56‑1 (MQ=255)
             cGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:335789/56‑1 (MQ=255)
             cGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:1737280/56‑1 (MQ=255)
             cGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:196849/56‑1 (MQ=255)
             cGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:2398773/56‑1 (MQ=255)
             cGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:398976/56‑1 (MQ=255)
             cGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:1131517/56‑1 (MQ=255)
             cGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:1104944/56‑1 (MQ=255)
             cGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACa                                            >  1:1417297/1‑69 (MQ=255)
             cGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACa                                            >  1:2167466/1‑69 (MQ=255)
             cGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGAAACa                                            >  1:1001810/1‑69 (MQ=255)
              gTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGcc                                                   <  1:971103/61‑1 (MQ=255)
               tAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAg                                                         <  1:2500205/54‑1 (MQ=255)
                                     tCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAAc                     <  1:1013036/68‑1 (MQ=255)
                                     tCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAAc                     <  1:2400471/68‑1 (MQ=255)
                                     tCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAAc                     <  1:2722131/68‑1 (MQ=255)
                                     tCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAAc                     <  1:3082744/68‑1 (MQ=255)
                                     tCCGCATTGTGTAAGGTTGCGATACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAAc                     <  1:637642/68‑1 (MQ=255)
                                       cGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGtt                                  <  1:1128290/53‑1 (MQ=255)
                                       cGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGtt                                  <  1:900523/53‑1 (MQ=255)
                                       cGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGtt                                  <  1:129965/53‑1 (MQ=255)
                                       cGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGtt                                  <  1:3379230/53‑1 (MQ=255)
                                       cGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGtt                                  <  1:1821364/53‑1 (MQ=255)
                                       cGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGtt                                  <  1:2305688/53‑1 (MQ=255)
                                       cGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGtt                                  <  1:2359583/53‑1 (MQ=255)
                                       cGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGtt                                  <  1:1640569/53‑1 (MQ=255)
                                       cGCATTGTGTAAGGTTGCGTTACACACCAGCACGCCTGCAACATGGACAGGtt                                  <  1:994402/53‑1 (MQ=255)
                                        gCATTGTGTAAGGTTGCGTTACTCACCATCACGCCTGCAACATGGACAGGtt                                  <  1:1549934/52‑1 (MQ=255)
                                        gCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGtt                                  <  1:379608/52‑1 (MQ=255)
                                           ttGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTa              <  1:146472/69‑1 (MQ=255)
                                           ttGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTa              <  1:440979/69‑1 (MQ=255)
                                           ttGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTa              <  1:678785/69‑1 (MQ=255)
                                           ttGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTa              <  1:3181124/69‑1 (MQ=255)
                                           ttGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTa              <  1:2548250/69‑1 (MQ=255)
                                           ttGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTa              <  1:3021579/69‑1 (MQ=255)
                                           ttGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTa              <  1:2619999/69‑1 (MQ=255)
                                                    gTTGCGTTACTCACCAGCACGCCTGTAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGcc      >  1:3012179/1‑68 (MQ=255)
                                                    gTTGCGTTACTCACCAGCACGCCTGCAACATGGACAgg                                    >  1:2573754/1‑38 (MQ=255)
                                                    gTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGcc      >  1:2260235/1‑68 (MQ=255)
                                                    gTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGcc      >  1:2730428/1‑68 (MQ=255)
                                                    gTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGcc      >  1:974907/1‑68 (MQ=255)
                                                    gTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGcc      >  1:1441110/1‑68 (MQ=255)
                                                    gTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGcc      >  1:2831898/1‑68 (MQ=255)
                                                    gTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGcc      >  1:1051513/1‑68 (MQ=255)
                                                    gTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGcc      >  1:2523739/1‑68 (MQ=255)
                                                    gTTGAGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGcc      >  1:1877196/1‑68 (MQ=255)
                                                       gCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTa  <  1:1019647/69‑1 (MQ=255)
                                                       gCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTa  <  1:933728/69‑1 (MQ=255)
                                                       gCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTa  <  1:438976/69‑1 (MQ=255)
                                                       gCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTa  <  1:2492369/69‑1 (MQ=255)
                                                       gCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTa  <  1:3414893/69‑1 (MQ=255)
                                                       gCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTa  <  1:2385946/69‑1 (MQ=255)
                                                       gCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTa  <  1:3250951/69‑1 (MQ=255)
                                                       gCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTa  <  1:13931/69‑1 (MQ=255)
                                                       gCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTa  <  1:2984594/69‑1 (MQ=255)
                                                       gCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTa  <  1:2872905/69‑1 (MQ=255)
                                                       gCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTa  <  1:2138876/69‑1 (MQ=255)
                                                           tACTCACCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTaa             <  1:2302127/54‑1 (MQ=255)
                                                           tACTCACCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTaa             <  1:2769330/54‑1 (MQ=255)
                                                                aCCAACACGCCAGCATCATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:2202732/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCaa                      >  1:833751/1‑40 (MQ=38)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAgga                 >  1:402630/1‑44 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:78043/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:1884844/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:1408241/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:681775/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:1052437/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:1543638/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:552911/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:495960/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:1644792/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:390327/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:1768031/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:3231059/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:1897208/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:2026947/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:3008506/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:2690425/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:2485480/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:2363223/1‑45 (MQ=255)
                                                                aCCAACACGCCAGCAACATGGACAGGTTCCAGACGCGCAACAGGc                 >  1:2424496/1‑45 (MQ=255)
                                                                    |                                                       
TTTATCGCCAATGCGTAAACCAAGACGTTCGATTTCATCCGCATTGTGTAAGGTTGCGTTACTCACCAGCACGCCTGCAACATGGACAGGTTCCAGACGCGCAACAGGCGTAATCGCGCCAGTA  >  minE/1582113‑1582236

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: