Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,655 A→G 19.7% E579E (GAA→GAG rsxC → SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,708,6550AG19.7% 32.5 / 3.2 23E579E (GAA→GAGrsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base A (7/10);  new base G (2/2);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGC  >  NC_000913/1708508‑1708741
                                                                                                                                                   |                                                                                      
ccGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAc                                                                                       >  7:239734/1‑149 (MQ=255)
 cGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAc                                                                                       >  4:149960/1‑148 (MQ=255)
  gTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCa                                                                                     <  5:172418/149‑1 (MQ=255)
     aaGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGaa                                                                                  <  3:149960/149‑1 (MQ=255)
       gcgcAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAAGGCGGACCCAGcac                                                                                 >  1:10396/1‑148 (MQ=255)
                  agCGAAACAGGCAGAACGGCAGCAAACTCCCGCCGCAGCCACCGTTGCTGAGCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGGCGATc                                                                     <  6:138723/149‑1 (MQ=25)
                                     cagcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTc                                                    <  8:179602/147‑1 (MQ=255)
                                         aaaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACGAGTACGACAGGTCGAGCCGCCGACAGACGCCGTcaacgag                                              >  3:161942/1‑143 (MQ=21)
                                          aaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCAAATGCGAAGCCCGCACAACAGGTCTATCCGCACAAAACCCCCGGCCAAACCg                                             >  7:108113/1‑149 (MQ=21)
                                                      agcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGACGAACCAGAACAACACTTCCATCCGCTCAACGCAGGCGTCGTAGCAGCGAGTGTCCGTg                                 >  1:136000/1‑149 (MQ=18)
                                                           cccctgatGATCCTCGAAAAACTGCCGTAGAAGCAGCTATTTCCCGGGACTAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaaa                            <  3:54250/142‑1 (MQ=11)
                                                                   gATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaaa                            <  1:104393/141‑1 (MQ=255)
                                                                   gATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGAGCAAGCTGGAAAAGAAACAGGCTAATGCGGATACAGAACAAAAGTTCGATCCGCGCTAAGACTCAGTAGAAGACGCTATTGCCCGTGCCaaa                            >  2:104393/1‑141 (MQ=255)
                                                                        aCGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGCACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTg                >  4:107535/1‑148 (MQ=255)
                                                                             aaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAacaaca                                                                             <  2:27859/82‑1 (MQ=14)
                                                                             aaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACa           <  3:107535/148‑1 (MQ=255)
                                                                             aaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAacaaca                                                                             >  1:27859/1‑82 (MQ=255)
                                                                              aaCTGCCGTTGAAGCAGCTATCGCCCGCGCCACATCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGc         <  8:110624/149‑1 (MQ=255)
                                                                                  ctcgatGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACa     <  3:91848/144‑1 (MQ=255)
                                                                                   ccGTCGAAGCAGCTAGCGCCCGCGCCAAAGCGCGCAAGCAGGAACACCAACAGGCTAATCCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAg    <  5:16695/149‑1 (MQ=14)
                                                                                     gTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATACTCTCAAAGCCGCCGTCGAAGCCGCTATTGCGCGTGCCAAAGCGAGCAAGCTGGGACGGCAACAgag  >  4:91848/1‑147 (MQ=14)
                                                                                      ttGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTg                                       <  7:69676/111‑1 (MQ=38)
                                                                                      ttGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTg                                       >  8:69676/1‑111 (MQ=38)
                                                                                                                                                   |                                                                                      
CCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGC  >  NC_000913/1708508‑1708741

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