Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,472,797 C→T 7.2% G8S (GGC→AGC)  rsfS ← ribosomal silencing factor RsfS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,472,7970CT 7.2% 158.3 / 2.4 55G8S (GGC→AGC) rsfSribosomal silencing factor RsfS
Reads supporting (aligned to +/- strand):  ref base C (28/23);  new base T (2/2);  total (30/25)
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.94e-01

GCGTTGATCTGGCGGTTGGAGGTACCGGTGGCAATGATCATGTAGTCGGTCAGGCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCG  >  NC_002947/5472655‑5472943
                                                                                                                                               |                                                                                                                                                   
gcgTTGATCTGGCGGTTGGAGGTACCGGTGGCAATGATCATGTAGTCGGTCAGGCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAAt                                                                                                                                               >  1:272370/1‑149 (MQ=255)
  gTTGATCTGGCGGTTGGAGGTACCGGTGGCAATGATCATGTAGTCGGTCAGGCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAAttt                                                                                                                                             >  6:143861/1‑149 (MQ=255)
         tGGCGGTTGGAGGTACCGGTGGCAATGATCATGTAGTCGGTCAGGCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCtt                                                                                                                                      >  3:82632/1‑149 (MQ=255)
           gCGGTTGGAGGTACCGGTTGCAATGATCATGTAGTCGGTCAGGCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGc                                                                                                                                        <  8:248326/145‑1 (MQ=255)
           gCGGTTGGAGGTACCGGTGGCAATGATCATGTAGTCGGTCAGGCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTgg                                                                                                                                    >  6:14335/1‑149 (MQ=255)
                       aCCGGTGGCAATGATCATGTAGTCGGTCAGGCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTNTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCTGTAAATTTTCTGCTTGGTCATATAAAACt                                                                                                                        >  3:134244/1‑149 (MQ=255)
                       aCCGGTGGCAATGATCATGTAGTCGGTCAGGCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCTGTAAATTTTCTGCTTGGTCATATAAAACt                                                                                                                        >  3:129345/1‑149 (MQ=255)
                        ccGGTGGCAATGATCATGTAGTCGGTCAGGCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTc                                                                                                                       <  5:14335/149‑1 (MQ=255)
                           gTGGCAATGATCATGTAGTCGGTCAGGCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGtt                                                                                                                    >  1:130157/1‑149 (MQ=255)
                                          tAGTCGGTCAGGCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCTGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGAt                                                                                                       <  1:242971/148‑1 (MQ=255)
                                             tCGGTCAGGCTGTGCTTCTCGCGCACNTCGATGACCNGGATNTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCTGTAAATTTTCTGCTTGGTCATATAANACTCGTTCAACTCGTTGGATGAGgc                                                                                                  <  4:134244/149‑1 (MQ=255)
                                             tCGGTCAGGCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCTGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGgc                                                                                                  <  4:129345/149‑1 (MQ=255)
                                                gTCAGGCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGgcgc                                                                                                >  5:71127/1‑148 (MQ=255)
                                                   aGGCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCTCTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCaga                                                                                            <  4:56843/149‑1 (MQ=255)
                                                     gCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCagag                                                                                           >  5:179517/1‑148 (MQ=255)
                                                     gCTGGGCCTCTCGCGCACGTCGATGACCTGTATGTCCTGGGCCTTTACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCagag                                                                                           <  6:179517/148‑1 (MQ=255)
                                                           cTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGt                                                                                    >  5:68071/1‑149 (MQ=255)
                                                            ttCTCGCGCACGTCGATGACCTGCATGTCCTGGGCATTGACGTCTTCCAGCGCTGACTTGGTCACTCCGACCTGTACTTCGCCGTAAATTTTCTGGTTGGTCATATTAAACTCGTTCAACTCGTTCGATGAGGCGCCAGAGGCACCGt                                                                                    >  2:186393/1‑148 (MQ=255)
                                                                 gcgcACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGtt                                                                               <  7:6206/148‑1 (MQ=255)
                                                                         cGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCAc                                                                       >  6:148861/1‑148 (MQ=255)
                                                                                   gctGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTc                                                             <  6:71127/146‑1 (MQ=255)
                                                                                     tGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTg                                                         >  8:102351/1‑149 (MQ=255)
                                                                                      gTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTc                                                                                                                       <  3:195554/87‑1 (MQ=255)
                                                                                      gTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTc                                                                                                                       >  4:195554/1‑87 (MQ=255)
                                                                                          tGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCt                                                    >  3:179501/1‑149 (MQ=255)
                                                                                          tGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCt                                                    >  3:26329/1‑149 (MQ=255)
                                                                                           gggCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTc                                                   >  7:276659/1‑149 (MQ=255)
                                                                                                ttGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCAc                                                                       <  3:152319/125‑1 (MQ=255)
                                                                                                ttGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCAc                                                                       >  4:152319/1‑125 (MQ=255)
                                                                                                 tGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGata                                               >  1:159139/1‑147 (MQ=255)
                                                                                                  gACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGatat                                              <  3:72497/147‑1 (MQ=255)
                                                                                                   aCGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAg                                            <  5:143861/148‑1 (MQ=255)
                                                                                                   aCGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAg                                            <  5:148861/148‑1 (MQ=255)
                                                                                                    cGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGc                                          >  1:176810/1‑149 (MQ=255)
                                                                                                    cGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGc                                          <  2:272370/149‑1 (MQ=255)
                                                                                                      tCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCa                                        <  2:159139/149‑1 (MQ=255)
                                                                                                       cTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCata                                                                                                                               >  4:46893/1‑62 (MQ=255)
                                                                                                       cTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCata                                                                                                                               <  3:46893/62‑1 (MQ=255)
                                                                                                       cTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGTGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAg                                       >  2:86264/1‑149 (MQ=255)
                                                                                                           cAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACg                                   >  2:163319/1‑149 (MQ=255)
                                                                                                            aGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTg                                                         >  1:96809/1‑126 (MQ=255)
                                                                                                            aGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTg                                                         <  2:96809/126‑1 (MQ=255)
                                                                                                                cTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAATTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTc                              >  4:99154/1‑149 (MQ=255)
                                                                                                                   ccTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTcc                             <  4:82632/147‑1 (MQ=255)
                                                                                                                       ggTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAATTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCAc                        <  3:99154/148‑1 (MQ=255)
                                                                                                                           aCTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAgg                    >  4:299859/1‑148 (MQ=255)
                                                                                                                           aCTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAgg                    >  1:21148/1‑148 (MQ=255)
                                                                                                                           aCTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAgg                    >  8:25502/1‑148 (MQ=255)
                                                                                                                            cTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCAc                                                                       <  2:159732/97‑1 (MQ=255)
                                                                                                                            cTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCAc                                                                       <  4:156003/97‑1 (MQ=255)
                                                                                                                            cTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGaa                  <  2:176810/149‑1 (MQ=255)
                                                                                                                                  ccAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGc                                 <  8:127270/128‑1 (MQ=255)
                                                                                                                                    aGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGata                                                                   >  6:33960/1‑93 (MQ=255)
                                                                                                                                    aGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCg          <  4:26329/149‑1 (MQ=255)
                                                                                                                                            cggcGTAAATTTTCTTCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGTGCGCAAGATATAGTTCTGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGAtt       <  1:86264/143‑1 (MQ=255)
                                                                                                                                             gCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCg  >  3:284815/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GCGTTGATCTGGCGGTTGGAGGTACCGGTGGCAATGATCATGTAGTCGGTCAGGCTGTGCTTCTCGCGCACGTCGATGACCTGGATGTCCTGGGCCTTGACGTCTTCCAGCGCTGCCTTGGTCACTGCGACCAGTTCTTCGCCGTAAATTTTCTGCTTGGTCATATAAAACTCGTTCAACTCGTTGGATGAGGCGCCAGAGGCACCGTCAGTTGGGCGCACGATATAGTTCGTGCGCCTCGATATAGGCCAGTACGGCGTCCGGCACCAGGAACCTCACCGATTTGCCG  >  NC_002947/5472655‑5472943

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: