Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,806,760 C→G 8.0% P459R (CCC→CGC)  PP_3364 → two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,806,7600CG 8.0% 153.6 / 4.2 50P459R (CCC→CGC) PP_3364two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base G (2/2);  total (25/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 = 3.24e-01

TCTGGTGTGCTGATTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGCGAACTGCTCGCCAGTGAACGCGC  >  NC_002947/3806618‑3806907
                                                                                                                                              |                                                                                                                                                   
tCTGGTGTGCTGATTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCaa                                                                                                                                               <  8:42117/149‑1 (MQ=255)
tCTGGTGTGCTGATTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCaa                                                                                                                                               <  7:32065/149‑1 (MQ=255)
 cTGTTGTGCTGACTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCaa                                                                                                                                               <  8:47319/148‑1 (MQ=255)
       tgCTGATTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTc                                                                                                                                        <  3:330149/149‑1 (MQ=255)
            attattCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCAt                                                                                                                                   >  8:173405/1‑149 (MQ=255)
               attCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCt                                                                                                                                 <  3:104827/148‑1 (MQ=255)
                           cAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGG‑CGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGtt                                                                                                                   <  2:166994/149‑1 (MQ=255)
                            aaCACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGtt                                                                                                                   >  6:39010/1‑149 (MQ=255)
                                 ccAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGcc                                                                                                                                                     >  5:131360/1‑110 (MQ=255)
                                 ccAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGcc                                                                                                                                                     <  6:131360/110‑1 (MQ=255)
                                  cAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCTATCTGGAAAGAAACCGTTCACcag                                                                                                             >  4:92637/1‑149 (MQ=255)
                                        aTTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCg                                                                                                       >  4:169003/1‑149 (MQ=255)
                                            cTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTAc                                                                                                   >  7:205785/1‑149 (MQ=255)
                                            cTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTAc                                                                                                   >  5:133504/1‑149 (MQ=255)
                                             ttGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTAcc                                                                                                  >  5:308124/1‑149 (MQ=255)
                                               gATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGgaaa                                                                                                                            >  2:228277/1‑121 (MQ=255)
                                               gATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGgaaa                                                                                                                            <  1:228277/121‑1 (MQ=255)
                                                      ttcCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGa                                                                                         >  8:148369/1‑149 (MQ=255)
                                                              ggaggTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGAcc                                                                                 <  6:133504/149‑1 (MQ=255)
                                                              ggaggTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGAcc                                                                                 <  8:205785/149‑1 (MQ=255)
                                                                    ggtCGAGACACTCGATTGGGCCGGGGACGCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACcagc                                                                                                            >  1:194861/1‑116 (MQ=255)
                                                                    ggtCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACcagc                                                                                                            <  2:194861/116‑1 (MQ=255)
                                                                      tCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGttt                                                                         >  7:325026/1‑149 (MQ=255)
                                                                        gagaCACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTgg                                                                       >  5:201152/1‑149 (MQ=255)
                                                                             acTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATg                                                                  >  2:148005/1‑149 (MQ=255)
                                                                                 gATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCAt                                                              <  7:254576/149‑1 (MQ=255)
                                                                                    tGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCg                                                           >  8:22349/1‑149 (MQ=255)
                                                                                        ccGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCAcc                                                       <  3:169003/149‑1 (MQ=255)
                                                                                            ggACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTAcc                                                                                                  <  4:206443/102‑1 (MQ=255)
                                                                                            ggACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTAcc                                                                                                  >  3:206443/1‑102 (MQ=255)
                                                                                              aCCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTAtc                                                 <  8:325026/149‑1 (MQ=255)
                                                                                                   aaCAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGa                                            <  6:308124/149‑1 (MQ=255)
                                                                                                        gACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCCGGACCGCCAGTTTGGCGATGCGATCCGCGCCGCGATCGTCGAggtgg                                       >  8:178512/1‑149 (MQ=255)
                                                                                                           cTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCAc                                                        >  5:282818/1‑129 (MQ=255)
                                                                                                           cTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCAc                                                        <  6:282818/129‑1 (MQ=255)
                                                                                                           cTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGtgc                                    >  4:320612/1‑149 (MQ=255)
                                                                                                           cTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGtgc                                    >  3:26480/1‑149 (MQ=255)
                                                                                                            tATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGc                                                                                       <  2:40776/97‑1 (MQ=255)
                                                                                                            tATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGc                                                                                       >  1:40776/1‑97 (MQ=255)
                                                                                                              ctcttgcggggCGTCAGGCGGTCGCCCAACGCGCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCg                                                      <  1:13082/123‑1 (MQ=37)
                                                                                                              ctcttgcggggCGTCAGGCGGTCGCCCAACGCGCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCg                                                      >  2:13082/6‑128 (MQ=37)
                                                                                                                aCAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGca                                <  7:148369/148‑1 (MQ=255)
                                                                                                                aCAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGca                                >  6:156501/1‑148 (MQ=255)
                                                                                                                aCAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGca                                >  4:214263/1‑148 (MQ=255)
                                                                                                                  tgcggggCGTCAGGCGGTCGCCCAACGCGCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGc                                                                      >  2:245933/2‑108 (MQ=37)
                                                                                                                  tgcggggCGTCAGGCGGTCGCCCAACGCGCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGc                                                                      <  1:245933/107‑1 (MQ=37)
                                                                                                                    cggcgCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACaa                             >  5:52823/1‑147 (MQ=255)
                                                                                                                     ggcgCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAg                          <  3:320612/149‑1 (MQ=255)
                                                                                                                      gcgcGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGc                         <  7:173405/149‑1 (MQ=255)
                                                                                                                      gcgcGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGc                         <  4:26480/149‑1 (MQ=255)
                                                                                                                         cgTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGCGaa                      <  5:39010/149‑1 (MQ=255)
                                                                                                                                   ccGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGCGAACTGCTCGCCa            <  6:201152/149‑1 (MQ=255)
                                                                                                                                              cccGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCCCCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTAGGCGATGCCATCCGCACCGCTATCGTCGAGGGGGTGCTGCACAATAGCGAACTGCTCTCCACTGAACgcgc  >  6:273314/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TCTGGTGTGCTGATTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGCGAACTGCTCGCCAGTGAACGCGC  >  NC_002947/3806618‑3806907

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: ‑