Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,428,065 A→T 8.8% S672C (AGC→TGC)  PP_3042 → terminase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,428,0650AT 8.8% 161.9 / 3.2 56S672C (AGC→TGC) PP_3042terminase large subunit
Reads supporting (aligned to +/- strand):  ref base A (24/27);  new base T (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 = 1.00e+00

GTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGACGAACTGATCAAG  >  NC_002947/3427921‑3428206
                                                                                                                                                |                                                                                                                                             
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                                                                                                                                  <  7:86598/125‑1 (MQ=255)
                                 ccGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGaa                                                                                                                                  >  8:86598/1‑125 (MQ=255)
                                          gaTCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCacac                                                                                                  >  3:302321/1‑148 (MQ=255)
                                            tCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGACTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTc                                                                                               >  2:188625/1‑149 (MQ=255)
                                                gTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCgg                                                                                           <  8:67181/149‑1 (MQ=255)
                                                gTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCgg                                                                                           <  5:301807/149‑1 (MQ=255)
                                                     aCCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTc                                                                                      <  2:200314/149‑1 (MQ=255)
                                                       cGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGa                                                                                    <  1:171961/149‑1 (MQ=255)
                                                          cAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCt                                                                                                                          <  3:43746/108‑1 (MQ=255)
                                                          cAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCt                                                                                                                          >  4:43746/1‑108 (MQ=255)
                                                          cAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCa                                                                                                                 <  4:43497/117‑1 (MQ=255)
                                                           aaCACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCTGCTGGAGCTGCGGc                                                                                                                                          >  3:33512/1‑91 (MQ=255)
                                                           aaCACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCTGCTGGAGCTGCGGc                                                                                                                                          <  4:33512/91‑1 (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                                                                    >  8:320567/1‑133 (MQ=255)
                                                                                       tccactccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGt                                                                    <  7:320567/133‑1 (MQ=255)
                                                                                         cactccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGa                                                  <  7:234283/149‑1 (MQ=255)
                                                                                         cactccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGCCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGa                                                  >  5:303808/1‑149 (MQ=255)
                                                                                         cactccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCAGATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGa                                                  >  5:157329/1‑149 (MQ=255)
                                                                                          actccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAACGTGCGATTGCGCGTGGCGAACGGgttgtt                                                     <  5:45457/145‑1 (MQ=255)
                                                                                          actccacCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAACGTGCGATTGCGCGTGGCGAACGGgttgtt                                                     >  6:45457/1‑145 (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                                                                                                                 >  7:300178/1‑58 (MQ=255)
                                                                                                                     ccACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCa                                                                                                                 <  8:300178/58‑1 (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)
                                                                                                                                                |                                                                                                                                             
GTCGTTGAGCACCAACACGATGACGACGATGACCCGGACGACGATCCAGTTTCACCGCCAACACCACCGGCTCGGCCCGCTCCACGCTCCACTCCACCGCCGCCTCGCCCGGCACCTCCACCCATGCAACGCCGCAGCTCCAGCAGCGGCTATCTGAAGAGACGCTGACATGGCATATACCCAGGCACACCTCGCGGCCGTCGAGCGTGCGATTGCGCGTGGCGAACGGGTTGTTCGATACAGCGACCGCACCGTGGAGTATCGGACGGTAGACGAACTGATCAAG  >  NC_002947/3427921‑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: ‑