Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,268,487 (T)7→8 intergenic (+180/+322) tyrB → / ← yjbS tyrosine aminotransferase, tyrosine‑repressible, PLP‑dependent/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,268,4801.T89.8% 172.4 / 11.5 59intergenic (+173/+329)tyrB/yjbStyrosine aminotransferase, tyrosine‑repressible, PLP‑dependent/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base T (24/29);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

GTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGA  >  NC_000913/4268355‑4268613
                                                                                                                              |                                                                                                                                      
gtcgtcGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGttttttttggg                                                                                                                            <  1:1429344/137‑5 (MQ=255)
gtcgtcGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTtggg                                                                                                                            >  2:1429344/1‑136 (MQ=255)
 tcgtcgCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGtttttt                                                                                                                                 >  2:930184/1‑131 (MQ=255)
 tcgtcgCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGtttttt                                                                                                                                 <  1:930184/131‑1 (MQ=255)
       cTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                                    <  2:1400915/139‑1 (MQ=255)
          tCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTaataat                                                                                                                 <  1:525152/139‑1 (MQ=255)
           cACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAAtt                                                                                                                >  2:853571/1‑139 (MQ=255)
               ctttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  2:100651/139‑1 (MQ=255)
                 ttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca                                                                                                          <  2:13487/139‑1 (MQ=255)
                      tataGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGttttttttggg                                                                                                                            >  1:525764/1‑111 (MQ=255)
                      tataGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTtggg                                                                                                                            <  2:525764/115‑2 (MQ=255)
                      tataAATTATTTTTTTGATGCCATCGTTCTACGTGAGAGATACTAAACGTTGTTAGTTCTTATATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCt                                                                                                     <  2:1188060/139‑1 (MQ=255)
                             gatTTTTTTGATGCCATCGTTCTACGTGAGAGGTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCTATTTATCTGGAAGTCCTTATCCTGTTTTTTTTGGCTGAAGTACTAATTCCTGCGCATCTAACAatt                                                                                              >  2:985183/1‑139 (MQ=255)
                              atTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCatta                                                                                             >  2:1346812/1‑139 (MQ=255)
                                ttttttGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGttttttttggg                                                                                                                            <  2:1381070/105‑5 (MQ=255)
                                ttttttGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTtggg                                                                                                                            >  1:1381070/1‑104 (MQ=255)
                                            aTCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataat                                                                               <  1:711257/139‑1 (MQ=255)
                                                     cGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTCTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCgg                                                                      >  2:359172/1‑139 (MQ=255)
                                                       tGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  2:1567436/1‑99 (MQ=255)
                                                       tGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  1:1567436/99‑1 (MQ=255)
                                                           agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAAc                                                                                                  <  2:815571/105‑1 (MQ=255)
                                                           agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAAc                                                                                                  >  1:815571/1‑105 (MQ=255)
                                                               aataaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTg                                                                                                                           <  1:1492627/76‑1 (MQ=255)
                                                               aataaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTg                                                                                                                           >  2:1492627/1‑76 (MQ=255)
                                                                     cGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGaa                                                      <  1:729537/139‑1 (MQ=255)
                                                                      gttgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  2:1278659/1‑84 (MQ=255)
                                                                      gttgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  1:1278659/84‑1 (MQ=255)
                                                                      gttgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAat                                                     >  2:885914/1‑139 (MQ=255)
                                                                       ttgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAata                                                    >  2:1195833/1‑139 (MQ=255)
                                                                         gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca                                                                                                          >  1:1248661/1‑83 (MQ=255)
                                                                         gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca                                                                                                          <  2:1248661/83‑1 (MQ=255)
                                                                         gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAAc                                                                                                  >  2:334201/1‑91 (MQ=255)
                                                                         gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAAc                                                                                                  <  1:334201/91‑1 (MQ=255)
                                                                         gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTtggg                                                                                                                            <  1:1585799/64‑2 (MQ=255)
                                                                         gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTtggg                                                                                                                            >  2:1585799/1‑63 (MQ=255)
                                                                                    tATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTc                                       >  2:257593/1‑139 (MQ=255)
                                                                                           aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  1:1166592/1‑63 (MQ=255)
                                                                                           aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  2:1166592/63‑1 (MQ=255)
                                                                                               ttATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCt                            <  1:885914/139‑1 (MQ=255)
                                                                                                tATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCtt                           <  1:359172/139‑1 (MQ=255)
                                                                                                 aTCCCAATTCTCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTCTATCCTCAAGAATATGCttt                          <  2:974588/139‑1 (MQ=255)
                                                                                                  tCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                                    >  2:84788/1‑48 (MQ=255)
                                                                                                  tCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                                    <  1:84788/48‑1 (MQ=255)
                                                                                                   cccAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  1:27231/1‑55 (MQ=255)
                                                                                                   cccAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  2:27231/55‑1 (MQ=255)
                                                                                                   cccAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCtat                                             <  2:1208843/118‑1 (MQ=255)
                                                                                                   cccAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCtat                                             >  1:1208843/1‑118 (MQ=255)
                                                                                                       aTTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  2:1343365/1‑51 (MQ=255)
                                                                                                       aTTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  1:1343365/51‑1 (MQ=255)
                                                                                                       aTTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTa                    <  1:853571/139‑1 (MQ=255)
                                                                                                        ttATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAAtt                                                                                      <  2:883154/72‑1 (MQ=255)
                                                                                                        ttATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAAtt                                                                                      >  1:883154/1‑72 (MQ=255)
                                                                                                           tCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  2:145061/47‑1 (MQ=255)
                                                                                                           tCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  1:145061/1‑47 (MQ=255)
                                                                                                            cTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCattatt                                                                                           >  2:1576917/1‑63 (MQ=255)
                                                                                                            cTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCattatt                                                                                           <  1:1576917/63‑1 (MQ=255)
                                                                                                             tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCtata                                            <  1:1459713/109‑1 (MQ=255)
                                                                                                             tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCtata                                            >  2:1459713/1‑109 (MQ=255)
                                                                                                                aaTTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCggg                                                                     <  2:543910/81‑1 (MQ=255)
                                                                                                                aaTTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCggg                                                                     >  1:543910/1‑81 (MQ=255)
                                                                                                                aaTTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCtgtg           >  1:954927/1‑139 (MQ=255)
                                                                                                                     cTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCt                                               >  2:803309/1‑98 (MQ=255)
                                                                                                                     cTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCt                                               <  1:803309/98‑1 (MQ=255)
                                                                                                                         tCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTga  <  1:977254/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                      
GTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGA  >  NC_000913/4268355‑4268613

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