Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,573,801 C→T 15.3% G278G (GGC→GGT rsgA → ribosome biogenesis GTPase RsgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,573,8010CT15.3% 89.0 / 8.1 39G278G (GGC→GGTrsgAribosome biogenesis GTPase RsgA
Reads supporting (aligned to +/- strand):  ref base C (18/15);  new base T (3/3);  total (21/18)
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.79e-01

AGCCTGCTGCCAGACGCTGGCACCCGCGTCGGCGACTTGTCGGAATGGTCAGGCCAGGGTACCCACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAGCATGACCGCGAACCGGGCTGCGCGCTGCTCAAGGCACTGGAAGAGGGGCGCATCAAGCCGCAG  >  NC_002947/5573652‑5573945
                                                                                                                                                      |                                                                                                                                                
aGCCTGCTGCCAGACGCTGGCACCCGCGTCGGCGACTTGTCGGAATGGTCAGGCCAGGGTACCCACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGc                                                                                                                                                  <  1:675530/150‑1 (MQ=255)
aGCCTGCTGCCAGACGCTGGCACCCGCGTCGGCGACTTGTCGGAATGGTCAGGCCAGGGTACCCACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGc                                                                                                                                                  <  1:675451/150‑1 (MQ=255)
   ctgctgCCAGACGCTGGCACCCGCGTCGGCGACTTGTCGGACTGGCGAGGGGAGGATTCCCACACCACCACAACGCCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCAt                                                                                                                                               <  1:406065/150‑1 (MQ=255)
     gctgcCAGACGCTGGCACCCGCGTCGGCGACTTGTCGGAATGGTCAGGCCAGGGTACCCACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGt                                                                                                                                             <  2:541544/150‑1 (MQ=255)
                 tGGCACCCGCGTCGGCGACTTGTCGGAATGGTCAGGCCAGGGTACCCACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGcgacga                                                                                                                                 <  2:28657/150‑1 (MQ=255)
                       ccGCGTCGGCGACTTGTCGGAATGGTCAGGCCAGGGTACCCACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTgga                                                                                                                           >  1:403514/1‑150 (MQ=255)
                           gTCGGCGACTTGGCGGAAGGGTCAGGCCAGGGCACCCAACCCACCACCACGGCGCGGCTTTACCACTGCCCCAATGGCGGGGACCTGATAGACTCGCCGGGGATTCGCGAGCTCGGCCTTGGCCATGTCAGCCGTGACGACGTGGAGGAt                                                                                                                       <  2:95582/150‑1 (MQ=255)
                               gCGACTTGTCGGAATGGTCAGGCCAGGGTACCCACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCGATGTGAGCCGCGACGGCGTGGAGGATGGCt                                                                                                                   >  2:624077/1‑150 (MQ=255)
                                                  aGGCCAGGGTACCCACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGTAGGATGGCTTCATCGAGTTCCGCGACCt                                                                                                >  1:525492/1‑150 (MQ=255)
                                                      cAGGGTACCCACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTc                                                                                            >  1:591704/1‑150 (MQ=255)
                                                       atggTACCCACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCg                                                                                           <  1:563783/148‑1 (MQ=255)
                                                        gggTACCCACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCgg                                                                                          >  1:302882/1‑150 (MQ=255)
                                                              ccACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCt                                                                                                                   <  2:641131/119‑1 (MQ=255)
                                                              ccACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCt                                                                                                                   >  1:641131/1‑119 (MQ=255)
                                                              ccACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTc                                                                                                                 >  1:258937/1‑121 (MQ=255)
                                                              ccACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTc                                                                                                                 <  2:258937/121‑1 (MQ=255)
                                                                acacCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCAGATAACTGCACCTGGTATTCCGGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGcc                                                                                  >  1:531943/1‑150 (MQ=255)
                                                                  accaccaccacGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGc                                                                                >  1:310101/1‑150 (MQ=255)
                                                                              gcgcGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGTCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTc                                                                                            >  2:135905/1‑126 (MQ=255)
                                                                              gcgcGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGTCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTc                                                                                            >  2:135998/1‑126 (MQ=255)
                                                                              gcgcGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGTCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTc                                                                                            <  1:135998/126‑1 (MQ=255)
                                                                              gcgcGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGTCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTc                                                                                            <  1:135905/126‑1 (MQ=255)
                                                                              gcgcGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGTCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGCTCCGCGACCTGTTc                                                                                            <  1:135908/126‑1 (MQ=255)
                                                                              gcgcGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGTCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGCTCCGCGACCTGTTc                                                                                            >  2:135908/1‑126 (MQ=255)
                                                                                 cgGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAg                                                                 <  2:302882/150‑1 (MQ=255)
                                                                                    cTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAGCAt                                                              <  2:403514/150‑1 (MQ=255)
                                                                                                 ccAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAGCATGACCGCGAACCgg                                                 <  1:624077/150‑1 (MQ=255)
                                                                                                          gggACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGTAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGATTGCAAGCATGACCGCGAACCGGGCTgcgcgc                                        <  2:525492/150‑1 (MQ=255)
                                                                                                                 gATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAGCATGACCGCGAACCggg                                                >  2:354798/1‑135 (MQ=255)
                                                                                                                 gATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAGCATGACCGCGAACCggg                                                <  1:354798/135‑1 (MQ=255)
                                                                                                                              gggTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGAt                                                                                                                       >  2:668936/1‑52 (MQ=255)
                                                                                                                              gggTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGAt                                                                                                                       <  1:668936/52‑1 (MQ=255)
                                                                                                                              gggTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAGCATGACCGCGAACCGGGCTGCGCGCTGCTCAAGGCACTGGAaga                     >  2:291798/1‑150 (MQ=255)
                                                                                                                                   ttCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAGCATGACCGCGAACCGGGCTGCGCGCTGCTCAAGGCACTGGAAGAGGGgc                >  2:560140/1‑150 (MQ=255)
                                                                                                                                   ttCGCGAGTTCGGCCTTCGCAATGTCCGGCGCCACTACCTGGAAGATGCCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAGCATGACCGCGAACCGGGCTGCGCGCTGCTCAAGGCACTGGAAGAGGGgc                <  2:310101/150‑1 (MQ=255)
                                                                                                                                      gcgAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAGCATGACCGCGAACCGGGCTGCGCGCTGCTCAAGGCACTGGAAGAGGGGCGCa             >  1:237484/1‑150 (MQ=255)
                                                                                                                                       cgAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAGCATGACCGCGAACCGGGCTGCGCGCTGCTCAAGGCAc                            <  2:757669/134‑1 (MQ=255)
                                                                                                                                       cgAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAGCATGACCGCGAACCGGGCTGCGCGCTGCTCAAGGCAc                            >  1:757669/1‑134 (MQ=255)
                                                                                                                                        gAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAGCATGACCGCGAACCGGGCTGCGCGCTGCTCAAGGCACTGGAAGAGGGGCGCATc           >  1:71767/1‑150 (MQ=255)
                                                                                                                                              ggCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAGCATGACCGCGAACCGGGCTGCGCGCTGCTCAAGGCACTGGAAGAGGGGCGCATCAAGCCg     >  2:71700/1‑150 (MQ=255)
                                                                                                                                                 cTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAGCATGACCGCGAACCGGGCTGCGCGCTGCTCAAGGCACTGGAAGAGGGGCGCATCAAGCCGCAg  >  1:487355/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
AGCCTGCTGCCAGACGCTGGCACCCGCGTCGGCGACTTGTCGGAATGGTCAGGCCAGGGTACCCACACCACCACCACGGCGCGGCTTTACCACTTCCCCAATGGCGGGGACCTGATCGACTCGCCGGGTATTCGCGAGTTCGGCCTTGGCCATGTCAGCCGCGACGACGTGGAGGATGGCTTCATCGAGTTCCGCGACCTGTTCGGCACCTGCCGCTTCCGCGACTGCAAGCATGACCGCGAACCGGGCTGCGCGCTGCTCAAGGCACTGGAAGAGGGGCGCATCAAGCCGCAG  >  NC_002947/5573652‑5573945

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