Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 696,779 C→T 26.3% 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,7790CT26.3% 59.4 / 22.4 38P262L (CCC→CTC) appBcytochrome bd‑II oxidase, subunit II
Reads supporting (aligned to +/- strand):  ref base C (17/11);  new base T (2/8);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 6.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01

AATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTG  >  minE/696715‑696847
                                                                |                                                                    
aaTCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCTCGtttt                                                                <  1:1755544/71‑1 (MQ=255)
   ccGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCTCGTTTtgtg                                                             >  1:1415494/1‑71 (MQ=255)
      ttAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGAt                                                          <  1:560085/71‑1 (MQ=255)
      ttAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGAt                                                          <  1:191080/71‑1 (MQ=255)
        aaTGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCt                                                        >  1:111229/1‑71 (MQ=255)
          tGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGtt                                                                  <  1:983619/59‑1 (MQ=255)
              aCTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCgc                                                  >  1:1116840/1‑71 (MQ=255)
                tggtggCAGTGCTACCTTGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCgct                                                 <  1:1968826/70‑1 (MQ=255)
                  gtggCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCgctgct                                              <  1:1882876/71‑1 (MQ=255)
                         tGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCGCTGCTgg                                            <  1:154075/66‑1 (MQ=255)
                         tGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCGCTGCTgg                                            <  1:410895/66‑1 (MQ=255)
                         tGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCt                                                        >  1:566546/1‑54 (MQ=255)
                         tGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGAttc                                       <  1:909597/71‑1 (MQ=255)
                          gCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGAttct                                      >  1:330649/1‑71 (MQ=255)
                           cTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCGCTGCTGGGAttct                                      >  1:149723/1‑70 (MQ=255)
                           cTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGt                                                                   <  1:1237206/41‑1 (MQ=255)
                           cTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTtgtg                                                             >  1:236109/1‑47 (MQ=255)
                           cTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTtgtg                                                             >  1:177544/1‑47 (MQ=255)
                                 ggTGGCTGGATGAATAATTTTGTCGAGTCGCCCGttt                                                                 <  1:80364/37‑1 (MQ=255)
                                  gTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCAGCTGCTGGGATTCTTCTGCCCa                              >  1:418860/1‑71 (MQ=255)
                                   tGCCTGGATGAATAATTTTGTCGAGTCGCTCGtttt                                                                <  1:358514/36‑1 (MQ=255)
                                      cTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGcc                                >  1:24533/1‑65 (MQ=255)
                                       tGGATGAATAATTTTGTCGAGGCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGAttcttc                                    <  1:499943/60‑1 (MQ=255)
                                         gATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCgct                                                 >  1:610023/1‑45 (MQ=255)
                                         gATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGa                                          >  1:847069/1‑52 (MQ=255)
                                         gATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGa                                          >  1:1359513/1‑52 (MQ=255)
                                         gATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTg                                  >  1:1681984/1‑60 (MQ=255)
                                          aTGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGa                       >  1:76834/1‑70 (MQ=255)
                                             aataatTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGt                   <  1:1985798/71‑1 (MQ=255)
                                             aataatTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGt                   <  1:1553532/71‑1 (MQ=255)
                                             aataatTTTGTCGAGTCGCCCGTTTTGTGGATCtt                                                       >  1:1172008/1‑35 (MQ=255)
                                             aataatTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTggg                                           >  1:1502084/1‑47 (MQ=255)
                                              ataatTTTGTCGAGTCGCTCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGt                   <  1:474674/70‑1 (MQ=255)
                                                         gAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTa       <  1:2022636/71‑1 (MQ=255)
                                                          aGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTAt      <  1:205264/71‑1 (MQ=255)
                                                            tCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCg    >  1:1390111/1‑71 (MQ=255)
                                                              gCCCGTGTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTg  >  1:1303052/1‑71 (MQ=255)
                                                              gCCAGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTg  <  1:189379/71‑1 (MQ=255)
                                                                |                                                                    
AATCCGTTAATGAAACTGGTGGCAGTGCTACCTGGTGCCTGGATGAATAATTTTGTCGAGTCGCCCGTTTTGTGGATCTTCCCGCTGCTGGGATTCTTCTGCCCATTGCTGACGGTGATGGCGATTTATCGTG  >  minE/696715‑696847

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: