Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,847 G→A 24.3% E643E (GAG→GAA 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,8470GA24.3% 11.2 / 5.9 22E643E (GAG→GAArsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base G (12/3);  new base A (3/2);  total (15/7)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-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: