Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,276,067 C→T 15.4% M226I (ATG→ATA cioA ← cyanide insensitive ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,276,0670CT15.4% 52.7 / 4.3 26M226I (ATG→ATAcioAcyanide insensitive ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (13/9);  new base T (2/2);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

GAACAGGATCAGTGGAGTCGGCTCACCGGGCACGTTCTCCCAGTGGCCTTCGATCGCAGCAATTTTCACCGGCTGGTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAGGCACCGACGAAGAACGCGGTGGACACGAACGCGGCAGTAGCCATGTGCATCAGGCGGTAGGGGAAGGAGGG  >  NC_002947/5275926‑5276213
                                                                                                                                             |                                                                                                                                                  
gAACAGGATCAGTGGAGTCGGCTCACCGGGCACGTTCTCCCAGTGGCCTTCGATCGCAGCAATTTTCACCGGCTGGTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACAt                                                                                                                                            >  2:838338/1‑150 (MQ=255)
      gATCAGTGGAGTCGGCTCACCGGGCACGTTCTCCCAGTGGCCTTCGATCGCAGCAATTTTCACCGGCTGGTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCAt                                                                                                                                      >  1:337088/1‑150 (MQ=255)
      gATCAGTGGAGTCGGCTCACCGGGCACGTTCTCCCAGTGGCCTTCGATCGCAGCAATTTTCACCGGCTGGTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCAt                                                                                                                                      >  1:93879/1‑150 (MQ=255)
            tGGAGTCGGCTCACCGGGCACGTTCTCCCAGTGGCCTTCGATCGCAGCAATTTTCACCGGCTGGTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAg                                                                                                                                >  2:340265/1‑150 (MQ=255)
                    gCTCACCGGGCACGTTCTCCCAGTGGCCTTCGATCGCAGCAATTTTCACCGGCTGGTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGc                                                                                                                        <  2:337088/150‑1 (MQ=255)
                      tCACCGGGCACGTTCTCCCAGTGGCCTTCGATCGCAGCAATTTTCACCGGCTGGTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGt                                                                                                                      >  1:495688/1‑150 (MQ=255)
                                                ttCGATCGCAGCAATTTTCACCGGCTGGTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCgcag                                                                                            >  1:39244/1‑150 (MQ=255)
                                                     tCGCAGCAATTTTCACCGGCTGGTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGt                                                                                       <  1:838338/150‑1 (MQ=255)
                                                         agcaATTTTCACCGGCTGGTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCa                                                                                   >  2:785319/1‑150 (MQ=255)
                                                         agcaATTTTCACCGGCTGGTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTGACGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCa                                                                                   >  1:8500/1‑150 (MQ=255)
                                                              ttttCACCGGCGGGTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGgccg                                                                              >  2:33370/1‑150 (MQ=255)
                                                                   aCCGGCTGGTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAg                                                                         >  1:717364/1‑150 (MQ=255)
                                                                      ggCTGGTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGCTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAgtcc                                                                      <  2:495688/150‑4 (MQ=255)
                                                                          ggTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAGGCACcga                                                                  <  2:342320/150‑1 (MQ=255)
                                                                              cTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAGGCACCGACgaa                                                              <  1:785319/150‑1 (MQ=255)
                                                                                  aGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAGGCACCGACGAAGAAc                                                          >  2:682945/1‑150 (MQ=255)
                                                                                     gTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAGGCACCGACGAAGAAcgcg                                                       <  1:340265/150‑1 (MQ=255)
                                                                                           aGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAGGCACCGACGAAGAACGCGGTGGac                                                 <  1:682945/150‑1 (MQ=255)
                                                                                                       tCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAGGCACCGACGAAGAACGCGGTGGACACGAACGCGGCa                                     <  2:717364/150‑1 (MQ=255)
                                                                                                                        aaTCGGGGCGACGATCAATGCTATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAGGCACCGACgaa                                                              <  1:118724/108‑1 (MQ=255)
                                                                                                                        aaTCGGGGCGACGATCAATGCTATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAGGCACCGACgaa                                                              >  2:118724/1‑108 (MQ=255)
                                                                                                                        aaTCGGGGCGACGATCAATGCTATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAGGCACCGACgaa                                                              >  2:118789/1‑108 (MQ=255)
                                                                                                                        aaTCGGGGCGACGATCAATGCTATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAGGCACCGACgaa                                                              <  1:118789/108‑1 (MQ=255)
                                                                                                                          tCGTGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAGGCACCGACGAAGAACGCGGTGGACACGAACGCGGCAGTAGCCATGTGCATCAGGc                  <  1:33370/150‑1 (MQ=255)
                                                                                                                                         aTGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAGGCACCGACGAAGAACGCGGTGGACACGAACGCGGCAGTAGCCATGTGCATCAGGCGGTAGGGGAaggagg   >  1:709735/1‑150 (MQ=255)
                                                                                                                                          tGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAGGCACCGACGAAGAACGCGGTGGACACGAACGCGGCAGTAGCCATGTGCATCAGGCGGTAGGGGAAGGAggg  >  1:122418/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GAACAGGATCAGTGGAGTCGGCTCACCGGGCACGTTCTCCCAGTGGCCTTCGATCGCAGCAATTTTCACCGGCTGGTGCTTGAGGGTATTGAGGCCGTGGAAATCGCCGATCACTGCCTGAATCGGGGCGACGATCAATGCCATCCACATGGCCATCGACAGCATCTTGCGTACAGCCGGGTTGTCACGCCCGCGCAGCAGGTGCCAGGCCGCCGAGGCACCGACGAAGAACGCGGTGGACACGAACGCGGCAGTAGCCATGTGCATCAGGCGGTAGGGGAAGGAGGG  >  NC_002947/5275926‑5276213

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