Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,662 C→G 14.3% 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,6620CG14.3% 110.4 / 5.9 41Q582E (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/22);  new base G (3/3);  total (15/26)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGAT  >  NC_000913/1708518‑1708769
                                                                                                                                                |                                                                                                           
cAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGaacaac                                                                                                          <  3:47630/148‑1 (MQ=255)
cAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGaacaac                                                                                                          <  2:317723/148‑1 (MQ=255)
  aGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAg                                                                                                        <  5:228804/148‑1 (MQ=255)
        agCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAACAACAGGTCGAt                                                                                                  >  5:363250/1‑148 (MQ=255)
         gCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATc                                                                                                 >  6:319643/1‑148 (MQ=255)
                      aaCTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAgccgc                                                                                    <  4:197288/148‑1 (MQ=255)
                        cTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGccgccg                                                                                  <  3:122690/148‑1 (MQ=255)
                          tcagcaAACTGCCGAAGCAGAACCCGTGGCTGATCCCCGTAAATATGCCGTTGAAGCAGCTATAGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTc                                                                                <  5:298061/147‑1 (MQ=14)
                           cagcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTc                                                                                >  1:17394/1‑147 (MQ=255)
                             gcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTACTGCGGAACCAGAACAACAGGTCGATCCGCGCAACGCCGCCGTCGAAg                                                                            >  3:248927/1‑148 (MQ=255)
                                 catcccgaccCCGCACCCGATGCTTATCCTCTCACACCCGCCGTCGAAGCAGCTATCGCCCTCGCCAAAGCGCGCAAGCTGCAACAGCAACAGGCTAATGCGGAACCAGAACACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCg                                                                         <  4:214420/143‑1 (MQ=11)
                                 aCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCg                                                                         <  4:149831/148‑1 (MQ=255)
                                   ttccgacgCAGAAACACTTGTTGATCCTCGAAAAACTGCAGTTGAAGCAGCTATTGCCCGTGACAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGTTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCt                                                                       <  3:315065/145‑1 (MQ=11)
                                   tAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGCGGAGACAGAGGAACAGGACGCTCCGCGCAACGCCGCAGTCGGAGCCGCt                                                                       >  1:197350/1‑148 (MQ=18)
                                     acgacgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCt                                                                       <  5:319643/146‑1 (MQ=255)
                                        acgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGACCAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCg                                                                         >  5:90375/1‑141 (MQ=255)
                                        acgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCg                                                                         <  6:90375/141‑1 (MQ=255)
                                        acgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTagagg                                                                  >  3:214420/1‑144 (MQ=255)
                                                 ccGTTGCTGATCCACGTAAAACTGCCGTTGAAGCATCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCaa                                                         <  8:251670/148‑1 (MQ=255)
                                                      gCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcg                                                     <  7:315197/147‑1 (MQ=255)
                                                           tCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                               <  6:63165/148‑1 (MQ=255)
                                                                gTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTc                                                                                >  6:325952/1‑110 (MQ=255)
                                                                gTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTc                                                                                <  5:325952/110‑1 (MQ=255)
                                                                  aaaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAgc                                                                                       >  5:268342/1‑101 (MQ=255)
                                                                  aaaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAgc                                                                                       <  6:268342/101‑1 (MQ=255)
                                                                   aaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGcc                                                                 <  1:148997/122‑1 (MQ=255)
                                                                   aaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGcc                                                                 >  2:148997/1‑122 (MQ=255)
                                                                   aaaCAGCCTTTGCAGCACCTATCGCCCTCGCCACACCGGGCAACCTGGCAGACCATCAGCATAATCCGGACGCAGCAAAACAGGCCGATCCTCGCAAACCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAc                                        <  6:363250/147‑1 (MQ=16)
                                                                    aaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAACCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAg                                      <  5:99176/148‑1 (MQ=255)
                                                                         ccGTTGACTCAGCTATCTCCCGCGCCACAGCGCGCAAGCTGGAACCGAAACATGCTAATGCTGAACCCTAACAACTGGGCTATCCGCGGAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCaa                                   <  4:248927/146‑1 (MQ=14)
                                                                           gTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAAGAACAGGta                                                                                                     >  7:364658/1‑77 (MQ=255)
                                                                             agccgcagcTATCTAACGCGCCAAATCTCGGAATCAGGAACACCTACAGGCTACTGAGGAACCAGCACCACATGCCGATCCGCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCt                             <  4:89519/144‑1 (MQ=14)
                                                                              gcagcagcTTTCTCCCTCGCCCAAGCGCGCAAGCTGGAACTCCATCAGGCTATTGCGGACCCAGCACAACAGGTCCTACCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTa                            <  3:199122/146‑1 (MQ=14)
                                                                                      tATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgccg                                                                                     >  3:310165/1‑83 (MQ=18)
                                                                                      tATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgccg                                                                                     <  4:310165/83‑1 (MQ=18)
                                                                                           cccGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCaga               <  3:350392/148‑1 (MQ=255)
                                                                                             cgcgCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGaa                                                                             <  3:263032/84‑1 (MQ=255)
                                                                                             cgcgCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGaa                                                                             >  4:263032/1‑84 (MQ=38)
                                                                                                        gctcAAGCTTGAACAGCATCAGGCTAATTCCGAACCAGAACAACAGGTCGATCCTCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGAt  <  2:17394/145‑1 (MQ=14)
                                                                                                               cTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGc                                                            <  5:345384/83‑1 (MQ=33)
                                                                                                               cTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGc                                                            >  6:345384/1‑83 (MQ=33)
                                                                                                                     cAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgc                                                    >  7:19265/1‑85 (MQ=37)
                                                                                                                     cAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgc                                                    <  8:19265/85‑1 (MQ=37)
                                                                                                                                                |                                                                                                           
CAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGAT  >  NC_000913/1708518‑1708769

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: