Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,275,342 T→G 17.4% E468A (GAG→GCG)  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,275,3420TG17.4% 39.5 / 6.3 23E468A (GAG→GCG) cioAcyanide insensitive ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base T (11/8);  new base G (2/2);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-01

GGAAGAGCATGCCGATCCCCAGGTCGAATCCGTCCATCACCACGTACATCATGACGCCGAAGATGATGATCACGGCCCAGATCAGTGGAAGGTCGATACCCATGGTTCAGTTCCCCTCGTTCAGGCTGGCAGAGGTGGCCTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCACCCTCACCGGTTTGCGGCCCTTTGCGCACCAGGCGCATCATGTAGCCAAGGCCGGTGCCAAACAGGGCGAA  >  NC_002947/5275202‑5275490
                                                                                                                                            |                                                                                                                                                    
ggAAGAGCATGCCGATCCCCAGGTCGAATCCGTCCATCACCACGTACATCATGACGCCGAAGATGATGATCACGGCCCAGATCAGTGGAAGGTCGATACCCATGGTTCAGTTCCCCTCGTTCAGGCTGGCAGAGGTGGCCTCGCGGCcat                                                                                                                                             >  1:454851/1‑150 (MQ=255)
                                        cacGTACATCATGACGCCGAAGATGATGATCACGGCCCAGATCAGTGGAAGGTCGATACCCATGGTTCAGTTCCCCTCGTTCAGGCTGGCAGAGGTGGCCTCGCGGCCAGCGGCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTg                                                                                                     >  1:468634/1‑150 (MQ=255)
                                             aCATCATGACGCCGAAGATGATGATCACGGCCCAGATCAGTGGAAGGTCGATACCCATGGTTCAGTTCCCCTCGTTCAGGCTGGCAGAGGTGGCCTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTg                                                                                                <  2:454851/150‑1 (MQ=255)
                                             aCATCATGACGCCGAAGATGATGATCACGGCCCAGATCAGTGGAAGGTCGATACCCATGGTTCAGTTCCCCTCGTTCAGGCTGGCAGAGGTGGCCTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTg                                                                                                <  2:39991/150‑1 (MQ=255)
                                                          gAAGATGATGATCACGGCCCAGATCAGTGGAAGGTCGATACCCATGGTTCAGTTCCCCTCGTTCAGGCTGGCAGAGGTGGCCGCGAGGcc                                                                                                                                               <  1:111960/90‑1 (MQ=255)
                                                          gAAGATGATGATCACGGCCCAGATCAGTGGAAGGTCGATACCCATGGTTCAGTTCCCCTCGTTCAGGCTGGCAGAGGTGGCCGCGAGGcc                                                                                                                                               >  2:111960/1‑90 (MQ=255)
                                                                    atCACGGCCCAGATCAGTGGAAGGTCGATACCCATGGTTCAGTTCCCCTCGTTCAGGCTGGCAGAGGTGGCCTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCa                                                                         >  2:243954/1‑150 (MQ=255)
                                                                      cACGGCCCAGATCAGTGGAAGGTCGATACCCATGGTTCAGTTGCCCTCGTTCAGGCTGGCAGAGGTGGCCTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCAcc                                                                       >  2:215580/1‑150 (MQ=255)
                                                                         ggCCCAGATCAGTGGAAGGTCGATACCCATGGTTCAGTTCCCCTCGTTCAGGCTGGCAGAGGTGGCCTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCACCCTc                                                                    >  2:142750/1‑150 (MQ=255)
                                                                                 tCAGTGGAAGGTCGATACCCATGGTTCAGTTCCCCTCGTTCAGGCTGGCAGAGGTGGCCTCGCgg                                                                                                                                                 >  2:73142/1‑65 (MQ=255)
                                                                                 tCAGTGGAAGGTCGATACCCATGGTTCAGTTCCCCTCGTTCAGGCTGGCAGAGGTGGCCTCGCgg                                                                                                                                                 <  1:73142/65‑1 (MQ=255)
                                                                                       gAAGGTCGATACCCATGGTTCAGTTCCCCTCGTTCAGGCTGGCAGAGGTGGCCTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCACCCTCACCGGTTTGCGGcc                                                      <  1:142750/150‑1 (MQ=255)
                                                                                                      tGGTTCAGTTCCCCTCGTTCAGGCTGGCAGAGGTGGCCTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCACCCTCACCGGTTTGCGGCCCTTTGCGCACCAGgc                                       <  1:243954/150‑1 (MQ=255)
                                                                                                         ttCAGTTCCCCTCGTTCAGGCTGGCAGAGGTGGCCTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCACCCTCACCGGTTTGCGGCCCTTTGCGCACCAGGCGca                                    >  2:603826/1‑150 (MQ=255)
                                                                                                               tCCCCTCGTTCAGGCTGGCAGAGGTGGCCTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCACCCTCACCGGTTTGCGGCCCTTTGCGCACCAGGCGCATCATGt                              <  1:364136/150‑1 (MQ=255)
                                                                                                                      gTTCAGGCTGGCAGAGGTGGCCTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCACCCTCACCGGTTTGCGGCCCTTTGCGCACCAGGCGCATCATGTAGCCAAg                       >  1:319376/1‑150 (MQ=255)
                                                                                                                           ggCTGGCAGAGGTGGCCTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCACCCTCACCGGTTTGCGGCCCTTTGCGCACCAGGCGCATCATGTAGCCAAGGCCgg                  <  1:603826/150‑1 (MQ=255)
                                                                                                                             cTGGCAGAGGTGGCCTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCACCCTCACCGGTTTGCGGCCCTTTGCGCACCAGGCGCATCATGTAGCCAAGGCCGGTg                <  2:346266/150‑1 (MQ=255)
                                                                                                                             cTGGCAGAGGTGGCCTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCACCCTCACCGGTTTGCGGCCCTTTGCGCACCAGGCGCATCATGTAGCCAAGGCCGGTg                >  1:205129/1‑150 (MQ=255)
                                                                                                                                     gatGGCCGCGAGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCACCCTCACCGGTTTGCGGCCCTTTGCGCACCAGGCGcatcat                                >  1:283029/3‑126 (MQ=255)
                                                                                                                                     gatGGCCGCGAGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCACCCTCACCGGTTTGCGGCCCTTTGCGCACCAGGCGcatcat                                <  2:283029/124‑1 (MQ=255)
                                                                                                                                           cTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCACCCTCACCGGTTTGCGGCCCTTTGCGCACCAGGCGCATCATGTAGCCAAGGCCGGTGCCAAACAGGGCGaa  >  1:363020/1‑150 (MQ=255)
                                                                                                                                           cTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCACCCTCACCGGTTTGCGGCCCTTTGCGCACCAGGCGCATCATGTAGCCAAGGCCGGTGCCAAACAGGGCGaa  >  1:19470/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
GGAAGAGCATGCCGATCCCCAGGTCGAATCCGTCCATCACCACGTACATCATGACGCCGAAGATGATGATCACGGCCCAGATCAGTGGAAGGTCGATACCCATGGTTCAGTTCCCCTCGTTCAGGCTGGCAGAGGTGGCCTCGCGGCCATCATCGGCGGCGGACAAAGGCCGTGCCGGCGTGCGTTTCTGGCCTGGGCCTCCCGGTGTGTGCTCGTCACCCTCACCGGTTTGCGGCCCTTTGCGCACCAGGCGCATCATGTAGCCAAGGCCGGTGCCAAACAGGGCGAA  >  NC_002947/5275202‑5275490

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: