Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,351,974 G→A 8.0% R79H (CGC→CAC)  PP_2947 → transcriptional regulator MvaT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,351,9740GA 8.0% 136.5 / 2.1 50R79H (CGC→CAC) PP_2947transcriptional regulator MvaT
Reads supporting (aligned to +/- strand):  ref base G (23/22);  new base A (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

TGAGCTTCAGAAAGAGATCGAGTTCGAGACCAAGTTGCGTGACCTGCTAAGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTG  >  NC_002947/3351828‑3352119
                                                                                                                                                   |                                                                                                                                                 
tGAGCTTCAGAAAGAGATCGAGTTCGAGACCAAGTTGCGTGACCTGCTAAGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGcc                                                                                                                                                 <  3:109997/149‑1 (MQ=255)
    cTTCAGAAAGAGATCGAGTTCGAGACCAAGTTGCGTGACCTGCTAAGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGc                                                                                                                                             <  4:191000/149‑1 (MQ=255)
            agagATCGAGTTCGAGACCAAGTTGCGTGACCTGCTAAGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAg                                                                                                                                     <  2:324490/149‑1 (MQ=255)
                aTCGAGTTCGAGACCAAGTTGCGTGACCTGCTAAGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCaa                                                                                                                                      <  1:134445/144‑1 (MQ=255)
                aTCGAGTTCGAGACCAAGTTGCGTGACCTGCTAAGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCaa                                                                                                                                      >  2:134445/1‑144 (MQ=255)
                   gAGTTCGAGACCAAGTTGCGTGACCTGCTAAGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCa                                                                                                                              <  1:6081/149‑1 (MQ=255)
                         gagaCCAAGTTGCGTGACCTGCTAAGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGt                                                                                                                                    >  5:94019/1‑137 (MQ=255)
                         gagaCCAAGTTGCGTGACCTGCTAAGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGt                                                                                                                                    <  6:94019/137‑1 (MQ=255)
                           gaCCAAGTTGCGTGACCTGCTAAGCAAATACGGTTACAGCCTGCGCGTATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACACCGTGCCCGGCAAGTGAAGCAGTAc                                                                                                                          >  1:259905/1‑146 (MQ=255)
                               aaGTTGCGTGACCTGCTAAGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAAccc                                                                                                                  <  1:330312/149‑1 (MQ=255)
                                       tGACCTGCTAAGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATg                                                                                                          >  6:134313/1‑149 (MQ=255)
                                           cTGCTAAGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACACCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGga                                                                                                        >  3:15465/1‑147 (MQ=255)
                                             gCTAAGCAAATACGGTTAAAGACTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCACTCGCCGTTCCGCAGCTCCTGAGGCTACAGAGAAAGCCCCACGCCTTGCCCGCCAAGTGAAGCAGTACAAGAACCCCCAATATGGAGAgg                                                                                                    >  5:63125/1‑149 (MQ=255)
                                              cTAAGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACACCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGagag                                                                                                     <  2:259905/147‑1 (MQ=255)
                                                aaGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGgaga                                                                                                      <  7:256181/144‑1 (MQ=255)
                                                aaGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGgaga                                                                                                      >  8:256181/1‑144 (MQ=255)
                                                 aGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCAt                                                                                                >  2:308082/1‑149 (MQ=255)
                                                  gCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATc                                                                                               >  6:3928/1‑149 (MQ=255)
                                                  gCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATc                                                                                               <  6:114840/149‑1 (MQ=255)
                                                          ggTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGAcc                                                                                         >  3:303068/1‑147 (MQ=255)
                                                              aCAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAggg                                                                                    <  7:119710/148‑1 (MQ=255)
                                                                      cgcg‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAAc                                                                             <  5:134313/147‑1 (MQ=255)
                                                                         g‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACaa                                                                        <  1:308082/149‑1 (MQ=255)
                                                                           atatCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAg                                                                       <  5:3928/149‑1 (MQ=255)
                                                                                atcaACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAAcc                                                                            <  8:30186/139‑1 (MQ=255)
                                                                                atcaACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAAcc                                                                            >  7:30186/1‑139 (MQ=255)
                                                                                 tcaACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAgctgct                                                                  >  6:41050/1‑148 (MQ=255)
                                                                                   aaCATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACTCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCGACCACAAGCTGCTc                                                                 >  2:214930/1‑147 (MQ=255)
                                                                                       tCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCTTCGAGACCAAGGGTGGCAACAACAAGCTGCGCAAGGaa                                                           >  6:251291/1‑149 (MQ=255)
                                                                                         cTGGATGCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAATGCCCCACGCCGTGCCCTGCAAGTGAAGCAGTACAAGAACCAGCATAATGGAGAGCTCATCGAGAACAAGGGTGGCAACCACAAGCTGCTCAAGGAAt                                                          <  6:63125/148‑1 (MQ=255)
                                                                                          tGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATg                                                         >  2:51601/1‑148 (MQ=255)
                                                                                                 tCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCga                                                                                             >  8:320296/1‑105 (MQ=255)
                                                                                                 tCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCga                                                                                             <  7:320296/105‑1 (MQ=255)
                                                                                                 tCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCTAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGcc                                                  >  8:166606/1‑148 (MQ=255)
                                                                                                 tCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGcc                                                  >  6:53814/1‑148 (MQ=255)
                                                                                                 tCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAACAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGcc                                                  >  3:175039/1‑148 (MQ=255)
                                                                                                            gCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCa                                      <  1:51601/149‑1 (MQ=255)
                                                                                                              cGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACACCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCaga                                    <  4:15465/149‑1 (MQ=255)
                                                                                                                        cTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGaa                          >  6:103989/1‑149 (MQ=255)
                                                                                                                         tCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAAt                         >  3:198854/1‑149 (MQ=255)
                                                                                                                                      gagaAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCt             >  1:72615/1‑148 (MQ=255)
                                                                                                                                       agaAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGa           >  4:99867/1‑149 (MQ=255)
                                                                                                                                        gaAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGAt          <  5:41050/149‑1 (MQ=255)
                                                                                                                                         aaaGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGGATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGcc                                                  <  2:298530/108‑1 (MQ=255)
                                                                                                                                         aaaGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGGATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGcc                                                  >  1:298530/1‑108 (MQ=255)
                                                                                                                                            gCCCCACGCCGTGCCCGGCAAGTGAAACAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTtgt      <  4:175039/149‑1 (MQ=255)
                                                                                                                                              cccACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTgg                    >  5:201261/1‑133 (MQ=255)
                                                                                                                                              cccACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTgg                    <  6:201261/133‑1 (MQ=255)
                                                                                                                                               ccACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTTCGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCt   >  7:314845/1‑149 (MQ=255)
                                                                                                                                                cACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTg  <  5:103989/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGAGCTTCAGAAAGAGATCGAGTTCGAGACCAAGTTGCGTGACCTGCTAAGCAAATACGGTTACAGCCTGCGCG‑ATATCATCAACATCCTGGATCCTCAGGCCAGTCGCCGTTCCGCAGCTCCTGCGGCTACAGAGAAAGCCCCACGCCGTGCCCGGCAAGTGAAGCAGTACAAGAACCCGCATAATGGAGAGGTCATCGAGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTG  >  NC_002947/3351828‑3352119

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