Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 355,392 G→A 36.5% R79R (CGC→CGT ybcI ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE355,3920GA36.5% 73.9 / 68.0 112R79R (CGC→CGTybcIconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (39/32);  new base A (10/31);  total (49/63)
Fisher's exact test for biased strand distribution p-value = 2.79e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 115 total aligned reads displayed.

GCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGAT  >  minE/355331‑355459
                                                             |                                                                   
gCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCa                                                               >  1:2154970/1‑68 (MQ=255)
gCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCa                                                               >  1:825328/1‑68 (MQ=255)
gCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCa                                                               >  1:2375138/1‑68 (MQ=255)
gCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCa                                                               >  1:2717887/1‑68 (MQ=255)
gCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCa                                                               >  1:1574160/1‑68 (MQ=255)
gCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCa                                                               >  1:1553083/1‑68 (MQ=255)
gCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCa                                                               >  1:1272379/1‑68 (MQ=255)
gCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCa                                                               >  1:2871851/1‑68 (MQ=255)
     aGCAACTAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCTGCGTCCAATa                                                            >  1:2110406/1‑66 (MQ=255)
     aGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAAt                                                             >  1:2967291/1‑65 (MQ=255)
     aGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATa                                                            >  1:1052548/1‑66 (MQ=255)
     aGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATa                                                            >  1:3279862/1‑66 (MQ=255)
     aGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATa                                                            >  1:3378580/1‑66 (MQ=255)
     aGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATa                                                            >  1:354367/1‑66 (MQ=255)
     aGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATa                                                            >  1:56339/1‑66 (MQ=255)
     aGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATa                                                            >  1:595032/1‑66 (MQ=255)
     aGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATa                                                            >  1:1151225/1‑66 (MQ=255)
     aGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCAAGCCCTGAACCATCGGCGTCCAATa                                                            >  1:367962/1‑66 (MQ=255)
                 gTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGc                                             <  1:299510/69‑1 (MQ=255)
                 gTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGc                                             <  1:365872/69‑1 (MQ=255)
                 gTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGc                                             <  1:2684555/69‑1 (MQ=255)
                         cagccagcAGCGAATCAGCCCAGCCCTGAACCATCTGCGTCCAATAAACACACATAACAGCGGGacaa                                      >  1:945255/1‑68 (MQ=255)
                         cagccagcAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGacaa                                      >  1:821537/1‑68 (MQ=255)
                         cagccagcAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGacaa                                      >  1:1823988/1‑68 (MQ=255)
                         cagccagcAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGacaa                                      >  1:22232/1‑68 (MQ=255)
                         cagccagcAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGacaa                                      >  1:2281081/1‑68 (MQ=255)
                         cagccagcAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGacaa                                      >  1:295481/1‑68 (MQ=255)
                         cagccagcAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGAGGGacaa                                      >  1:2717976/1‑68 (MQ=255)
                          agccagcAGCGAATCAGCCCAGCCCTGAACCATCGGCGTc                                                                 >  1:2647563/1‑40 (MQ=255)
                          agccagcAGCGAATCAGCCCAGCCCTGAACCATCGGCGTc                                                                 >  1:3407467/1‑40 (MQ=255)
                          agccagcAGCGAATCAGCCCAGCCCTGAACCATCGGCGTc                                                                 >  1:970866/1‑40 (MQ=255)
                          agccagcAGCGAATCAGCCCAGCCCTGAACCATCGGCGTc                                                                 >  1:1978156/1‑40 (MQ=255)
                          agccagcAGCGAATCAGCCCAGCCCTGAACCATCGGCGTc                                                                 >  1:1694797/1‑40 (MQ=255)
                                cagcGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCa                              <  1:2793965/69‑1 (MQ=255)
                                cagcGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCa                              <  1:1581880/69‑1 (MQ=255)
                                 agcGAATCAGCCCAGCCCTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCaa                             <  1:1551874/69‑1 (MQ=255)
                                 agcGAATCAGCCCAGCCCTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCaa                             <  1:2759016/69‑1 (MQ=255)
                                 agcGAATCAGCCCAGCCCTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCaa                             <  1:2913458/69‑1 (MQ=255)
                                 agcGAATCAGCCCAGCCCTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCaa                             <  1:3410138/69‑1 (MQ=255)
                                 agcGAATCAGCCCAGCCCTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCaa                             <  1:1507517/69‑1 (MQ=255)
                                 agcGAATCAGCCCAGCCCTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCaa                             <  1:3178106/69‑1 (MQ=255)
                                 agcGAATCAGCCCAGCCCTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCaa                             <  1:640676/69‑1 (MQ=255)
                                 agcGAATCAGCCCAGCCCTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCaa                             <  1:2271315/69‑1 (MQ=255)
                                 agcGAATCAGCCCAGCCCTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCaa                             <  1:1356096/69‑1 (MQ=255)
                                 agcGAATCAGCCCAGCCCTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCaa                             <  1:253660/69‑1 (MQ=255)
                                 agcGAATCAGCCCAGCCCTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCaa                             <  1:241006/69‑1 (MQ=255)
                                  gcGAATCAGCCCAGCCCTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCaa                             <  1:400782/68‑1 (MQ=255)
                                   cGAATCAGCCCAGCCCCGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCaaa                            <  1:1151339/68‑1 (MQ=255)
                                           cccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCg                    <  1:1322520/68‑1 (MQ=255)
                                           cccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCg                    <  1:1318531/68‑1 (MQ=255)
                                           cccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCg                    <  1:3187108/68‑1 (MQ=255)
                                           cccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCg                    <  1:1500659/68‑1 (MQ=255)
                                           cccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCg                    <  1:3305325/68‑1 (MQ=255)
                                           cccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCg                    <  1:3348888/68‑1 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  <  1:1550844/69‑1 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  <  1:3274121/69‑1 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  <  1:396411/69‑1 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  <  1:1944302/69‑1 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  <  1:2760963/69‑1 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  >  1:9342/1‑69 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  >  1:1808454/1‑69 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  <  1:809405/69‑1 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  >  1:660295/1‑69 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  <  1:17616/69‑1 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  >  1:638704/1‑69 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  >  1:2335183/1‑69 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  <  1:2578229/69‑1 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  <  1:117119/69‑1 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  <  1:1236808/69‑1 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  <  1:1634465/69‑1 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  <  1:1609668/69‑1 (MQ=255)
                                            ccagcccTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGaa                  >  1:448348/1‑69 (MQ=255)
                                              agcgcTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACATACGCAAACACCAGCGAATg                <  1:381270/69‑1 (MQ=255)
                                              agcccTGAACCATCGACGTCCAATCAGGACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATg                <  1:300607/69‑1 (MQ=255)
                                              agcccTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATg                <  1:1057201/69‑1 (MQ=255)
                                              agcccTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATg                <  1:1008453/69‑1 (MQ=255)
                                              agcccTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATg                <  1:1386381/69‑1 (MQ=255)
                                              agcccTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATg                <  1:326118/69‑1 (MQ=255)
                                              agcccTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATg                <  1:3052067/69‑1 (MQ=255)
                                              agcccTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATg                <  1:3050557/69‑1 (MQ=255)
                                              agcccTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATg                <  1:2117549/69‑1 (MQ=255)
                                              agcccTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATg                <  1:2898222/69‑1 (MQ=255)
                                              agcccTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATg                <  1:1663231/69‑1 (MQ=255)
                                              agcccTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATg                <  1:2176844/69‑1 (MQ=255)
                                              agcccTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATg                <  1:1704614/69‑1 (MQ=255)
                                                 ccTGAACCATCGGCGTCCAATAAACACACATAACTGc                                             <  1:388916/37‑1 (MQ=255)
                                                  cTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCg                    <  1:2259550/61‑1 (MQ=255)
                                                  cTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTa            >  1:2052772/1‑69 (MQ=255)
                                                  cTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTa            >  1:1409956/1‑69 (MQ=255)
                                                  cTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTa            >  1:1143590/1‑69 (MQ=255)
                                                  cTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTa            >  1:2153381/1‑69 (MQ=255)
                                                  cTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTa            >  1:635614/1‑69 (MQ=255)
                                                  cTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTa            >  1:2764792/1‑69 (MQ=255)
                                                  cTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTa            >  1:2782062/1‑69 (MQ=255)
                                                  cTGAACCATCGACGTCCAATCAGCACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTa            >  1:2986999/1‑69 (MQ=255)
                                                    gAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTaaa          <  1:2603722/69‑1 (MQ=255)
                                                    gAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTaaa          <  1:1502922/69‑1 (MQ=255)
                                                            ggCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGAt  <  1:2210913/69‑1 (MQ=255)
                                                            ggCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGAt  <  1:587569/69‑1 (MQ=255)
                                                            ggCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGAt  <  1:1892025/69‑1 (MQ=255)
                                                             |                                                                   
GCGATAGCAACGAGACGGTTAAAAACAGCCAGCAGCGAATCAGCCCAGCCCTGAACCATCGGCGTCCAATAAACACACATAACAGCGGGACAACAAACGCAAACACCAGCGAATGGGTAAACCCGCGAT  >  minE/355331‑355459

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: