Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,854 C→G 22.0% Q646E (CAA→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,8540CG22.0% 41.6 / 12.9 27Q646E (CAA→GAA) rsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base C (12/9);  new base G (3/3);  total (15/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 = 1.00e+00

CCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCG  >  NC_000913/1708711‑1708998
                                                                                                                                               |                                                                                                                                                
ccAAACCCCGCAATCTTTAAAAGCAACCGGCTAATCCGTAACCACCACACCAGCTCCATCCCTGCCAATCCCCCTTCGAATCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAg                                                                                                                                             <  4:380861/149‑1 (MQ=2)
                 gAACAGCAACAGGCTAATGCGGAACCATAAGAACAGTTCGATCCGCTCAAAGCCGCCTTCGCTGCCGCTATTGCCCCTGCCACAGCACGCTAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGt                                                                                                                                           <  7:358536/134‑1 (MQ=14)
                               tAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc                                                                                                               >  6:139660/1‑148 (MQ=255)
                                                  cAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgc                                                                                           <  5:139660/149‑1 (MQ=255)
                                                     gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaa                                                                                                                                                <  4:239649/93‑1 (MQ=255)
                                                                                          cccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCTCCCGGCTAGTGCGGGGCCag                                                     >  6:354043/1‑147 (MQ=17)
                                                                                              tGCCCAACCACGCCAGCTGGACACCCAACAGCCTAACGCAGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAAc                                               <  3:155718/149‑1 (MQ=14)
                                                                                                 cAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGa                                                                               >  5:412554/1‑114 (MQ=255)
                                                                                                 cAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGa                                                                               <  6:412554/114‑1 (MQ=255)
                                                                                                    aGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                                                                   <  2:419932/91‑1 (MQ=255)
                                                                                                    aGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGa                                                                                                   >  1:419932/1‑91 (MQ=255)
                                                                                                        cGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGAGCAAACGGGCTCCTCAACCGGCTAATGCGGAGCCAGAAGAACAGGATGagcg                                     >  8:387764/1‑146 (MQ=14)
                                                                                                                  aaCAGCAACAGGCTAATGCGGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGc                                                                          <  8:235584/102‑1 (MQ=255)
                                                                                                                     aGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGACGCCGCTATTGCCCGCGCACAAGCCAGCAAACGCGACCAGCACCAGGATAATTCCGAGCAAGAAGAACACGTTCAGCACAGCAAAGCTgac                        >  1:135264/1‑147 (MQ=9)
                                                                                                                            ggCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                                                            >  8:130103/1‑90 (MQ=11)
                                                                                                                            ggCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                                                            <  7:130103/90‑1 (MQ=11)
                                                                                                                                aaTGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                                                            <  8:118008/86‑1 (MQ=255)
                                                                                                                                aaTGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                                                            >  7:118008/1‑86 (MQ=255)
                                                                                                                                  tGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAACCGGGAACAGCAACCGGCTAATGCGGAGCAAGAACAACAGGGTGATCACCGCAAACATGCCGGCGAAGCgga             >  3:66394/1‑146 (MQ=17)
                                                                                                                                     cGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAACCGGGAACAGCAACAGCGATATGCCGAGCCACAAGAACAGGTGCATCCGCCCAACTCTCCCGACGAAGCCGCGAtat        >  5:85253/1‑147 (MQ=11)
                                                                                                                                        accAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGc                                                                   >  1:115915/2‑87 (MQ=255)
                                                                                                                                        accAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGc                                                                   <  2:115915/86‑1 (MQ=255)
                                                                                                                                          cAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                                                            >  6:18749/1‑76 (MQ=255)
                                                                                                                                          cAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACa                                                                            <  5:18749/76‑1 (MQ=255)
                                                                                                                                          cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAAc                                                                       >  7:82931/1‑81 (MQ=255)
                                                                                                                                          cAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAAc                                                                       <  8:82931/81‑1 (MQ=255)
                                                                                                                                           aGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCACCAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCCGCGGCTATTGCAcgcg  >  1:269433/1‑149 (MQ=14)
                                                                                                                                            gAACAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGCGCAAAAGCCAAAAAACCCGAACAGCAACAGGAAAATCCGGAGCCAGAAGAACAGGTTAATCCGCGCAACGCGTCAGTAGAAGCGGCTATTGCCcgcg  >  2:388072/1‑148 (MQ=2)
                                                                                                                                             aacaacAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAGCGGGAACAGCAACCGGCTAATGCGGa                                                          >  5:336677/1‑91 (MQ=37)
                                                                                                                                             aacaacAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGa                                                          <  6:336677/91‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCG  >  NC_000913/1708711‑1708998

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: