Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,708,662 C→G 16.7% 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.7% 97.1 / 9.7 46Q582E (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 (5/32);  new base G (5/3);  total (10/36)
Fisher's exact test for biased strand distribution p-value = 8.25e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCG  >  NC_000913/1708513‑1708782
                                                                                                                                                     |                                                                                                                        
aaGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGaac                                                                                                                          <  1:754356/150‑1 (MQ=255)
 aGCGCAAGCCAGAGCGAAACAGGCAGATCTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGaaca                                                                                                                         <  1:278110/150‑1 (MQ=255)
  gcgcAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGaacaa                                                                                                                        >  2:234480/1‑150 (MQ=255)
     cAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAg                                                                                                                     >  2:279044/1‑150 (MQ=255)
      aaGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAgg                                                                                                                    <  1:908450/150‑1 (MQ=255)
       aGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGt                                                                                                                   <  1:281280/150‑1 (MQ=255)
       aGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTAGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGt                                                                                                                   <  1:39322/150‑1 (MQ=255)
        gCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCCGCTATCGCCCGCGCCAAAGCGCGCACGCAGGAACAGAAGAATGCTCAAGCGGCACCAGCACAGCACGta                                                                                                                  >  1:670096/1‑149 (MQ=255)
        gCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTc                                                                                                                  <  1:177828/150‑1 (MQ=255)
        gCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTc                                                                                                                  <  2:572836/150‑1 (MQ=255)
            tagCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAATCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATc                                                                                                              <  1:612667/149‑1 (MQ=255)
            gagCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATc                                                                                                              <  1:1241699/150‑1 (MQ=255)
            gagCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATc                                                                                                              <  2:805063/150‑1 (MQ=255)
             agCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGCGGAACCAGAAGAACAggacaaagc                                                                                                             >  2:1259977/1‑143 (MQ=255)
                gAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCcgcg                                                                                                          <  1:752525/150‑1 (MQ=255)
                 aaaCAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCgcgc                                                                                                         <  1:430296/150‑1 (MQ=255)
                  aaCAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCa                                                                                                        <  2:509499/150‑1 (MQ=255)
                     aGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAg                                                                                                     <  2:165927/150‑1 (MQ=255)
                      ggCAGATCTGCAGCATACTCAGGACCCAGAAAACGTTGCTGATCCTCGTAAAACTGCCGTTGAAGCAGCTATTGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAgc                                                                                                    <  2:208937/150‑1 (MQ=14)
                       gCAGATCTTCAGCAAACTTCCGACGCAGCAACAGATGCTGATCCACGTAAATCTGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAgcg                                                                                                   <  1:759336/150‑2 (MQ=14)
                          gAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAACGCGCGCAAGCTGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTGGATCCGCGCAAGgcagcg                                                                                                >  1:1227671/1‑149 (MQ=25)
                           aaCTGCAGCAAACTAACGACGCCGCAACCGTTGCTGATCCTCGTAAACCAGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAACCGGGAACCGCAACCGGCTAAGGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGctgccg                                                                                               <  1:234480/150‑1 (MQ=17)
                               gcagcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCAGtagc                                                                                           >  2:56229/1‑147 (MQ=255)
                                     aaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                                     <  1:238581/150‑1 (MQ=255)
                                     aaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGc                                                                                     <  1:997362/150‑1 (MQ=255)
                                     aaCTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTAGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACCGCAACAGGCTAATGCGGACCCAGAAGAACAGGGAGATCCGCGCAAAGCCGCCGTAGCAGCCGc                                                                                     >  1:60705/1‑150 (MQ=18)
                                          acgacgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTg                                                                                <  2:692141/150‑1 (MQ=255)
                                           cgacgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAATGCGGAACCAGAAGAACAGGTAGATACGCGCAAACCACCCGTCGAAGCAGCTAttta                                                                               >  2:1188501/1‑148 (MQ=18)
                                             acgCAGCTACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCTACAGGCTTATGCGGATCCAGAACAAGAGGTCGATCCGCGCAATGCCGCCTTCGAAGCCGCTTTTGccc                                                                             <  1:897307/150‑1 (MQ=18)
                                              cgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCTACAGGCTAATGCGGAACCAGAACATCAGGTCGATCCGCGCAAATCCGCCGTCGAAGCCGCTATTGCCCg                                                                            <  2:932823/150‑1 (MQ=255)
                                               gcagcaACCGTTGCTGATCCACGTAAATCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                           <  2:304535/150‑1 (MQ=255)
                                               gcagcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGt                                                                           <  2:1077873/150‑1 (MQ=255)
                                                  gcaACCGTTGCTGATCCACGTAAAACTGCGGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGcc                                                                        <  2:1137454/150‑1 (MQ=255)
                                                   caACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCa                                                                       <  2:1178680/150‑1 (MQ=255)
                                                            cTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGACCCAGAACAACAGGTGGATCCGCGCAAAGCCGCCGTGGCAGCGGCTATTGCCCGCGCGAAGgcccgaa                                                              >  2:887965/1‑148 (MQ=17)
                                                              taTCCACGTAAAACTGCCGTTGATGCAGCAATCGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGACCACCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAg                                                            <  2:231702/149‑1 (MQ=17)
                                                               aTCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAATCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGc                                                           <  2:743329/150‑1 (MQ=255)
                                                                tCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTCGGGAGGGCCCAGGACAACAGGACGATCCGCGCCAAGCCGCCGTCGCAGGCGCTATTGGCCGTGCCAAAGCGCGCACGcg                                                          >  1:379172/1‑149 (MQ=17)
                                                                         aaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCa                                                 <  1:533651/150‑1 (MQ=255)
                                                                           cTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAAc                                               <  1:1111979/150‑1 (MQ=255)
                                                                            tGCCGGCGAAGCAGCTAGCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTACTGCGGAACCAGAACAACAGCGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACa                                              <  1:686178/150‑1 (MQ=14)
                                                                                  tGAAGCAGCTATCGCCCGCGCCAACGCGCGCAAGCTGGAACAGCAACAGGCAAATCCGGAACCAGAACACCAGGTCGATCCCCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTaa                                        <  2:280447/150‑1 (MQ=14)
                                                                                     agcagcTATCTCACGGGACTAAGCGGGCAAGCTGGAAGAGCAACAGGCTAATTCGGAACCAGAACACCAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGc                                     <  2:239346/150‑1 (MQ=11)
                                                                                      gcagcTATCTCCCGCGCAAAAGCGCGCAAGCTGGAACATCAACAGGATACTGCGGAACCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGGGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCg                                    <  2:1107452/150‑1 (MQ=14)
                                                                                                     gcCAAAGCGCGGAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGCACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAg                     <  2:600394/150‑1 (MQ=255)
                                                                                                      ccAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgc                                                                   >  2:510601/1‑103 (MQ=17)
                                                                                                      ccAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgc                                                                   <  1:510601/103‑1 (MQ=17)
                                                                                                                cAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCcgcg          <  1:609879/150‑1 (MQ=255)
                                                                                                                        aacagcaacagcAAAATGAGGAAGCAGAAACCGTTGTCGATGCCCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGAGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgccg  <  2:625596/150‑1 (MQ=11)
                                                                                                                                                     |                                                                                                                        
AAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCG  >  NC_000913/1708513‑1708782

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: