Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,363,056 G→C 47.7% T182T (ACG→ACC metG → methionyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,363,0560GC47.7% 5.8 / 107.2 88T182T (ACG→ACCmetGmethionyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (2/44);  new base C (16/26);  total (18/70)
Fisher's exact test for biased strand distribution p-value = 1.08e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGTGATTCTGAACACTTCTTCTTTGATCTGCCCTCTTTCAGCGAAAT  >  minE/1362993‑1363112
                                                               |                                                        
cTGCGGCGCGACCTACAGGCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:1239841/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGTTGGTTTCTGGCGCTACGCCGGt                                                     <  1:2105692/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:300636/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:1845061/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:1923409/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:2200670/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:2260426/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:2420250/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:1000471/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:2574603/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:2744965/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:2475219/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:3013620/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:3279622/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:3337403/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:365730/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:426631/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:477754/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:995402/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:1255912/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:1632170/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:1097194/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:1700008/69‑1 (MQ=255)
cTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:1796054/69‑1 (MQ=255)
 tGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:1097352/68‑1 (MQ=255)
   cggcgCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:770810/66‑1 (MQ=255)
    ggcgCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     <  1:3380526/65‑1 (MQ=255)
       gcgACCTACATCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGt                                              <  1:2595265/69‑1 (MQ=255)
       gcgACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGt                                              <  1:2599212/69‑1 (MQ=255)
       gcgACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGt                                              <  1:1917503/69‑1 (MQ=255)
       gcgACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGt                                              <  1:2103363/69‑1 (MQ=255)
       gcgACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGt                                              <  1:1423475/69‑1 (MQ=255)
       gcgACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGt                                              <  1:2951292/69‑1 (MQ=255)
       gcgACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGt                                              <  1:359786/69‑1 (MQ=255)
               cAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGt                                                     >  1:2081037/1‑54 (MQ=255)
                          aacttaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGATTCTGAACActtcttct                           >  1:452043/1‑69 (MQ=255)
                          aactgaTTGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGATTCTGAACActtcttct                           >  1:1381956/1‑69 (MQ=255)
                          aactgaTCGGGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:2748140/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:1133250/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:2972525/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:1334228/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:3007219/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:1321128/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:3158675/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:3243471/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:1827512/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:425711/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:586328/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:783064/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:837966/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:923956/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:1993186/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:2180444/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:2162823/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:2740907/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:2349182/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:2154145/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:245512/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:186479/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:2731908/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:2689065/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:1505004/54‑1 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGATTCTGAACActtcttct                           >  1:3119221/1‑69 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGATTCTGAACActtcttct                           >  1:2904074/1‑69 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGATTCTGAACActtcttct                           >  1:1980258/1‑69 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGATTCTGAACActtcttct                           >  1:776221/1‑69 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGATTCTGAACActtcttct                           >  1:2246760/1‑69 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGATTCTGAACActtcttct                           >  1:2531096/1‑69 (MQ=255)
                          aactgaTCGAGCCGAAATCGGTGGTATCTGGTGCTACCCCGGTAATGCGTGAtt                                          <  1:3412289/54‑1 (MQ=255)
                          aactgaTCGAGCAGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGATTCTGAACActtcttct                           >  1:2777959/1‑69 (MQ=255)
                               aTCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGTGa                                            <  1:3013803/47‑1 (MQ=255)
                               aTCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGTGa                                            <  1:3268336/47‑1 (MQ=255)
                               aTCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGTGa                                            <  1:1256834/47‑1 (MQ=255)
                               aTCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGTGa                                            <  1:1588169/47‑1 (MQ=255)
                               aTCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGTGATTCTGAACACTTCTTCTTTGAt                      >  1:2267110/1‑69 (MQ=255)
                                tCGAGCCGAAATCGGTGGTTTCTGGCTCTACGCCGGTAATGCGt                                              <  1:445892/44‑1 (MQ=255)
                                tCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGTGATTCTGAaca                                   <  1:2395739/55‑1 (MQ=255)
                                   aGCCGAAATCGGTGGTTTCTTGTGCTACCCCGGTAATGCGTGATTCTGAACACTTCTTCTTTGATCTGc                  >  1:2431530/1‑69 (MQ=255)
                                   aGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGATTCTGAACACTTCTTCTTTGATCTGc                  >  1:3079380/1‑69 (MQ=255)
                                   aGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGATTCTGAACACTTCTTCTTTGATCTGc                  >  1:1556549/1‑69 (MQ=255)
                                   aGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGATTCTGAACACTTCTTCTTTGATCTGc                  >  1:617242/1‑69 (MQ=255)
                                   aGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGATTCTGAACACTTCTTCTTTGATCTGc                  >  1:1887362/1‑69 (MQ=255)
                                   aGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGATTCTGAACACTTCTTCTTTGATCTGc                  >  1:2872782/1‑69 (MQ=255)
                                   aGCCGAAATCGGTGGTTTCTGGTGCTACCCCGGTAATGCGTGATTCTGAACACTTCTTCTTTGATCTGc                  >  1:1824474/1‑69 (MQ=255)
                                    gCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGTGATTCTGAACACTTCTTCTTTGATCTGc                  <  1:1565139/68‑1 (MQ=255)
                                                 gtTTCTGGCGCTACGCCGGTAATGCGTGATTCTGAACACTTCTTCTTTGATCTGCCCTCTTTCAGCGaa    <  1:1882608/69‑1 (MQ=255)
                                                   ttCTGGCGCTACGCCGGTAATGCGTGATTCTGAACACTTCTTCTTTGATCTGCCCTCTTTCAGCGAAAt  <  1:1069698/69‑1 (MQ=255)
                                                        gcgcTACGCCGGTAATGCGTGATTCTGAACACTTCTTCTTTGATCTGCCCTCTTTCAGCGAAAt  <  1:465746/64‑1 (MQ=255)
                                                               |                                                        
CTGCGGCGCGACCTACAGCCCGACTGAACTGATCGAGCCGAAATCGGTGGTTTCTGGCGCTACGCCGGTAATGCGTGATTCTGAACACTTCTTCTTTGATCTGCCCTCTTTCAGCGAAAT  >  minE/1362993‑1363112

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: