Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,428,061 C→G 8.9% S670S (TCC→TCG) ‡ PP_3042 → terminase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,428,0610CG 8.9% 171.5 / 7.1 56S670S (TCC→TCG) ‡PP_3042terminase large subunit
Reads supporting (aligned to +/- strand):  ref base C (24/27);  new base G (3/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

GCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGACGAACTGATCAAG  >  NC_002947/3427914‑3428206
                                                                                                                                                   |                                                                                                                                                 
gCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCca                                                                                                                                                  <  4:5426/149‑1 (MQ=255)
  aaCCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCcagc                                                                                                                                                <  1:93042/149‑1 (MQ=255)
   aCCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCcag                                                                                                                                                 >  2:7435/1‑147 (MQ=255)
       gTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCgg                                                                                                                                           <  1:7435/149‑1 (MQ=255)
               gCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTg                                                                                                                                    >  6:301807/1‑148 (MQ=255)
                cACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTg                                                                                                                                    >  6:278171/1‑147 (MQ=255)
                 acgaaCACGATGACGACGATGACCCGGACGACGATCCAGTTTCAGCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGaa                                                                                                                                  <  1:142249/145‑1 (MQ=255)
                  ccAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAaga                                                                                                                                <  5:278171/149‑1 (MQ=255)
                  ccAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAaga                                                                                                                                <  4:150110/149‑1 (MQ=255)
                    aaCACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAgaga                                                                                                                              <  2:313871/149‑1 (MQ=255)
                     acacGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCTGCTGGAGCTGCGGcgtt                                                                                                                                       >  2:172643/1‑136 (MQ=255)
                     acacGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCTGCTGGAGCTGCGGcgtt                                                                                                                                       <  1:172643/139‑4 (MQ=255)
                         gATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCt                                                                                                                          >  7:67181/1‑148 (MQ=255)
                               gacgaTGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATgg                                                                                                                   >  3:10629/1‑149 (MQ=255)
                                  gaTGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCa                                                                                                                 >  2:171961/1‑148 (MQ=255)
                                        ccGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGaa                                                                                                                                  >  8:86598/1‑125 (MQ=255)
                                        ccGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGaa                                                                                                                                  <  7:86598/125‑1 (MQ=255)
                                              gacgaTCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTcc                                                                                                                                                   >  8:72904/1‑102 (MQ=255)
                                              gacgaTCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTcc                                                                                                                                                   <  7:72904/102‑1 (MQ=255)
                                                 gaTCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCacac                                                                                                  >  3:302321/1‑148 (MQ=255)
                                                   tCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGACTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTc                                                                                               >  2:188625/1‑149 (MQ=255)
                                                       gTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCgg                                                                                           <  5:301807/149‑1 (MQ=255)
                                                       gTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCgg                                                                                           <  8:67181/149‑1 (MQ=255)
                                                            aCCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTc                                                                                      <  2:200314/149‑1 (MQ=255)
                                                              cGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGa                                                                                    <  1:171961/149‑1 (MQ=255)
                                                                 cAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCt                                                                                                                          >  4:43746/1‑108 (MQ=255)
                                                                 cAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCt                                                                                                                          <  3:43746/108‑1 (MQ=255)
                                                                 cAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCa                                                                                                                 <  4:43497/117‑1 (MQ=255)
                                                                  aaCACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCTGCTGGAGCTGCGGc                                                                                                                                          <  4:33512/91‑1 (MQ=255)
                                                                  aaCACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCTGCTGGAGCTGCGGc                                                                                                                                          >  3:33512/1‑91 (MQ=255)
                                                                  aaCACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGt                                                                                >  6:216872/1‑149 (MQ=255)
                                                                    caccacCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGc                                                                              >  8:234283/1‑149 (MQ=255)
                                                                    caccacCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGc                                                                              <  4:302321/149‑1 (MQ=255)
                                                                                        ccacgctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAAc                                                             >  8:149073/1‑146 (MQ=255)
                                                                                        ccacgctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAAc                                                             <  7:149073/146‑1 (MQ=255)
                                                                                         cacgctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGgt                                                         >  2:129353/1‑149 (MQ=255)
                                                                                            gctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCata                                                                                                               <  7:94210/92‑1 (MQ=255)
                                                                                            gctccacTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCata                                                                                                               >  8:94210/1‑92 (MQ=255)
                                                                                              tccactccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGt                                                                    <  7:320567/133‑1 (MQ=255)
                                                                                              tccactccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGt                                                                    >  8:320567/1‑133 (MQ=255)
                                                                                                cactccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGa                                                  <  7:234283/149‑1 (MQ=255)
                                                                                                cactccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGCCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGa                                                  >  5:303808/1‑149 (MQ=255)
                                                                                                cactccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCAGATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGa                                                  >  5:157329/1‑149 (MQ=255)
                                                                                                 actccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAACGTGCGATTGCGCGTGGCGAACGGgttgtt                                                     >  6:45457/1‑145 (MQ=255)
                                                                                                 actccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAACGTGCGATTGCGCGTGGCGAACGGgttgtt                                                     <  5:45457/145‑1 (MQ=255)
                                                                                                  ctccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATa                                                <  1:129353/149‑1 (MQ=255)
                                                                                                  cacccccgccgCCTCGCCCCGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGACACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATa                                                <  4:10629/144‑1 (MQ=255)
                                                                                                          ccgccTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTgcgc                                                                      >  3:310482/1‑119 (MQ=255)
                                                                                                          ccgccTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTgcgc                                                                      <  4:310482/119‑1 (MQ=255)
                                                                                                             ccTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGAc                                         >  3:303120/1‑145 (MQ=255)
                                                                                                             ccTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGAc                                         <  4:303120/145‑1 (MQ=255)
                                                                                                                    cGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGa                              <  6:303808/149‑1 (MQ=255)
                                                                                                                          cTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCg                        >  6:172209/1‑149 (MQ=255)
                                                                                                                           tCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCg                        >  5:220378/1‑148 (MQ=255)
                                                                                                                            ccACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCa                                                                                                                 <  8:300178/58‑1 (MQ=255)
                                                                                                                            ccACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCa                                                                                                                 >  7:300178/1‑58 (MQ=255)
                                                                                                                                 cATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACg                    >  8:212080/1‑146 (MQ=255)
                                                                                                                                 cATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACg                    <  7:212080/146‑1 (MQ=255)
                                                                                                                                     cAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGAc              <  5:216872/148‑1 (MQ=255)
                                                                                                                                      atagccgcTGCTGGAGCTGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGACg             >  6:158058/4‑148 (MQ=255)
                                                                                                                                                cTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGACGAACTGATCAAg  <  6:157329/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCAACCGGTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGACGAACTGATCAAG  >  NC_002947/3427914‑3428206

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: