Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,512,601 A→G 9.5% intergenic (‑257/+33) yhfZ ← / ← trpS putative DNA‑binding transcriptional regulator/tryptophanyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,512,6010AG 9.5% 120.2 / 4.7 42intergenic (‑257/+33)yhfZ/trpSputative DNA‑binding transcriptional regulator/tryptophanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (19/19);  new base G (2/2);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

CTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCTTTTTCCGCGCCATCTTTCATCACCTGTTGCAGGAA  >  NC_000913/3512526‑3512738
                                                                           |                                                                                                                                         
cTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAAt                                                                                                                 <  4:1115606/102‑1 (MQ=255)
cTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAAt                                                                                                                 <  4:1048351/102‑1 (MQ=255)
cTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAAt                                                                                                                 >  3:1048351/1‑102 (MQ=255)
cTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAAt                                                                                                                 >  3:1115606/1‑102 (MQ=255)
  tATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGt                                                                           >  2:534642/1‑138 (MQ=255)
  tATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGt                                                                           <  1:534642/138‑1 (MQ=255)
       ggCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAg                                                                                                               >  2:662178/1‑97 (MQ=255)
       ggCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAg                                                                                                               <  1:662178/97‑1 (MQ=255)
         ccTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCtt                                                                                                  <  4:302911/108‑1 (MQ=255)
         ccTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCtt                                                                                                  >  3:302911/1‑108 (MQ=255)
          cTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGc                                                                                               <  2:354492/110‑1 (MQ=255)
          cTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGc                                                                                               >  1:354492/1‑110 (MQ=255)
           tACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAAt                                                                                                                 >  3:142285/1‑91 (MQ=255)
           tACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAAt                                                                                                                 <  4:142285/91‑1 (MQ=255)
                              tgCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCg                                                                                                >  3:155805/1‑89 (MQ=255)
                              tgCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCg                                                                                                <  4:155805/89‑1 (MQ=255)
                                              ccGGATGCGGCGTGAACGCCTTATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCtt                                                                              >  1:128731/1‑91 (MQ=25)
                                              ccGGATGCGGCGTGAACGCCTTATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCtt                                                                              <  2:128731/91‑1 (MQ=25)
                                              ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCtttt                                                                >  2:884630/1‑105 (MQ=255)
                                              ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCtttt                                                                <  1:884630/105‑1 (MQ=255)
                                              ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAg                                                              >  4:577440/1‑107 (MQ=255)
                                              ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAg                                                              <  3:577440/107‑1 (MQ=255)
                                                  aTAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGa                                                                                                            <  2:377978/57‑1 (MQ=255)
                                                  aTAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGa                                                                                                            >  1:377978/1‑57 (MQ=255)
                                                       ttgTTCATGCCGGATGCGGCGTGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCtttt                                  >  3:812419/3‑126 (MQ=255)
                                                       ttgTTCATGCCGGATGCGGCGTGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCtttt                                  <  4:812419/124‑1 (MQ=255)
                                                       gCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTa                                                          >  3:898656/1‑102 (MQ=255)
                                                       gCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTa                                                          <  4:898656/102‑1 (MQ=255)
                                                        cGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCg                                                 >  3:359258/1‑110 (MQ=255)
                                                        cGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTgg                                       >  1:1041316/1‑120 (MQ=255)
                                                        cGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTgg                                       <  2:1041316/120‑1 (MQ=255)
                                                        cGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTGGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCg                                                 <  4:359258/110‑1 (MQ=255)
                                                           tcatgccggATGCGGCGTGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCTTTTTCCGCGCCATCTTTCATCACCTGTTg        >  2:806981/5‑148 (MQ=255)
                                                                ccgcATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGTCACAAAACTAATCGCTTCGTACCCCGCTTTTAGCGTACGGGAAg                                                   <  3:1066915/100‑1 (MQ=255)
                                                                ccgcATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGTCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAg                                                   >  4:1066915/1‑100 (MQ=255)
                                                                ccgcATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGa                                                     <  2:903371/98‑1 (MQ=255)
                                                                ccgcATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGa                                                     >  1:903371/1‑98 (MQ=255)
                                                                ccgcATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCttttt                                 <  3:996202/118‑1 (MQ=255)
                                                                ccgcATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCttttt                                 >  4:996202/1‑118 (MQ=255)
                                                                 cggatgcGGCGTGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCTTTTTCCGCGCCATCTTTCATCACCTGTTGCAGGaa  <  1:806981/142‑1 (MQ=255)
                                                                 cgcATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCt                                         >  2:758188/1‑109 (MQ=255)
                                                                 cgcATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCt                                         <  1:758188/109‑1 (MQ=255)
                                                                           |                                                                                                                                         
CTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCTTTTTCCGCGCCATCTTTCATCACCTGTTGCAGGAA  >  NC_000913/3512526‑3512738

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