Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,662 C→G 16.9% Q582E (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,6620CG16.9% 54.9 / 6.2 29Q582E (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 (7/15);  new base G (3/2);  total (10/19)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCG  >  NC_000913/1708514‑1708774
                                                                                                                                                    |                                                                                                                
aGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGaac                                                                                                                  <  1:296726/149‑1 (MQ=255)
    cAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAacaaca                                                                                                              <  5:350001/149‑1 (MQ=255)
       gCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCCAGCAGCAACAGCAACAGGCTAATGCGCAACCAGAACAACAgga                                                                                                           >  4:359358/1‑148 (MQ=255)
        ccAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAg                                                                                                             >  8:261972/1‑146 (MQ=255)
        ccAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAg                                                                                                             <  7:261972/146‑1 (MQ=255)
             gCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTAGAACAGCAACAGGCTAATGCGGAACCAGAACACCAGGTCGATcc                                                                                                     >  1:92834/1‑149 (MQ=255)
                        aGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAgccg                                                                                          >  8:261176/1‑149 (MQ=255)
                        aGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAgccg                                                                                          <  4:186987/149‑1 (MQ=255)
                        aGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAgccg                                                                                          <  6:6807/149‑1 (MQ=255)
                               cagcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGa                                                                                   <  2:40634/149‑1 (MQ=255)
                                   aaaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGtt                                                                                                          >  8:291820/1‑121 (MQ=255)
                                   aaaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGtt                                                                                                          <  7:291820/122‑2 (MQ=18)
                                   aaaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATAGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTACTGCGGACGCAGAACAACAGGTCGATCCGCGAAAAGCCGCCGTCGAAGcc                                                                               >  3:282463/1‑149 (MQ=21)
                                           gacgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAt                                                                          <  8:104392/146‑1 (MQ=255)
                                           gacgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAt                                                                          >  7:104392/1‑146 (MQ=255)
                                               cagcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                   <  6:245152/149‑1 (MQ=255)
                                                     ccGTTGCTGATCCACGTAAAACTGCCGTTGAAGCTGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAAAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaaa                                                             <  3:359358/149‑1 (MQ=255)
                                                       gTTGCTGATCCATGTAAAACTGCCGTTGAAGCATCTATAGCCCGAGCCAAAGATCGCAATCTGGAACAGCATCAGGATAATGAGGAACCAGAACAACATGTCGATCCTCGCAAAGCTGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgc                                                           <  2:92834/149‑1 (MQ=17)
                                                          gCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGATAATGCGGAACCAGAGGAAAAGGTAGATCAGCGAAAAGAAGCCGTCGAAGCCGGTATTGCCCGAGCCAAAgc                                                           >  8:110166/1‑146 (MQ=17)
                                                                    gTAAAACTGCCGTTGCAGCAGCTATCTCCCGCGCCACAGCGCGCAAGCTGGAACAGCATCAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGTAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGaa                                              <  2:364444/149‑1 (MQ=17)
                                                                    gTAAAACTGCCGGTGCCGCAGCGATCGCCCGCGCCACCGCGCGCAAGCTGGAACAGCACCAGGCTACTGCGGAACCAGAACACCAGGCCGAGCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGaa                                              <  8:391790/149‑1 (MQ=17)
                                                                      aaaaCTGCAGCTGAAGAAGCTATTGACCGGGCCCAAGCGCGGAACCTGGAACAGCAACAGGCTAATCCGGAACCAGAACAACAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACa                                            <  1:265568/149‑1 (MQ=11)
                                                                             cctttGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAg                                     <  7:110166/146‑1 (MQ=255)
                                                                                   aaGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTaa                                <  7:261176/148‑1 (MQ=255)
                                                                                     gcagcTCTCACCCGCGCCAACCCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGCACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGc                             <  5:52841/149‑1 (MQ=14)
                                                                                            tCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgc                                                                                                   <  7:336725/72‑1 (MQ=17)
                                                                                            tCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgc                                                                                                   >  8:336725/1‑72 (MQ=17)
                                                                                                ccGCGACAAATCGCGGAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCAGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCag                     <  5:329938/146‑1 (MQ=255)
                                                                                                              gcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaa                                                    <  4:223921/101‑2 (MQ=14)
                                                                                                              gcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaa                                                    >  3:223921/1‑100 (MQ=14)
                                                                                                                aatcTGGAACAGCAACAGGAAACTTCGGAACCAGAAAAACAGGCCGATCCGCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgcg  <  7:140251/146‑1 (MQ=14)
                                                                                                                                                    |                                                                                                                
AGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCG  >  NC_000913/1708514‑1708774

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: