Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 709,461 T→A 41.3% A138A (GCA→GCT yceJ ← predicted cytochrome b561

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE709,4610TA41.3% 19.4 / 43.9 46A138A (GCA→GCTyceJpredicted cytochrome b561
Reads supporting (aligned to +/- strand):  ref base T (10/17);  new base A (6/13);  total (16/30)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACC  >  minE/709393‑709529
                                                                    |                                                                    
aaTCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCt                                                                      >  1:1447252/1‑69 (MQ=255)
   ccGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGcc                                                                   >  1:1179173/1‑69 (MQ=255)
   ccGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGcc                                                                   >  1:2625744/1‑69 (MQ=255)
   ccGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGcc                                                                   >  1:702801/1‑69 (MQ=255)
          tAACGGACAGTACTACAACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTgcg                                                             <  1:958061/68‑1 (MQ=255)
          tAACGGACAGTACTACAACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTgcg                                                             <  1:1245504/68‑1 (MQ=255)
          tAACGGACAGTACTACAACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTgcg                                                             <  1:3068394/68‑1 (MQ=255)
          tAACGGACAGTACTACAACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTgcg                                                             <  1:3051126/68‑1 (MQ=255)
          tAACGGACAGTACTACAACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTgcg                                                             <  1:1102905/68‑1 (MQ=255)
          tAACGGACAGTACTACAACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTgcg                                                             <  1:2967476/68‑1 (MQ=255)
          tAACGGACAGTACTACAACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTgcg                                                             <  1:2751189/68‑1 (MQ=255)
          tAACGGACAGTACTACAACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTgcg                                                             <  1:2701216/68‑1 (MQ=255)
          tAACGGACAGTACTACAACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTgcg                                                             <  1:1822345/68‑1 (MQ=255)
          tAACGGACAGTACTACAACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTgcg                                                             <  1:2461190/68‑1 (MQ=255)
          tAACGGACAGTACTACAACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTgcg                                                             <  1:2443049/68‑1 (MQ=255)
          tAACGGACAGTACTACAACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTgcg                                                             <  1:2290649/68‑1 (MQ=255)
           aaCGGACAGTACTACAACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTgcg                                                             <  1:700334/67‑1 (MQ=255)
                  aGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTTCCTGTGCGCCAGCGTCg                                                    <  1:962457/69‑1 (MQ=255)
                  aGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCg                                                    <  1:74133/69‑1 (MQ=255)
                  aGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCg                                                    <  1:1731318/69‑1 (MQ=255)
                  aGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCg                                                    <  1:2934752/69‑1 (MQ=255)
                   gTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCg                                                    <  1:2029164/68‑1 (MQ=255)
                   gTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCg                                                    <  1:1298941/68‑1 (MQ=255)
                        acgacgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCATgg                                              >  1:353620/1‑69 (MQ=255)
                        acgacgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAgg                                              >  1:637707/1‑69 (MQ=255)
                        acgacgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAgg                                              >  1:1473962/1‑69 (MQ=255)
                        acgacgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAgg                                              >  1:575131/1‑69 (MQ=255)
                        acgacgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAgg                                              >  1:3279562/1‑69 (MQ=255)
                            cgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTGCGCCAGCGTCGGAGAGGATcg                                          >  1:2244246/1‑69 (MQ=255)
                            cgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTGCGCCAGCGTCGGAGAGGATcg                                          >  1:724456/1‑69 (MQ=255)
                            cgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTGCGCCAGCGTCGGAGAGGATcg                                          >  1:1348251/1‑69 (MQ=255)
                            cgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTGCGCCAGCGTCGGAGAGGATcg                                          >  1:623123/1‑69 (MQ=255)
                            cgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTGCGCCAGCGTCGGAGAGGATcg                                          >  1:3415028/1‑69 (MQ=255)
                            cgCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCAGCCTGTGCGCCAGCGTCGGAGAGGATcg                                          >  1:2678647/1‑69 (MQ=255)
                                       aGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCa                               <  1:126850/69‑1 (MQ=255)
                                                     ggCACCGGCAAAGTCTGCCTGTGCTCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAAc                 <  1:1690657/69‑1 (MQ=255)
                                                     ggCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAAc                 <  1:3091128/69‑1 (MQ=255)
                                                     ggCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAAc                 <  1:3425659/69‑1 (MQ=255)
                                                     ggCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAAc                 <  1:1414987/69‑1 (MQ=255)
                                                     ggCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAAc                 <  1:672129/69‑1 (MQ=255)
                                                     ggCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAAc                 <  1:2648935/69‑1 (MQ=255)
                                                     ggCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAAc                 <  1:2641088/69‑1 (MQ=255)
                                                     ggCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAAc                 <  1:2367157/69‑1 (MQ=255)
                                                     ggCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAAc                 <  1:787343/69‑1 (MQ=255)
                                                     ggCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAAc                 <  1:834216/69‑1 (MQ=255)
                                                      gCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCCGCCAAAAACg                >  1:1195024/1‑69 (MQ=255)
                                                      gCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACg                >  1:1851820/1‑69 (MQ=255)
                                                                    tGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAACCGCTGATCGGTTTAcc  >  1:14251/1‑69 (MQ=255)
                                                                    |                                                                    
AATCCGTGCATAACGGACAGTACTACGACGCTCCACGCGAGCCAAAAATGCAGGGCACCGGCAAAGTCTGCCTGTGCGCCAGCGTCGGCAAGGGTCGCGGGGACGTCAAACCAGCCAAAAACGCTGATCGGTTTACC  >  minE/709393‑709529

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