Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 122,856 G→A *255* (TAG→TAA pdhR → pyruvate dehydrogenase complex repressor, autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913122,8560GA97.4% 115.5 / ‑6.3 38*255* (TAG→TAApdhRpyruvate dehydrogenase complex repressor, autorepressor
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (18/19);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAA  >  NC_000913/122727‑122989
                                                                                                                                 |                                                                                                                                     
tAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTc                                                                                                                              >  2:765336/1‑139 (MQ=255)
                  cgATGCATCGCATCGCCATCTGGCCATTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCaga                                                                                                            <  1:693622/139‑1 (MQ=255)
                  cgAAGCATCGCTTCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCaga                                                                                                            <  1:695923/139‑1 (MQ=255)
                  cgAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCaga                                                                                                            >  1:946581/1‑139 (MQ=255)
                         tcgcatcgcTATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGAtgtt                                                                                                     <  1:30162/139‑1 (MQ=255)
                                 ccATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTaaaaa                                                                                                                     <  1:36475/115‑1 (MQ=255)
                                 ccATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTaaaaa                                                                                                                     >  2:36475/1‑115 (MQ=255)
                                           tttATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCAtata                                                                                   >  2:533367/1‑139 (MQ=255)
                                                 gaagaaATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAgcg                                                                             <  1:227382/139‑1 (MQ=255)
                                                                  cAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTc                                                                                                                              >  1:307785/1‑73 (MQ=255)
                                                                  cAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTc                                                                                                                              <  2:307785/73‑1 (MQ=255)
                                                                    gAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAAc                                                          <  2:534511/139‑1 (MQ=255)
                                                                        tCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGt                                                      <  2:940101/139‑1 (MQ=255)
                                                                                  agCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGa                                                            >  2:750188/1‑123 (MQ=255)
                                                                                  agCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGa                                                            <  1:750188/123‑1 (MQ=255)
                                                                                  agCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGc                                            >  2:429459/1‑139 (MQ=255)
                                                                                  agCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGc                                            >  2:197749/1‑139 (MQ=255)
                                                                                  agCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGc                                            >  1:1044274/1‑139 (MQ=255)
                                                                                    ccgccgTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATc                                                                                                                <  2:508557/69‑1 (MQ=255)
                                                                                    ccgccgTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATc                                                                                                                >  1:508557/1‑69 (MQ=255)
                                                                                    ccgccgTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCtt                                          >  2:728013/1‑139 (MQ=255)
                                                                                            aGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCga                                  >  1:282577/1‑139 (MQ=255)
                                                                                                 tctcTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTa                                                                    <  1:704466/100‑1 (MQ=255)
                                                                                                 tctcTGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTa                                                                    >  2:704466/1‑100 (MQ=255)
                                                                                                   tctGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCa                           <  1:533367/139‑1 (MQ=255)
                                                                                                    ctGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAAc                                                          >  2:1060164/1‑107 (MQ=255)
                                                                                                    ctGCGTCGTCTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAAc                                                          <  1:1060164/107‑1 (MQ=255)
                                                                                                           gtcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAgc                                                                                          <  1:918247/68‑1 (MQ=255)
                                                                                                           gtcTGGAGCAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAgc                                                                                          >  2:918247/1‑68 (MQ=255)
                                                                                                                   cAACGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAg           <  1:58218/139‑1 (MQ=255)
                                                                                                                      cGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGtt                             >  2:99809/1‑118 (MQ=255)
                                                                                                                      cGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGtt                             <  1:99809/118‑1 (MQ=255)
                                                                                                                      cGAAAGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGaaa        >  1:902291/1‑139 (MQ=255)
                                                                                                                          aGAATTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGTGAGAGTTCAATGGGACAGGTTCCAGAAAACTc    >  2:159891/1‑139 (MQ=255)
                                                                                                                            aaTTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGTGAGAGTTCAATGGGACAGGTTCCAGAAAACTCaa  <  1:159891/139‑1 (MQ=255)
                                                                                                                            aaTTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCaa  <  2:839891/139‑1 (MQ=255)
                                                                                                                            aaTTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCaa  <  2:902291/139‑1 (MQ=255)
                                                                                                                            aaTTAATGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCaa  <  1:521568/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                     
TAAGCCGGAAGAAGCGCGCGAAGCATCGCATCGCCATCTGGCCTTTATCGAAGAAATTTTGCTCGACAGAAGTCGTGAAGAGAGCCGCCGTGAGCGTTCTCTGCGTCGTCTGGAGCAACGAAAGAATTAGTGATTTTTCTGGTAAAAATTATCCAGAAGATGTTGTAAATCAAGCGCATATAAAAGCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAA  >  NC_000913/122727‑122989

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: