Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,388,284 A→T 20.0% intergenic (+23/+21) tyrR → / ← tpx aromatic amino acid biosynthesis and transport regulon transcriptional regulator; autorepressor; ATPase; phosphatase/lipid hydroperoxide peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,388,2840AT20.0% 57.4 / 10.5 25intergenic (+23/+21)tyrR/tpxaromatic amino acid biosynthesis and transport regulon transcriptional regulator; autorepressor; ATPase; phosphatase/lipid hydroperoxide peroxidase
Reads supporting (aligned to +/- strand):  ref base A (10/10);  new base T (2/3);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

CCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTT  >  NC_000913/1388141‑1388430
                                                                                                                                               |                                                                                                                                                  
ccAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAacac                                                                                                                                                <  4:312436/148‑1 (MQ=255)
  aGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCaa                                                                                                                                                   >  4:494013/1‑143 (MQ=255)
  aGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCaa                                                                                                                                                   <  3:494013/143‑1 (MQ=255)
            cAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCa                                                                                                                                           <  1:1098110/141‑1 (MQ=255)
            cAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATCa                                                                                                                                           >  2:1098110/1‑141 (MQ=255)
                      aGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCATc                                                                                                                                            <  4:246310/130‑1 (MQ=255)
                                aaCTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCa                                                                                                                                           >  1:346800/1‑121 (MQ=255)
                                aaCTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCa                                                                                                                                           <  2:346800/121‑1 (MQ=255)
                                           cGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCa                                                                                                                                                    >  4:833045/1‑101 (MQ=255)
                                           cGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCa                                                                                                                                                    <  3:833045/101‑1 (MQ=255)
                                                      ttCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTGCACCAt                                                                                                                                             >  3:441145/1‑97 (MQ=255)
                                                                                                                                      tGATGGTGTTGCACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaaaa                                                      <  3:218258/104‑1 (MQ=255)
                                                                                                                                      tGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCAtttt                                       <  3:718324/119‑1 (MQ=255)
                                                                                                                                      tGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCAtttt                                       >  4:904433/1‑119 (MQ=255)
                                                                                                                                      tGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCAtttt                                       >  4:718324/1‑119 (MQ=255)
                                                                                                                                      tGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCAtttt                                       <  3:904433/119‑1 (MQ=255)
                                                                                                                                      tGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCAtttt                                       <  1:1056930/119‑1 (MQ=255)
                                                                                                                                      tGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCAtttt                                       >  2:1056930/1‑119 (MQ=255)
                                                                                                                                        aTGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTc                                                                                      <  3:526112/70‑1 (MQ=255)
                                                                                                                                        aTGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTc                                                                                      >  4:526112/1‑70 (MQ=255)
                                                                                                                                        aTGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAg                                                                 >  4:306735/1‑91 (MQ=255)
                                                                                                                                        aTGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAg                                                                 <  3:306735/91‑1 (MQ=255)
                                                                                                                                         tGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCa                                                                         <  3:793273/82‑1 (MQ=255)
                                                                                                                                         tGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCa                                                                         >  4:793273/1‑82 (MQ=255)
                                                                                                                                          ggTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc                                                            <  2:806143/94‑1 (MQ=255)
                                                                                                                                          ggTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc                                                            >  1:806143/1‑94 (MQ=255)
                                                                                                                                             gCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaa                                                        >  3:927189/1‑95 (MQ=255)
                                                                                                                                             gCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaa                                                        <  4:927189/95‑1 (MQ=255)
                                                                                                                                              cAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCttt  >  1:350500/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAACACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTT  >  NC_000913/1388141‑1388430

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