Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,843,875 G→C 8.4% L241L (CTG→CTC)  PP_4259 → iron‑sulfur cluster‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,843,8750GC 8.4% 141.5 / 2.7 46L241L (CTG→CTC) PP_4259iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (22/20);  new base C (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 = 9.99e-01

GCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCTGC  >  NC_002947/4843735‑4844023
                                                                                                                                            |                                                                                                                                                      
gCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTc                                                                                                                                                <  4:120585/149‑1 (MQ=255)
   gATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGATCAGCGTGTCCtt                                                                                                                                             <  6:104401/149‑2 (MQ=255)
       gCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGa                                                                                                                                          <  6:268743/148‑1 (MQ=255)
                       ggTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCt                                                                                                                                              <  3:128937/128‑1 (MQ=255)
                       ggTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCt                                                                                                                                              >  4:128937/1‑128 (MQ=255)
                            gCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCcgc                                                                                                                    <  7:23016/149‑1 (MQ=255)
                                       cTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCgcgc                                                                                                           >  2:212637/1‑147 (MQ=255)
                                         cGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaag                                                                                                        >  2:247873/1‑148 (MQ=255)
                                           ccACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACACTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCaagaaa                                                                                                     <  3:209994/149‑1 (MQ=255)
                                               ttACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCcgcg                                                                                                            >  1:46772/1‑138 (MQ=255)
                                               ttACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCcgcg                                                                                                            <  2:46772/138‑1 (MQ=255)
                                                 aCGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATAc                                                                                                <  3:85947/148‑1 (MQ=255)
                                                                 gctACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAgg                                                                             <  4:5630/149‑1 (MQ=255)
                                                                   tACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGc                                                                           >  2:163394/1‑149 (MQ=255)
                                                                    aCGCGAACAGGTGTGCGTGTACATCTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGTCGAAACCCGCGGACCGCGCAAGAAAGATAACGACTACAAGG‑‑CCCAGGGcc                                                                          >  7:32699/1‑149 (MQ=255)
                                                                     cgcgAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCt                                                                         >  1:236284/1‑149 (MQ=255)
                                                                        gAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTggg                                                                      <  2:236284/149‑1 (MQ=255)
                                                                            aGGTGTGCGTGTACATGTGCCCGTATGCCCGCTTCCTGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGAc                                                                  <  6:139974/149‑1 (MQ=255)
                                                                              gtgtGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTg                                                                >  3:236742/1‑149 (MQ=255)
                                                                                gtgCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCa                                                              >  3:17028/1‑149 (MQ=255)
                                                                                       tACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAgg                                                                             >  6:236438/1‑127 (MQ=255)
                                                                                       tACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAgg                                                                             <  5:236438/127‑1 (MQ=255)
                                                                                        aCATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGAGAAACCCGCGGCCCGCGCAAGAAAGATACAGACTACAAGG‑‑CCCAGGGCCTGCGAGACTCCATCGACTg                                                       >  5:254102/1‑148 (MQ=37)
                                                                                         cATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGCCTGCATCGACTGCa                                                     >  6:133929/1‑149 (MQ=255)
                                                                                               ccccTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCAtg                                                <  6:74852/148‑1 (MQ=255)
                                                                                                 ccTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATgtgc                                             <  1:163394/149‑1 (MQ=255)
                                                                                                     tGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATgtgcgt                                           <  5:133929/147‑1 (MQ=255)
                                                                                                       cccGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGa                                                          >  1:187292/1‑130 (MQ=255)
                                                                                                       cccGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGa                                                          <  2:187292/130‑1 (MQ=255)
                                                                                                          gCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTc                                    >  7:91026/1‑149 (MQ=255)
                                                                                                                agagCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCg                              >  3:13358/1‑149 (MQ=255)
                                                                                                                       gATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCa                       >  7:102242/1‑149 (MQ=255)
                                                                                                                              gACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGAc                                                                                            >  2:247530/1‑75 (MQ=255)
                                                                                                                              gACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGAc                                                                                            <  1:247530/75‑1 (MQ=255)
                                                                                                                              gACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGGCCCCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGa                   <  1:212637/148‑1 (MQ=255)
                                                                                                                              gACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACAt                >  7:54489/1‑149 (MQ=255)
                                                                                                                                cAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATcc              <  1:247873/149‑1 (MQ=255)
                                                                                                                                 taGGACACGATCAGCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCg             >  8:245530/2‑149 (MQ=255)
                                                                                                                                 aaGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCg             >  7:68194/1‑149 (MQ=255)
                                                                                                                                  aGGCCACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGt            <  8:214142/149‑1 (MQ=255)
                                                                                                                                   ggACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGt            >  8:190861/1‑148 (MQ=255)
                                                                                                                                    gACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGa          <  8:102242/149‑1 (MQ=255)
                                                                                                                                    gACACGATCATGGTGTCCTACGACCCACGCCGTGGCGAAACCCTCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCGTTGGTGACTGCATCGACTGCACCATgtgc                                             <  8:6893/114‑1 (MQ=255)
                                                                                                                                    gACACGATCAGCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATgtgc                                             >  7:6893/1‑114 (MQ=255)
                                                                                                                                         gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGcc     >  7:68675/1‑149 (MQ=255)
                                                                                                                                         gCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGcc     >  3:114076/1‑149 (MQ=255)
                                                                                                                                            gATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCTGc  <  8:68194/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                      
GCCGATGGCTGGGCCTACTTCTGGGTCGGCTTCTTCACCCTCGCCACTTACGGCAACGCCGGCTGGCTACGCGAACAGGTGTGCGTGTACATGTGCCCCTATGCCCGCTTCCAGAGCGTGATGTTCGACAAGGACACGCTGATCGTGTCCTACGACCCACGCCGTGGCGAAACCCGCGGCCCGCGCAAGAAAGATACCGACTACAAGG‑‑CCCAGGGCCTGGGTGACTGCATCGACTGCACCATGTGCGTGCAGGTCTGCCCGACCGGCATCGACATCCGTGACGGCCTGC  >  NC_002947/4843735‑4844023

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