Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,844 C→G 30.6% A642A (GCC→GCG 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,8440CG30.6% 7.6 / 3.2 13A642A (GCC→GCGrsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base G (2/2);  total (6/7)
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.83e-01

GCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGG  >  NC_000913/1708772‑1708956
                                                                        |                                                                                                                
gcgcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGCGCCAGAACAACAGGTCGATCCGCCGAAAGCCGACGTCGAAGCCGCAATAGCACTAGCAAAAGCGCGCAAACgg                                       >  8:171712/1‑148 (MQ=14)
     aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACATCAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAg                                 <  7:171712/149‑1 (MQ=255)
                gAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTAGATCCGAGCAAAGCAGCCGTCGAATCCGCTATTGCCCGAGCATAGGCGCGCACACGGTCAGCGCAACAGGCTaa                      >  7:132860/1‑149 (MQ=14)
                       cTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACCGGCTAATGCCGAGCCAGAAGAACAGGTCGGTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGCACAGCAACCGGCTAATGCggcg               >  1:22600/1‑147 (MQ=255)
                                 cgCCAAATCACGAAACCTGGAACAGCCACACCCTACAGCCGACCCAGCACACCCTGTCGATCCCAACAACCCCGCCGTCGCAGCAGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAAc     <  2:22600/148‑1 (MQ=255)
                                    taaataaCGAAAGCGGGAACAGCAACATCAAAATCCCGACCCACAACAATAGCTTCAACCTCGCAAATCCGAAGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAgg  <  8:132860/143‑1 (MQ=9)
                                    cAAAGCGCGCTAGCTGGAGCTGCATCACGCTAATGCGGTTGCTCTTCAAATGGTAGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAgg  <  1:13006/149‑1 (MQ=11)
                                            gCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa                                                                                                         <  1:119415/38‑1 (MQ=255)
                                            gCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaa                                                                                                         >  2:119415/1‑38 (MQ=255)
                                                           aaCAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGaa                                    >  1:72659/1‑92 (MQ=255)
                                                           aaCAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGaa                                    <  2:72659/92‑1 (MQ=255)
                                                                  tAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGaa                                    <  1:87923/85‑1 (MQ=255)
                                                                  tAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGaa                                    >  2:87923/1‑85 (MQ=255)
                                                                        |                                                                                                                
GCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGG  >  NC_000913/1708772‑1708956

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: