Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,352,104 G→A 8.1% V122V (GTG→GTA)  PP_2947 → transcriptional regulator MvaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,352,1040GA 8.1% 137.3 / 2.5 49V122V (GTG→GTA) PP_2947transcriptional regulator MvaT
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base A (2/2);  total (24/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 = 8.04e-01

TACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGG  >  NC_002947/3351957‑3352240
                                                                                                                                                    |                                                                                                                                        
tACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGgtgt                                                                                                                                         >  6:194724/1‑149 (MQ=255)
tACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGgtgt                                                                                                                                         >  2:198835/1‑149 (MQ=255)
      gaAAGCCCCACGCCTTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGATAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGa                                                                                                                                    <  3:261440/148‑1 (MQ=255)
       aaaGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGAtt                                                                                                                                  >  4:68909/1‑149 (MQ=255)
        aaGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGAttt                                                                                                                                 <  5:194724/149‑1 (MQ=255)
          gCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTtgt                                                                                                                               <  3:112737/149‑1 (MQ=255)
           ccccACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCATGATTtgtg                                                                                                                              >  2:198787/1‑149 (MQ=255)
            cccAAGCCGTGCCCTTCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAAGCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGTGG‑CAGAAGAAGTGGAATCCTTGGTGTCCTGATTTGTGc                                                                                                                             <  4:81090/149‑1 (MQ=255)
                           gCAAGTGAAGCATTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAATGAATGGAAAGCCCAGTACGGGGCCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCa                                                                                                               <  1:198787/149‑1 (MQ=255)
                                  aaGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTcc                                                                                                                                       >  3:137397/1‑117 (MQ=255)
                                  aaGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTcc                                                                                                                                       <  4:137397/117‑1 (MQ=255)
                                     cAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGAAATCCTGGGTATCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCt                                                                                                    >  7:6149/1‑149 (MQ=255)
                                     cAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGAAATCCTGGGTATCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCt                                                                                                    >  7:10207/1‑149 (MQ=255)
                                                 ccGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGt                                                                                        >  8:159358/1‑149 (MQ=255)
                                                   gCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCGCTCATTGCGCTATCTGAAACGCCTTGTTc                                                                                      >  8:188257/1‑149 (MQ=255)
                                                                 tCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCt                                                                                           >  1:174540/1‑130 (MQ=255)
                                                                 tCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCt                                                                                           <  2:174540/130‑1 (MQ=255)
                                                                 tCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGc                                                                        <  3:155394/149‑1 (MQ=255)
                                                                           cAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCg                                                              >  3:87480/1‑149 (MQ=255)
                                                                           cAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGCAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCACTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCg                                                              >  1:31647/1‑149 (MQ=255)
                                                                              gggTGGAAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTAATGATTTGTGCTGATGTATCCCTCATTGCGCTATCTGAAACGACTTGTTCTATTTCTAGGAATCGAAAGCGTCGcct                                                           >  3:139746/1‑149 (MQ=255)
                                                                                   gCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAGGCGACGCTGTCGCttcctc                                                       <  3:68909/148‑1 (MQ=255)
                                                                                       ccACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGAAATCCTGGGTATCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAGGCGACGCTGTCGCTTCCTCGaaaa                                                  <  8:10207/149‑2 (MQ=255)
                                                                                       ccACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGAAATCCTGGGTATCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAGGCGACGCTGTCGCTTCCTCGaaaa                                                  <  8:6149/149‑2 (MQ=255)
                                                                                        cacaAGCTGCTCAAGGAATGGAAAGCCCAGTACGGTG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGcc                                                                                            <  6:232442/106‑1 (MQ=255)
                                                                                        cacaAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGcc                                                                                            >  5:232442/1‑106 (MQ=255)
                                                                                          caAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAATGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAgg                                               <  1:198835/149‑1 (MQ=255)
                                                                                          caAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCACCCTCATTGCGCTATCTGAAACGCCTTGTTCttt                                                                                   >  5:202631/1‑113 (MQ=255)
                                                                                          caAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCACCCTCATTGCGCTATCTGAAACGCCTTGTTCttt                                                                                   <  6:202631/113‑1 (MQ=255)
                                                                                                 ctCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGaaa                                                   >  6:232237/1‑138 (MQ=255)
                                                                                                 ctCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGaaa                                                   <  5:232237/138‑1 (MQ=255)
                                                                                                      ggAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgtg                                   >  8:185391/1‑149 (MQ=255)
                                                                                                      ggAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgtg                                   >  7:119367/1‑149 (MQ=255)
                                                                                                       gAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgtgg                                  <  7:185391/149‑1 (MQ=255)
                                                                                                           ggAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCg                                            <  5:33895/135‑1 (MQ=255)
                                                                                                           ggAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCg                                            >  6:33895/1‑135 (MQ=255)
                                                                                                            gAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGtgggtggg                             >  7:139930/1‑149 (MQ=255)
                                                                                                            gAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGtgggtggg                             <  7:159358/149‑1 (MQ=255)
                                                                                                              aaGCCCAGTACGGTG‑GAGAAGAAGTGGAATCGTTGGTGTTCTGATTTGTGCTTATGTCTCCTTCATTGCGCTATATGAAACTGCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAATCATTGGGTTCATCGTGGGTGGGc                            <  4:139746/148‑1 (MQ=255)
                                                                                                              aaGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCACTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGc                            <  2:31647/148‑1 (MQ=255)
                                                                                                                  ccAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCTTTCATATTGGGTGGGCGGAAt                       >  7:202266/1‑149 (MQ=255)
                                                                                                                  ccAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAAt                       <  8:119367/149‑1 (MQ=255)
                                                                                                                               agaagaAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCg                           <  8:9486/133‑1 (MQ=255)
                                                                                                                               agaagaAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCg                           >  7:9486/1‑133 (MQ=255)
                                                                                                                               agaagaAGTGGAATCCTGGCTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTgg           >  2:40043/1‑149 (MQ=255)
                                                                                                                                 aagaagTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGt          >  7:128229/1‑148 (MQ=255)
                                                                                                                                      gTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACa    <  8:128229/149‑1 (MQ=255)
                                                                                                                                       tGGAATCCTGGCTGTCCTGATTTGTGCTGATGACGCCATCAATGCGATATCTTAAACCCCTTCTTCTTTTTAGACGAAGCGAAAGATTCGCCTCCTCGAAACAGGGCGTTAATCGTGCGTGGGCGGAATCAGCTTACGATGTTAAACAg   >  2:183416/1‑149 (MQ=255)
                                                                                                                                         gAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCgg                          >  6:31389/1‑124 (MQ=255)
                                                                                                                                         gAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCgg                          <  5:31389/124‑1 (MQ=255)
                                                                                                                                         gAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAgg  <  8:139930/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                        
TACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGG‑CAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGG  >  NC_002947/3351957‑3352240

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑