Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,930 C→A 24.1% P671Q (CCG→CAG)  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,9300CA24.1% 44.3 / 2.6 30P671Q (CCG→CAG) rsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base C (15/7);  new base A (6/2);  total (21/9)
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

CCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGT  >  NC_000913/1708802‑1709074
                                                                                                                                |                                                                                                                                                
ccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGGAACGGGGGAGCGCACCAGGCTACTGCGGACCCaggag                                                                                                                              >  5:368122/1‑149 (MQ=14)
  gTGCCAAAGCCCGCAAGCCGGAACAGCAACAGGCTAATGCCGACCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAgaagaa                                                                                                                            <  7:64419/149‑1 (MQ=255)
                         cAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGc                                                                                                     >  7:65349/1‑149 (MQ=255)
                             aaCAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCACCCGGCTAATGCGGACCCAGAAGAACAGGGTGATCCGCGCAAAGATgactgc                                                                                                 >  1:21889/1‑144 (MQ=14)
                                       tGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCGGCCGTCGAAGCCGCTATTGCCCGCGCCCACGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGACGCGGCTa                                                                                       >  5:79387/1‑149 (MQ=14)
                                                acaacaacCGGTCGATCCGCGCAAATCCGCCGGCGTATCCGCTATTTCCCGATCCAATGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTAtt                                                                                     <  5:160854/140‑1 (MQ=21)
                                                aGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCCCGAAAACGGGAACAGCAACCGGCTAATGCGGACCCAGAAGAACAGGTTGATAAGCGCAAAGCTGCAGTGGCAGCGGATATTGCAagcg                                                                              >  8:193649/1‑149 (MQ=14)
                                                aGAACAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGCGCGAAAGCCCTCAAACCGGCACAGCAACAGGATAATGCGGAGCCATAAAAACAATTAAATCAGCGCAAAGCTACAGTGGCCTCGTCTATTGCAcgc                                                                               >  8:100758/1‑148 (MQ=2)
                                                aGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGGAAAAGGGAAAGGCAACAGGCAAATGAGGAGCAAGAAGAGCAGGTGGAACCGCGAAAAGATACCGTCGAAGCGGCGATTgagagcg                                                                              >  6:160854/1‑143 (MQ=14)
                                                  aacaacAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaa                                                                                                                               <  2:266449/98‑1 (MQ=255)
                                                  aacaacAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCCCGCAACCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCCCGAAAAGCAGCCGTCGAAGCGGCAATAGCAcgcac                                                                             >  7:266189/1‑146 (MQ=14)
                                                     aacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAACCAGAAGAACAGGTTGAGACGCGCAAAGCACCCGTCGGAGCTCCTCTTCCACGCGCCaat                                                                         >  5:68539/1‑148 (MQ=14)
                                                     aacaGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACAGGCTAATGCGGAGCAAGAATAACAGTTTGATCCGCGAAAAGCTGCCGTAGACGCGGCTATTGCCCGAGCaaat                                                                         >  8:104227/1‑148 (MQ=14)
                                                       caGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAAACGGCTAATGCGGAGCCAGAAGAACAGGTTGATCAGCGCAAAGCTGCAGTAGCCGCGGCTATTGCACGCGCAAAGGc                                                                       >  4:147394/1‑149 (MQ=17)
                                                          gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGt                                                                                                                       <  6:393016/98‑1 (MQ=255)
                                                          gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGt                                                                                                                       >  5:393016/1‑98 (MQ=255)
                                                          gTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCCGCGGCGATTGCCCGCGCCAAAGCACg                                                                    >  1:226164/1‑149 (MQ=255)
                                                             gATCCTCGCAACTCCGCCTTCGAATCCTTTATTGCCATAGCCAAATCGAGGATCTTGGTAGATCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCaa                                                                 <  2:21889/149‑1 (MQ=18)
                                                              aTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTaa                                                                                                                                            <  2:24531/73‑1 (MQ=255)
                                                              aTCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTaa                                                                                                                                            >  1:24531/1‑73 (MQ=255)
                                                                               gTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCAGTCGACGCGGCTATTGCCCGCGCACAAGCACTCAAGGGGGACCAGCAAAAgg                                                >  4:274266/1‑148 (MQ=14)
                                                                                       cGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTAtt                                                                                     >  7:405200/1‑103 (MQ=255)
                                                                                       cGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTAtt                                                                                     <  8:405200/103‑1 (MQ=255)
                                                                                        gCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGCAGCGGCTATTGCCCGCGCCAACGCCCGAAAGCGGGCACAGCACCAGGCTATTGCGGt                                      >  8:96467/1‑149 (MQ=21)
                                                                                               cccGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTAGACGCGGCTATTGCACGCGCCAATGCAAGAAAGCAGGACCAGCAAAAGGATAAAGAGGAACCaga                                >  7:312235/1‑148 (MQ=21)
                                                                                                       aaaGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTa                                                                                       <  7:319856/85‑1 (MQ=11)
                                                                                                                aaaCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAACGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCcgcg              >  5:9625/1‑149 (MQ=255)
                                                                                                                      gAACAGCAACAGGCTAATGCGGACCCAGAACAACAGGTAGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGAACGCGCCaaa                                                                         >  2:373044/1‑84 (MQ=11)
                                                                                                                           ccAACAGGCTAATGCGGAGCCAGTCTAACTCGTTGAGCCGCGCTAAGCTGCGGTCGAAGCGGCTATTGCACGAGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAg   <  3:361650/148‑1 (MQ=17)
                                                                                                                            cAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCCAGCAACCGCGAACAGCAACAGGCTAAAGCGGTACCAGAAGACCAGGTAAAAACGCCCAAAGCGGCAGt  >  4:2197/1‑149 (MQ=25)
                                                                                                                                |                                                                                                                                                
CCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGT  >  NC_000913/1708802‑1709074

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: