Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,806,754 T→C 7.5% L457P (CTG→CCG) ‡ 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,7540TC 7.5% 151.0 / 2.5 53L457P (CTG→CCG) ‡PP_3364two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base C (2/2);  total (26/27)
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.79e-01

AATGATGTGTCTGGTGTGCTGATTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGCGAACTGCTCGCCA  >  NC_002947/3806609‑3806897
                                                                                                                                                 |                                                                                                                                               
aaTGATGTGTCTGGTGTGCTGATTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGa                                                                                                                                               <  7:79658/148‑1 (MQ=255)
         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)
                  cTGATTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCt                                                                                                                                                 >  7:127474/1‑128 (MQ=255)
                  cTGATTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCt                                                                                                                                                 <  8:127474/128‑1 (MQ=255)
                     attattCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCAt                                                                                                                         >  8:173405/1‑149 (MQ=255)
                        attCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCt                                                                                                                       <  3:104827/148‑1 (MQ=255)
                                    cAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGtt                                                                                                         <  2:166994/149‑1 (MQ=255)
                                     aaCACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGtt                                                                                                         >  6:39010/1‑149 (MQ=255)
                                          ccAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGcc                                                                                                                                           <  6:131360/110‑1 (MQ=255)
                                          ccAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGcc                                                                                                                                           >  5:131360/1‑110 (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                                                                       <  8:205785/149‑1 (MQ=255)
                                                                       ggaggTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGAcc                                                                       <  6:133504/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                                              <  6:282818/129‑1 (MQ=255)
                                                                                                                    cTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCAc                                              >  5:282818/1‑129 (MQ=255)
                                                                                                                    cTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGtgc                          >  3:26480/1‑149 (MQ=255)
                                                                                                                    cTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGtgc                          >  4:320612/1‑149 (MQ=255)
                                                                                                                     tATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGc                                                                             >  1:40776/1‑97 (MQ=255)
                                                                                                                     tATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGc                                                                             <  2:40776/97‑1 (MQ=255)
                                                                                                                       ctcttgcggggCGTCAGGCGGTCGCCCAACGCGCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCg                                            >  2:13082/6‑128 (MQ=37)
                                                                                                                       ctcttgcggggCGTCAGGCGGTCGCCCAACGCGCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCg                                            <  1:13082/123‑1 (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)
                                                                                                                                                 |                                                                                                                                               
AATGATGTGTCTGGTGTGCTGATTATTCCGATGTCGCAACACCCAAGGGATTACTTGATCTTCTTCCGCAAGGAGGTGGTCGAGACACTCGATTGGGCCGGGGACCCGAACAAGACCTATGACAGCGGCGCGTTGGGCGACCGCCTGACGCCCCGCAAGAGCTTCGCCATCTGGAAAGAAACCGTTCACCAGCAGTCGCTACCCTGGACCGAGCAGGACCGCCAGTTTGGCGATGCCATCCGCACCGCTATCGTCGAGGTGGTGCTGCACAATAGCGAACTGCTCGCCA  >  NC_002947/3806609‑3806897

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: