Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 258,182 A→G 8.8% K27K (AAA→AAG)  PP_0206 → 4Fe‑4S ferredoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947258,1820AG 8.8% 126.0 / 3.2 46K27K (AAA→AAG) PP_02064Fe‑4S ferredoxin
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base G (2/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

TGGGGCCATCAGCCTCTGCAAAAAACAGAACGAGGGCCACCGGCCCCAAGGACTCCGTGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTACTGCATGCCATGCGAGAAGGACTGCCCCACCGGTGCGGTGAAGGTCGAC  >  NC_002947/258041‑258326
                                                                                                                                             |                                                                                                                                                
tGGGGCCATCAGCCTCTGCAAAAAACAGAACGAGGGCCACCGGCCCCAAGGACTCCGTGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCAc                                                                                                                                          >  2:1100681/1‑150 (MQ=255)
   ggCCATCAGCCTCTGCAAAAAACAGAACGAGGGCCACCGGCCCCAAGGACTCCGTGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCgt                                                                                                                                       >  2:697750/1‑150 (MQ=255)
                   aaaaaaCAGAACGAGGGCCACCGGCCCCAAGGACTCCGTGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAGGGCTGCACCGTGTGCGTCGAGGTCTGc                                                                                                                       >  2:1268308/1‑150 (MQ=255)
                          aGAACGAGGGCCACCGGCCCCAAGGACTCCGTGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCGGAAAAAGGCTGCACCGTGTGCGTCGAGGTCt                                                                                                                         >  1:461416/1‑141 (MQ=255)
                          aGAACGAGGGCCACCGGCCCCAAGGACTCCGTGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCt                                                                                                                         <  2:461416/141‑1 (MQ=255)
                           gAACGAGGGCCACCGGCCCCAAGGACTCCGTGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGa                                                                                                               >  2:1127894/1‑150 (MQ=255)
                               gAGGGCCACCGGCCCCAAGGACTCCGTGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACctg                                                                                                           >  2:147385/1‑150 (MQ=255)
                                  ggCCACCGGCCCCAAGGACTCCGTGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGCCctgctg                                                                                                        <  1:54011/150‑1 (MQ=255)
                                     cACCGGCCCCAAGGACTCCGTGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAGGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCg                                                                                                     <  1:1268308/150‑1 (MQ=255)
                                      aCCGGCCTCAAGGACTCCGTGAAATTGCCTACCAGGCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAACGCTGCACCGTGTGCGTCGAGGTCTGGCCCATGGACCTGCTGGCGa                                                                                                    <  2:191834/150‑1 (MQ=255)
                                               aaGGACTCCGTGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTgg                                                                                                       <  2:384621/138‑1 (MQ=255)
                                               aaGGACTCCGTGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTgg                                                                                                       >  1:384621/1‑138 (MQ=255)
                                                        gTGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGc                                                                                                      >  1:94644/1‑130 (MQ=255)
                                                        gTGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGc                                                                                                      <  2:94644/130‑1 (MQ=255)
                                                         tGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCaaaa                                                                                 <  2:313250/150‑1 (MQ=255)
                                                          gAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAg                                                                                <  1:9759/150‑1 (MQ=255)
                                                                     ccAGCCGCAAGAAATCTTCTTCCACAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGc                                                                     >  2:740052/1‑150 (MQ=255)
                                                                            cAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAAccc                                                                                             >  2:774841/1‑119 (MQ=255)
                                                                            cAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAAccc                                                                                             <  1:774841/119‑1 (MQ=255)
                                                                            cAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTcgac                                                              >  1:802577/1‑150 (MQ=255)
                                                                             aaGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCaa                                                                                                <  2:106750/115‑1 (MQ=255)
                                                                             aaGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCaa                                                                                                >  1:106750/1‑115 (MQ=255)
                                                                               gAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGaa                                                           >  1:535150/1‑150 (MQ=255)
                                                                                   tcttctTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACAt                                                                        >  1:490386/1‑133 (MQ=255)
                                                                                   tcttctTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACAt                                                                        <  2:490386/133‑1 (MQ=255)
                                                                                    cttcttCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTg                                                      <  1:697750/150‑1 (MQ=255)
                                                                                        ttcCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCAGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTAc                                                  >  2:237654/1‑150 (MQ=255)
                                                                                            gcagcaGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTACTGca                                              <  2:409671/150‑1 (MQ=255)
                                                                                                   cgcCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTACTGcatgccatg                                       >  1:466158/1‑150 (MQ=255)
                                                                                                    gcCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTACTGcatgccatgc                                      <  1:1100681/150‑1 (MQ=255)
                                                                                                      cccGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTACTGCATGCCATGCga                                    >  1:401368/1‑150 (MQ=255)
                                                                                                      cccGGTCACCATCGACGGAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCGGCTGGCGATCAACCCGGCTACGAAAACAGCCTTCATAGCTTTCGACGAATGCTGGTACTGCATGCCATGCga                                    >  1:1048829/1‑150 (MQ=255)
                                                                                                         ggTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTACTGCATGCCATGCGAGaa                                 <  2:986397/150‑1 (MQ=255)
                                                                                                          gTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTACTGCATGCCATGCGAGAAg                                <  1:1127894/150‑1 (MQ=255)
                                                                                                              ccATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTACTGCATGCCATGCGAGAAGGACt                            <  2:1186545/150‑1 (MQ=255)
                                                                                                               cATCGACGTAGACAAATGCATCGCCGAAAACGGCTGCACAGTGTCTGCTGCGGTCTGCCCCATGGACCTGCGGGCGATCAACCCGGCTACGCAAAAGGCCTACATGACTTTCGACGAAAGCTGATACTGCATGCCATGCGAGAAGGACTg                           >  1:1186545/1‑150 (MQ=255)
                                                                                                                aTCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTACTGCATGCCATGCGAGAAGGACTGc                          <  2:401368/150‑1 (MQ=255)
                                                                                                                    aCGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTACTGCATGCCATGCGAGAAGGACTGCCCCa                      >  1:812797/1‑150 (MQ=255)
                                                                                                                      gTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCaaaa                                                                                 <  2:923979/89‑1 (MQ=255)
                                                                                                                      gTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCaaaa                                                                                 >  1:923979/1‑89 (MQ=255)
                                                                                                                          aCAAATGCATCGCCGAAAAGGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCaa                                                                                                <  1:1108990/70‑1 (MQ=255)
                                                                                                                          aCAAATGCATCGCCGAAAAGGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCaa                                                                                                >  2:1108990/1‑70 (MQ=255)
                                                                                                                          aCAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTACTGCATGCCATGCGAGAAGGGCTGCCCCACcggtg                <  1:1204302/150‑1 (MQ=255)
                                                                                                                                  aTCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTACTGCATGCCATGCGAGAAGGACTGCCCCACCGGTGCGGTGAAg        <  1:147385/150‑1 (MQ=255)
                                                                                                                                  aTCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTACTGCATGCCATGCGAGAAGGACTGCCCCACCGGTGCGGTGAAg        >  1:454032/1‑150 (MQ=255)
                                                                                                                                        gAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTACTGCATGCCATGCGAGAAGGACTGCCCCACCGGTGCGGTGAAGGTCGAc  >  2:1188081/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                
TGGGGCCATCAGCCTCTGCAAAAAACAGAACGAGGGCCACCGGCCCCAAGGACTCCGTGAAATGGCCTACCAGCCGCAAGAAATCTTCTTCCGCAGCAGCGCCCCGGTCACCATCGACGTAGACAAATGCATCGCCGAAAAAGGCTGCACCGTGTGCGTCGAGGTCTGCCCCATGGACCTGCTGGCGATCAACCCGGCTACGCAAAAGGCCTACATGGCTTTCGACGAATGCTGGTACTGCATGCCATGCGAGAAGGACTGCCCCACCGGTGCGGTGAAGGTCGAC  >  NC_002947/258041‑258326

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑