Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,505,094 C→A 8.1% L52L (CTC→CTA tusE → sulfur transfer protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,505,0940CA 8.1% 150.4 / 2.7 50L52L (CTC→CTAtusEsulfur transfer protein
Reads supporting (aligned to +/- strand):  ref base C (19/26);  new base A (2/2);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

CGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAGGGCAACAGCCCGCACCTTAACCGCCTGTTCAACGGCACCCCCGCCAAACTTGCC  >  NC_002947/4504946‑4505232
                                                                                                                                                    |                                                                                                                                          
cGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGCTGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTcc                                                                                                                                           <  1:193815/150‑1 (MQ=255)
  cTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGt                                                                                                                                         <  1:1000141/150‑1 (MQ=255)
    cACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCa                                                                                                                                       <  2:521578/150‑1 (MQ=255)
         agggCGATCACCTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCt                                                                                                                                  >  2:677918/3‑150 (MQ=255)
                 cAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGCCCACTGGGAAATCCTCGAACTGCGCCGGCAGTTCTACCAGGaa                                                                                                                          <  2:182714/150‑1 (MQ=255)
                   ggTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATa                                                                                                                        >  1:1122504/1‑150 (MQ=255)
                         cGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCt                                                                                                                  <  1:779737/150‑1 (MQ=255)
                                     ggATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCAc                                                                                                      <  2:280552/150‑1 (MQ=255)
                                        tGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCAcccg                                                                                                   >  1:858253/1‑150 (MQ=255)
                                            ttCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCAcccgcccg                                                                                               >  2:847177/1‑150 (MQ=255)
                                            ttCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGCCAGtt                                                                                                                                    >  2:709119/1‑113 (MQ=255)
                                            ttCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGCCAGtt                                                                                                                                    <  1:709119/113‑1 (MQ=255)
                                             tCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACccgcccgc                                                                                              >  2:181972/1‑150 (MQ=255)
                                              ccTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCt                                                                                             <  1:677918/150‑1 (MQ=255)
                                                          cgcatcttGATCACACGACGCCGCCGAGGCTGTGGCCAAGGGTGAGGGTATCCCGTTGACCGCGGACCACTGGGAATGCCTCGAACTGCTCCGTCAGTATTACCAGGAATACCAGCTGTCGCCGTCCTCCCGCCCGCTGATCAAGTacac                                                                                 <  2:410701/144‑1 (MQ=255)
                                                                ttGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAg                                                                                                                             >  1:703408/1‑100 (MQ=255)
                                                                ttGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAg                                                                                                                             <  2:703408/100‑1 (MQ=255)
                                                                  ggTCTCACGACGCCGCCGATGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCGCCGGCCACCCGCCCGCTGATCAATTACACGGCCCTCa                                                                         <  2:296322/150‑1 (MQ=255)
                                                                   gTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCaa                                                                        >  2:153995/1‑150 (MQ=255)
                                                                    tctcACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAg                                                                       >  2:338139/1‑150 (MQ=255)
                                                                       cACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTTTTGGGCGGAGATGGTGCGGGTGTGAAATCCGCGAAATGCGCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTg                                                                    <  1:805629/150‑1 (MQ=255)
                                                                       cACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCGCTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTg                                                                    >  2:598851/1‑150 (MQ=255)
                                                                       cACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTg                                                                    <  2:12709/150‑1 (MQ=255)
                                                                        acgacgCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTgt                                                                   >  1:462373/1‑149 (MQ=255)
                                                                          gacgCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAg                                                                                                                    >  2:547673/1‑99 (MQ=255)
                                                                          gacgCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAg                                                                                                                    <  1:547673/99‑1 (MQ=255)
                                                                           acgCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATTCCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGcc                                                                <  2:427898/150‑1 (MQ=255)
                                                                           acgCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGcc                                                                <  1:338139/150‑1 (MQ=255)
                                                                                  cgAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGaaaa                                                         <  1:686597/150‑1 (MQ=255)
                                                                                    aGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAgg                                                       <  1:989065/150‑1 (MQ=255)
                                                                                      gCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAGGGc                                                     <  2:1010180/150‑1 (MQ=255)
                                                                                        tttGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAGGGCaa                                                   >  1:1089928/1‑150 (MQ=255)
                                                                                            gCCAAGCGCGAGGGTATCCCGTTGACGGCGGAGCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAGGGCAACAGc                                               >  2:626067/1‑150 (MQ=255)
                                                                                             ccAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAGGGCAACAGcc                                              <  2:441424/150‑1 (MQ=255)
                                                                                              cAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAGGGCAACAGccc                                             <  1:616223/150‑1 (MQ=255)
                                                                                              cAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAGGGCAACAGccc                                             <  1:181972/150‑1 (MQ=255)
                                                                                               aaGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAGGGCAACAGCCCg                                            <  1:377088/150‑1 (MQ=255)
                                                                                                             cccGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAGGGCAACAGCCCGCACCTTAACCGCCt                              <  2:858253/150‑1 (MQ=255)
                                                                                                                   gACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGaa                                                                                                                          <  2:969997/52‑1 (MQ=255)
                                                                                                                   gACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGaa                                                                                                                          <  2:399671/52‑1 (MQ=255)
                                                                                                                   gACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGaa                                                                                                                          >  1:969997/1‑52 (MQ=255)
                                                                                                                   gACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGaa                                                                                                                          >  1:399671/1‑52 (MQ=255)
                                                                                                                                    aaaTCTTCGAACTGCTTCGTCAGTGCTACCAGGAATACCAGCTGTCTCCGGCCACGCACCCGCTGATCAAGTACACGGCCCTCAAGCTACACCCTGCAAAGGGCCACAGCCCGCACCTTAACTGCCTGTTCAACATCCCCCTCGCCacaa       >  2:751785/1‑147 (MQ=255)
                                                                                                                                        ccTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAGGGCAACAGCCCGCACCTTAACCGCCTGTTCAACGGCACCCCCGCCAAACTTgc   <  1:598851/150‑1 (MQ=255)
                                                                                                                                        ccTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCATCAAGCTGGGCCCTGAAAAGGGCAACAGCCCGCACCTTAACCGCCTGTTCAACGGCACCCCCGCCAAACTTgc   >  1:826852/1‑150 (MQ=255)
                                                                                                                                         cTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAGGGCAACAGCCCGCACCTTAACCGCCTGTTCAACGGCACCCCCGCCAAACTTgcc  >  2:1037322/1‑150 (MQ=255)
                                                                                                                                         cTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGAAAATAAAGCTGGGCCCTGAAAAGGGCAACAGCCCGCACCTTAACCGCCTGTTCAACGGCACCCCCGCCAAACTTgcc  >  2:821033/1‑150 (MQ=255)
                                                                                                                                         cTCGAACTGCTACGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAggg                                                      <  2:481001/98‑1 (MQ=255)
                                                                                                                                         cTCGAACTGCTACGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAggg                                                      <  2:480844/98‑1 (MQ=255)
                                                                                                                                         cTCGAACTGCTACGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAggg                                                      >  1:481001/1‑98 (MQ=255)
                                                                                                                                         cTCGAACTGCTACGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAggg                                                      >  1:480844/1‑98 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
CGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATACCAGCTGTCTCCGGCCACCCGCCCGCTGATCAAGTACACGGCCCTCAAGCTGGGCCCTGAAAAGGGCAACAGCCCGCACCTTAACCGCCTGTTCAACGGCACCCCCGCCAAACTTGCC  >  NC_002947/4504946‑4505232

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