Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,843,231 C→A 18.2% R27R (CGG→AGG)  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,2310CA18.2% 73.0 / 5.2 36R27R (CGG→AGG) PP_4259iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (16/13);  new base A (2/5);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGATTGCGCGAGACTGCGGCCCGCACGCACAGGCCGCGGCGTTTCTCCACTGCCGTGGGACTTGATGATGAGCAAGCAAATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCCTGCTGTACTTCGGCACCGTGTGGCTGAACTGGGAGGGGCACCAGGCCGTGTGGTGGAACCTGC  >  NC_002947/4843084‑4843369
                                                                                                                                                    |                                                                                                                                          
gCAGATTGCGCGAGACTGCGGCCCGCACGCACAGGCCGCGGCGTTTCTCCACTGCCGTGGGACTTGATGATGAGCAAGCAAATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTc‑cc                                                                                                                                            <  6:233483/148‑1 (MQ=255)
 cAGATTGCGCGAGACTGCGGCCCGCACGCACAGGCCGCGGCGTTTCTCCACTGCCGTGGGACTTGATGATGAGCAAGCAAATTCCGGTACATGATGTGACCACGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTA‑CCgg                                                                                                                                          >  7:214934/1‑149 (MQ=255)
         gcgAGACTGCGGCCCGCACGCACAGGCCGCGGCGTTTCTCCACTGCCGTGGGACTTGATGATGAGCAAGCAAATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGAt                                                                                                                                  <  1:280639/149‑1 (MQ=255)
           gagaCTGCGGCCCGCACGCACAGGCCGCGGCGTTTCTCCACTGCCGTGGGACTTGATGATGAGCAAGCAAATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATc                                                                                                                                 >  6:128005/1‑148 (MQ=255)
             gaCTGCGGCCCGCACGCACAGGCCGCGGCGTTTCTCCACTGCCGTGGGACTTGATGATGAGCAAGCAAATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTac                                                                                                                              <  1:22531/149‑1 (MQ=255)
                tGCGGCCCGCACGCACAGGCCGCGGCGTTTCTCCACTGCCGTGGGACTTGATGATGAGCAAGCAAATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACAcc                                                                                                                           >  2:206389/1‑149 (MQ=255)
                      ccGCACGCACAGGCCGCGGCGTTTCTCCACTGCCGTGGGACTTGATGATGAGCAAGCAAATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGcc                                                                                                                     <  5:28630/149‑1 (MQ=255)
                          acgcacAGGCCGCGGCGTTTCTCCACTGCCGTGGGACTTGATGATGAGCAAGCAAATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CAGGGAAAAGATCTACACCCGCGCCTTCa                                                                                                                 >  2:133193/1‑149 (MQ=255)
                             cacaGGCCGCGGCGTATCTCCACTGCCGTGGGACGTGATGATGAGCAACCAAATTCCGGTACATGATGTGACCACGCCTGCCAACAAAGGGAAGGATTTTTTTTTTATTTATTCTTC‑GAGGTAATAGATCTACACCCGCGCCTTCACCg                                                                                                              <  8:214934/149‑1 (MQ=255)
                                aGGCCGCGGCGTTTCTCCACTGCCGTGGGACTTGATGATGAGCAAGCAAATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTCTCCGGGAAAAGATCTACACCCGCGCCTTCACCGGc                                                                                                            <  1:206389/149‑1 (MQ=255)
                                   ccGCGGCGTTTCTCCACTGCCGTGGGACTTGATGATGAGCAAGCAAATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGt                                                                                                        <  1:133193/149‑1 (MQ=255)
                                                ccACTGCCGTGGGACTTGATGATGAGCAAGCAAATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTCCGCTTC‑CCGGGAAAAGATCTACACCCCCACCTTCACCGGCCTGTTCCGTCGCCTgcg                                                                                           >  8:106596/1‑149 (MQ=255)
                                                      ccGTGGGACTTGATGATGAGCAAGCAAATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATgg                                                                                      >  8:24007/1‑148 (MQ=255)
                                                           ggACTTGATGATGAGCAAGCAAATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTcggcg                                                                                >  3:207788/1‑149 (MQ=255)
                                                                     atgaGCAAGCAAATTCCGGTACATGATGTCACCCCGCCTGCCAACAAAGGGAAGGATTTTTTTTTTATTTATTCTTC‑GAGGTAATAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGc                                                                       <  6:247207/148‑1 (MQ=255)
                                                                       gaGCAAGCAAATTCCGGTACATGATGTAACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCTTGCGCATGGTCGGCGGTGCAGTGCTGt                                                                    >  6:56366/1‑149 (MQ=255)
                                                                        agcaagcaAATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGt                                                                    <  7:24007/148‑1 (MQ=255)
                                                                          caagcaaATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CAGGGAAAAGATCTACACCCGCGCGTTCACCTGCCTGTTACGACGCGAGCGAAAGGTCGGCGGTCCAGTGCTGGTca                                                                 >  8:292490/1‑148 (MQ=255)
                                                                            agcaaATTCCGGTACATGATGTAACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCTTGCGCATGGTCGGCGGTGCAGTGCTGTTCctg                                                               <  5:56366/149‑1 (MQ=255)
                                                                              caaATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCctgct                                                             >  4:167893/1‑149 (MQ=255)
                                                                                        tACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTCTTCCTGCTGTACTTCgg                                                    >  6:301228/1‑148 (MQ=255)
                                                                                             gatgTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCCTGCTGTACTTCGGCACCgt                                              >  1:64008/1‑149 (MQ=255)
                                                                                              atgTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACAGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCCTGCTGTACTTCGGCAcc                                                >  3:41673/1‑146 (MQ=255)
                                                                                              atgTGACCCCGCCTGCCAACAAAGGGAAGGATTATGTTGATATTTATTCTTC‑GAGGTAATAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCCTGCTGTACTTCGGCAcc                                                <  4:41673/146‑1 (MQ=255)
                                                                                               tgtgACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCCTGCTGTACTTCGGCACCgtgt                                            <  1:153361/149‑1 (MQ=255)
                                                                                                            gCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGATGTTCCTGCTGTACTTCGGCACCGGGTGGCTGAACTgg                                 >  5:271234/1‑147 (MQ=255)
                                                                                                             ccAACAAAGGGAAGGATTCTGTTTATATTTATGCTTC‑CAGGTAATAGATCTACACCCGCGCCTTCACCGGCCTGTTCCg                                                                                                    <  2:52612/79‑1 (MQ=255)
                                                                                                             ccAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCg                                                                                                    >  1:52612/1‑79 (MQ=255)
                                                                                                               aacaaAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCCTGCTGTACTTCGGCACCGTGTGGCTGAACTGGGAggg                            >  6:112300/1‑149 (MQ=255)
                                                                                                                 caaAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCCTGCTGTACTTCGGCACCGTGTGGCTGAACTGGGAgggg                           >  5:76068/1‑148 (MQ=255)
                                                                                                                       gAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCCTGCTGTACTTCGGCACCGTGTGGCTGAACTGGGAGGGGCACCAgg                    <  3:167893/149‑1 (MQ=255)
                                                                                                                            aTTCCGTCGAACTCTACGCTTC‑GCGGGAATAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCCTGCTGTACTTCGGCACCGTGTGGCTGAACTGGGAGGGGCACCAGGCCgtg               <  4:207788/149‑1 (MQ=255)
                                                                                                                               ccGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCCTGCTGTACTTCGGCACCGTGTGGCTGAACTGGGAGGGGCACCAGGCCGtgtg             <  5:128005/148‑1 (MQ=255)
                                                                                                                                     aCCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCCTGCTGTACTTCGGCACCGTGTGGCTGAACTGGGAGGGGCACCAGGCCGTGTGGTGGAAc      <  6:76068/149‑1 (MQ=255)
                                                                                                                                            cGCTTC‑CAGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCCTGCTGTACTTCGGCACCGTGTGGCTGAACTGGGAGGGGCACCAGGCCGTGTGGTGGAACCTGc  <  2:64008/146‑1 (MQ=255)
                                                                                                                                              cTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCctgct                                                             <  7:295907/85‑1 (MQ=255)
                                                                                                                                              cTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCctgct                                                             >  8:295907/1‑85 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
GCAGATTGCGCGAGACTGCGGCCCGCACGCACAGGCCGCGGCGTTTCTCCACTGCCGTGGGACTTGATGATGAGCAAGCAAATTCCGGTACATGATGTGACCCCGCCTGCCAACAAAGGGAAGGATTCCGTCGACCTCTACGCTTC‑CCGGGAAAAGATCTACACCCGCGCCTTCACCGGCCTGTTCCGTCGCCTGCGCATGGTCGGCGGTGCAGTGCTGTTCCTGCTGTACTTCGGCACCGTGTGGCTGAACTGGGAGGGGCACCAGGCCGTGTGGTGGAACCTGC  >  NC_002947/4843084‑4843369

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: ‑