Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 696,779 C→T 32.0% P262L (CCC→CTC)  appB → cytochrome bd‑II oxidase, subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE696,7790CT32.0% 64.5 / 41.8 50P262L (CCC→CTC) appBcytochrome bd‑II oxidase, subunit II
Reads supporting (aligned to +/- strand):  ref base C (21/13);  new base T (4/12);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 3.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGA  >  minE/696714‑696838
                                                                 |                                                           
cAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCTCGttt                                                        >  1:534274/1‑71 (MQ=255)
cAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGttt                                                        >  1:3736118/1‑71 (MQ=255)
         aaTGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATc                                                <  1:1072101/70‑1 (MQ=255)
         aaTGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATc                                                <  1:2951636/70‑1 (MQ=255)
         aaTGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCt                                               >  1:1902396/1‑71 (MQ=255)
                  ggtggCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCgctgc                                      >  1:803713/1‑71 (MQ=255)
                  ggtggCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCgctgc                                      >  1:647421/1‑71 (MQ=255)
                  ggtggCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCgctgc                                      >  1:3813715/1‑71 (MQ=255)
                  ggtggCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCAgctgc                                      >  1:3311174/1‑71 (MQ=255)
                    tggCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCGctgctg                                    <  1:3131834/71‑1 (MQ=255)
                    tggCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCGctgctg                                    <  1:4068272/71‑1 (MQ=255)
                        aGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGAt                                >  1:2897874/1‑71 (MQ=255)
                         gTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCTCGTTTTGTgg                                                   <  1:336652/51‑1 (MQ=255)
                         gTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCTCGTTTTGTgg                                                   <  1:3053666/51‑1 (MQ=255)
                         gTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCTCGTTTTGTgg                                                   <  1:2440492/51‑1 (MQ=255)
                           gCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGAttct                             >  1:3340159/1‑71 (MQ=255)
                           gCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGAttct                             >  1:3142521/1‑71 (MQ=255)
                           gCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGAttct                             >  1:3152997/1‑71 (MQ=255)
                            cTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCTCCCGTTTTGTGGATCTTCCCGCTGCTGGGAttct                             <  1:914322/70‑1 (MQ=255)
                            cTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGAttctt                            >  1:1386705/1‑71 (MQ=255)
                            cTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGAttctt                            >  1:451505/1‑71 (MQ=255)
                               ccTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCgctgc                                      >  1:2537145/1‑58 (MQ=255)
                               ccTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTg                         >  1:2500965/1‑71 (MQ=255)
                               ccTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTg                         >  1:1903010/1‑71 (MQ=255)
                               ccTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTg                         >  1:1312747/1‑71 (MQ=255)
                                 tGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTtg                                                      >  1:1146965/1‑40 (MQ=255)
                                  ggTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGcc                       >  1:50426/1‑70 (MQ=255)
                                   gTGCCTGGATGAATAATTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGcc                       <  1:3163249/69‑1 (MQ=255)
                                   gTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGccc                      <  1:330891/70‑1 (MQ=255)
                                   gTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGccc                      <  1:3733673/70‑1 (MQ=255)
                                   gTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGccc                      <  1:2294894/70‑1 (MQ=255)
                                   gTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCa                     >  1:2723198/1‑71 (MQ=255)
                                    tGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGccc                      <  1:3150156/69‑1 (MQ=255)
                                    tGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCAt                    <  1:2673520/71‑1 (MQ=255)
                                      ccTGGATGAATAATTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTg                  <  1:385774/71‑1 (MQ=255)
                                      ccTGGATGAATAATTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTg                  <  1:747000/71‑1 (MQ=255)
                                        tGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTgg                                                   >  1:2153649/1‑36 (MQ=255)
                                        tGGATGAATAATTTTGTCGAGTCCCCCGTTTTGTgg                                                   >  1:3966245/1‑36 (MQ=255)
                                          gATGAATAATTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCt                <  1:2472040/69‑1 (MQ=255)
                                          gATGAATAATTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCt                <  1:1472773/69‑1 (MQ=255)
                                              aataatTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGt          <  1:3504232/71‑1 (MQ=255)
                                              aataatTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGt          <  1:475390/71‑1 (MQ=255)
                                              aataatTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGt          <  1:1625947/71‑1 (MQ=255)
                                                taatTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCt                <  1:3172848/63‑1 (MQ=255)
                                                  atTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATg      >  1:547371/1‑71 (MQ=255)
                                                     ttGTCGAGTCGCTCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGc                 >  1:2141396/1‑57 (MQ=255)
                                                       gTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGa  <  1:12716/70‑1 (MQ=255)
                                                          gAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGAttcttc                           >  1:1420662/1‑42 (MQ=255)
                                                              cGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGt          <  1:3790078/55‑1 (MQ=255)
                                                              cGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGt          <  1:1744993/55‑1 (MQ=255)
                                                                 |                                                           
CAATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGA  >  minE/696714‑696838

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: