Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,844 C→G 28.1% 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,8440CG28.1% 5.5 / 6.6 21A642A (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 (12/2);  new base G (3/3);  total (15/6)
Fisher's exact test for biased strand distribution p-value = 1.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

AGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGC  >  NC_000913/1708723‑1708993
                                                                                                                         |                                                                                                                                                     
aGCTGGAACAGCAACAGGCTAATGAGGAACCAGAAGAACAGGTCGATCCGCGCAACGCCGCCGTCCACGCCGCGATTGACCCGGCCACAGCCCGCAAGCGGGAGGGGCAACACACTAAGGCCGGGCGCGAACGTCAGGTCGGTCCGCGCa                                                                                                                           >  2:689804/1‑150 (MQ=9)
           cAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGACGCCGCTATTGCCCGAGCAAAAGCACTCAAGCTCGAACAGCAACAGGATAATGACTAGACAGAAGAACAGGTGGATCCGCGCAAAGCCGCAGta                                                                                                                >  2:561763/1‑149 (MQ=11)
                            aCCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCGGGAACAGCAACAGGCTAATGCCGAGCCAGAAGAACAGGTAGATCCGAGCAAAGCAGCCGTCGCAGCCGCTATTGcccc                                                                                               >  1:1725348/1‑149 (MQ=11)
                              cAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTGGAAGCGGCTATTGCCCGAGCCAAAGCCCGCAACCAGGACCAGCACCAGGCTATTGCCGTGGCAGAAGAAGAGGTCGATCAGCGCACAGCCCCATTTGACGCGGCTAGTGccccgg                                                                                             >  2:353689/1‑147 (MQ=2)
                                gaagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGCCGCCGGTATTGCCCGAGCCAAGGCCCGCAAGCGGGCACAGCACCAGGCTGATGCCGAGCCAGAACAACTGGTCGGTCCGCGCAACGCCGCCGTAGCAGCGGCTAATGCCCGCGcc                                                                                           >  1:741464/1‑150 (MQ=14)
                                gaagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAGAACAGGTCGATCAGCGCAAAGCCGCGGTAGACGCGGCGATTGCCCGCGcc                                                                                           >  1:977737/1‑150 (MQ=11)
                                    aaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCAAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAGAACAGGTGGATCCGCGCAAAGCTGCAGTAGAAGCCGCTATTGCAAGAGCAAAAg                                                                                       >  2:1412892/1‑150 (MQ=11)
                                                   gcAAAGCCGCCCTCGAAGCCGCTATAGACCGCACCAAATCACCCAATCTGGAACAGGAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAAc                                                                        <  1:1646700/150‑1 (MQ=17)
                                                      ccgccgccGTCGTAGCGGCTATTGCCCGTGCCGAAGCCCGCAAGCTGGAAGAGCAACAGGCTCATGCCGAGCCAGAACACCAGGTCGATCCGCGCAAATCCGCCGGCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGc                                                                     <  1:274420/148‑1 (MQ=255)
                                                       aGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAAGAACAGCGCTATCCCCCCAGACCCGCCGTCGAACCCGCATATGCCGGCGCGGAAGCCCGCAACAGGGaaaaaca                                                                    >  1:1618873/1‑145 (MQ=2)
                                                                        cTATTGCCCGTTCCACAGCACGCAAGCTCTAACAGCAACAGGCTAATTCAGATCCAGAAGACCCTTTAGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGc                                                   <  1:544633/150‑1 (MQ=14)
                                                                             gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCGGCCGTAGCCGCGGCTATTGCCCGCGCAACCGCGCGAAAACGCGCACAGCAGCAGGCTGAAGCGGAGccctaa                                              >  1:1208976/1‑146 (MQ=14)
                                                                                       aaGCCCGCAATATGTAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCTATCTGCTAAAATCCTCCGGCGTAGCCGCTATTGCCCGAGCCAAAGCGCGCAATCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTg                                    <  2:1568298/150‑1 (MQ=11)
                                                                                              cAAGCTGGAACAGCAACAGGCTAATGCGGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAAc                                                                        <  1:1726236/107‑1 (MQ=255)
                                                                                              cAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAAc                                                                        >  2:1726236/1‑107 (MQ=12)
                                                                                              cAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCGCGCAAAATCGCTCAACCGTCAACATCAAAAGACTAATGAGGAGACTAAAGAACAGTTTTATCcgaa                             >  2:666627/1‑148 (MQ=9)
                                                                                                aGCTGGAACAGCAACAGGCTAATGCGGTACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCGAACGCGCGCAAACCGGAACAGCAGCCGGCTAATGCGGAGCCCGAACATCAGGTGGATCAGCGCa                           >  1:364951/1‑150 (MQ=11)
                                                                                                          gCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGc                                                              <  1:73137/105‑1 (MQ=255)
                                                                                                          gCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGc                                                              >  2:73137/1‑105 (MQ=255)
                                                                                                                  cTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCCGTCGACGCgg         >  2:1500492/1‑150 (MQ=255)
                                                                                                                       gCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGATCAAACGCGCAAACCAGCCATCGAAGCGGCTAta    >  1:498512/1‑149 (MQ=14)
                                                                                                                         cGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAATCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGc  <  2:571847/150‑1 (MQ=255)
                                                                                                                         |                                                                                                                                                     
AGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGC  >  NC_000913/1708723‑1708993

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