Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,296,406 A→G 31.9% G555G (GGT→GGC katG ← catalase/hydroperoxidase HPI(I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,296,4060AG31.9% 59.3 / 37.7 75G555G (GGT→GGCkatGcatalase/hydroperoxidase HPI(I)
Reads supporting (aligned to +/- strand):  ref base A (24/27);  new base G (10/14);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

CCGTCAGCAATTGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTTCTCAACACCAACCACA  >  minE/2296338‑2296472
                                                                    |                                                                  
ccGTCAGCAATTGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGa                                                                    <  1:433191/69‑1 (MQ=255)
ccGTCAGCAATTGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGa                                                                    <  1:2447050/69‑1 (MQ=255)
ccGTCAGCAATTGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGa                                                                    <  1:149779/69‑1 (MQ=255)
ccGTCAGCAATTGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGa                                                                    <  1:2728784/69‑1 (MQ=255)
       cAATTGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGgc                                                              >  1:2158393/1‑68 (MQ=255)
       cAATTGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGgc                                                              >  1:2369012/1‑68 (MQ=255)
       cAATTGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGgc                                                              >  1:1872412/1‑68 (MQ=255)
       cAATTGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGgc                                                              >  1:1822170/1‑68 (MQ=255)
       cAATTGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGgc                                                              >  1:686972/1‑68 (MQ=255)
       cAATTGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGgc                                                              >  1:137669/1‑68 (MQ=255)
         aTTGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGGCGCa                                                           <  1:273754/69‑1 (MQ=255)
         aTTGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGGCGCa                                                           <  1:762956/69‑1 (MQ=255)
         aTTGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGGCGCa                                                           <  1:803238/69‑1 (MQ=255)
         aTTGGCTCCAGCAGCTCAAACATCTCAATGTCAGCCTGATCCTGACGCGCATCAACGCGACCCGGCGCa                                                           <  1:2314606/69‑1 (MQ=255)
          ttGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGACCTGGCGCa                                                           <  1:441465/68‑1 (MQ=255)
          ttGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGGCGCa                                                           <  1:1193638/68‑1 (MQ=255)
                      gcTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGGCGCAAACGGTACATGaa                                              <  1:1489996/69‑1 (MQ=255)
                                 tcAATGTCTGTCTGATCCTGACGCGCATCAACGCGACCCGGCGCAAACGGTACATg                                                <  1:680742/56‑1 (MQ=255)
                                   aaTGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGGCGCa                                                           >  1:2830020/1‑43 (MQ=255)
                                   aaTGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGGCGCa                                                           >  1:1074864/1‑43 (MQ=255)
                                   aaTGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGGCGCa                                                           >  1:385370/1‑43 (MQ=255)
                                   aaTGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGGCGCa                                                           >  1:1978040/1‑43 (MQ=255)
                                   aaTGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGGCGCa                                                           >  1:2300201/1‑43 (MQ=255)
                                       tCAGTCTGATCCTGACGCGCATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTgcggcg                             >  1:894973/1‑69 (MQ=255)
                                       tCAGTCTGATCCTGACGCGCATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTgcggcg                             >  1:2287279/1‑69 (MQ=255)
                                             tcaTCCTGACGCGCATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGc                        >  1:1117798/3‑68 (MQ=255)
                                             tGATCCTGACGCGCATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCg                       >  1:3033686/1‑69 (MQ=255)
                                             tGATCCTGACGCGCATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCg                       >  1:1741586/1‑69 (MQ=255)
                                                tCCTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGc                                           >  1:614155/1‑46 (MQ=255)
                                                tCCTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGc                                           >  1:3134530/1‑46 (MQ=255)
                                                tCCTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGc                                           >  1:3414982/1‑46 (MQ=255)
                                                tCCTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGc                                           >  1:690814/1‑46 (MQ=255)
                                                tCCTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGc                                           >  1:2613194/1‑46 (MQ=255)
                                                tCCTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGc                                           >  1:310717/1‑46 (MQ=255)
                                                tCCTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGc                                           >  1:870127/1‑46 (MQ=255)
                                                tCCTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGc                                           >  1:293760/1‑46 (MQ=255)
                                                tCCTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGc                                           >  1:906072/1‑46 (MQ=255)
                                                tCCTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGc                                           >  1:1103380/1‑46 (MQ=255)
                                                 cgtGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCtt                   <  1:2985457/67‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCtt                   <  1:2961765/69‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCtt                   <  1:3342105/69‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCtt                   <  1:3008953/69‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCtt                   <  1:516740/69‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCtt                   <  1:285820/69‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCtt                   <  1:1020425/69‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCtt                   <  1:2400823/69‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCtt                   <  1:224266/69‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCtt                   <  1:2100166/69‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCtt                   <  1:2045484/69‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCtt                   <  1:1927734/69‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCtt                   <  1:829789/69‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCt                    <  1:2808613/68‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCt                    <  1:994929/68‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCt                    <  1:1647578/68‑1 (MQ=255)
                                                 ccTGACGCGCATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCt                    <  1:1696277/68‑1 (MQ=255)
                                                  cTGACGCGCATCAACGCGGCCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCtt                   <  1:579049/68‑1 (MQ=255)
                                                  cTGACGCGCATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCt                    <  1:3241502/67‑1 (MQ=255)
                                                  cTGACGCGCATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCt                    <  1:187665/67‑1 (MQ=255)
                                                       gcgcATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTtctc               >  1:495223/1‑67 (MQ=255)
                                                       gcgcATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTTCTCaa             >  1:3013705/1‑69 (MQ=255)
                                                       gcgcATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTTCTCaa             >  1:665671/1‑69 (MQ=255)
                                                       gcgcATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTTCTCaa             >  1:2837961/1‑69 (MQ=255)
                                                       gcgcATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTTCTCaa             <  1:3022405/69‑1 (MQ=255)
                                                       gcgcATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTTCTCaa             >  1:2961763/1‑69 (MQ=255)
                                                              aaCGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCttt                  >  1:760459/1‑57 (MQ=255)
                                                              aaCGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTtctc               <  1:754371/60‑1 (MQ=255)
                                                              aaCGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTtctc               <  1:1073523/60‑1 (MQ=255)
                                                              aaCGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTtctc               <  1:1220481/60‑1 (MQ=255)
                                                              aaCGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTtctc               <  1:1359476/60‑1 (MQ=255)
                                                              aaCGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTtctc               <  1:2613097/60‑1 (MQ=255)
                                                              aaCGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTtctc               <  1:1664401/60‑1 (MQ=255)
                                                              aaCGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTtctc               <  1:2937849/60‑1 (MQ=255)
                                                              aaCGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTtctc               <  1:2145877/60‑1 (MQ=255)
                                                              aaCGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTTCTCAACaccaac      >  1:2111436/1‑69 (MQ=255)
                                                              aaCGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTTCTCAACaccaac      >  1:1237167/1‑69 (MQ=255)
                                                                  cgACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTTCTCAACACCAACcaca  <  1:3278934/69‑1 (MQ=255)
                                                                  cgACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTTCTCAACACCAACcaca  <  1:1137828/69‑1 (MQ=255)
                                                                  cgACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTTCTCAACACCAACcaca  <  1:112991/69‑1 (MQ=255)
                                                                  cgACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTTCTCAACACCAACcaca  <  1:2830964/69‑1 (MQ=255)
                                                                    |                                                                  
CCGTCAGCAATTGGCTCCAGCAGCTCAAACATCTCAATGTCAGTCTGATCCTGACGCGCATCAACGCGACCCGGCGCAAACGGTACATGAATGCTCAAACCTGCGGCGCTTGCGGCTTTCTCAACACCAACCACA  >  minE/2296338‑2296472

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: