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

GTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑TTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGA  >  NC_000913/4268355‑4268613
                                                                                                                              |                                                                                                                                      
gtcgtcGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑ttt‑tttttggg                                                                                                                            <  1:1429344/137‑5 (MQ=255)
gtcgtcGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑TTT‑TTTTtggg                                                                                                                            >  2:1429344/1‑136 (MQ=255)
 tcgtcgCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑ttt‑ttt                                                                                                                                 >  2:930184/1‑131 (MQ=255)
 tcgtcgCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑ttt‑ttt                                                                                                                                 <  1:930184/131‑1 (MQ=255)
       cTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGtaat                                                                                                                    <  2:1400915/139‑1 (MQ=255)
          tCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTaataat                                                                                                                 <  1:525152/139‑1 (MQ=255)
           cACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAAtt                                                                                                                >  2:853571/1‑139 (MQ=255)
               ctttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  2:100651/139‑1 (MQ=255)
                 ttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGca                                                                                                          <  2:13487/139‑1 (MQ=255)
                      tataGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑ttt‑tttttggg                                                                                                                            >  1:525764/1‑111 (MQ=255)
                      tataGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑TTT‑TTTTtggg                                                                                                                            <  2:525764/115‑2 (MQ=255)
                      tataAATTATTTTTTTGATGCCATCGTTCTACGTGAGAGATACTAAACGTTGTTAGTTCTTATATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCt                                                                                                     <  2:1188060/139‑1 (MQ=255)
                             gatTTTTTTGATGCCATCGTTCTACGTGAGAGGTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCTATTTATCTGGAAGTCCTTATCCTGTTTT‑TTTTGGCTGAAGTACTAATTCCTGCGCATCTAACAatt                                                                                              >  2:985183/1‑139 (MQ=255)
                              atTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCatta                                                                                             >  2:1346812/1‑139 (MQ=255)
                                ttttttGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑ttt‑tttttggg                                                                                                                            <  2:1381070/105‑5 (MQ=255)
                                ttttttGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑TTT‑TTTTtggg                                                                                                                            >  1:1381070/1‑104 (MQ=255)
                                            aTCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataat                                                                               <  1:711257/139‑1 (MQ=255)
                                                     cGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑TTTCTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCgg                                                                      >  2:359172/1‑139 (MQ=255)
                                                       tGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  2:1567436/1‑99 (MQ=255)
                                                       tGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  1:1567436/99‑1 (MQ=255)
                                                           agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAAc                                                                                                  <  2:815571/105‑1 (MQ=255)
                                                           agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAAc                                                                                                  >  1:815571/1‑105 (MQ=255)
                                                               aataaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTg                                                                                                                           <  1:1492627/76‑1 (MQ=255)
                                                               aataaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTg                                                                                                                           >  2:1492627/1‑76 (MQ=255)
                                                                     cGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGaa                                                      <  1:729537/139‑1 (MQ=255)
                                                                      gttgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  2:1278659/1‑84 (MQ=255)
                                                                      gttgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  1:1278659/84‑1 (MQ=255)
                                                                      gttgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAat                                                     >  2:885914/1‑139 (MQ=255)
                                                                       ttgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAata                                                    >  2:1195833/1‑139 (MQ=255)
                                                                         gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGca                                                                                                          >  1:1248661/1‑83 (MQ=255)
                                                                         gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGca                                                                                                          <  2:1248661/83‑1 (MQ=255)
                                                                         gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAAc                                                                                                  >  2:334201/1‑91 (MQ=255)
                                                                         gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAAc                                                                                                  <  1:334201/91‑1 (MQ=255)
                                                                         gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑TTT‑TTTTtggg                                                                                                                            <  1:1585799/64‑2 (MQ=255)
                                                                         gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑TTT‑TTTTtggg                                                                                                                            >  2:1585799/1‑63 (MQ=255)
                                                                                    tATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTc                                       >  2:257593/1‑139 (MQ=255)
                                                                                           aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  1:1166592/1‑63 (MQ=255)
                                                                                           aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  2:1166592/63‑1 (MQ=255)
                                                                                               ttATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCt                            <  1:885914/139‑1 (MQ=255)
                                                                                                tATCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCtt                           <  1:359172/139‑1 (MQ=255)
                                                                                                 aTCCCAATTCTCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTCTATCCTCAAGAATATGCttt                          <  2:974588/139‑1 (MQ=255)
                                                                                                  tCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGtaat                                                                                                                    >  2:84788/1‑48 (MQ=255)
                                                                                                  tCCCAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGtaat                                                                                                                    <  1:84788/48‑1 (MQ=255)
                                                                                                   cccAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  1:27231/1‑55 (MQ=255)
                                                                                                   cccAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  2:27231/55‑1 (MQ=255)
                                                                                                   cccAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCtat                                             <  2:1208843/118‑1 (MQ=255)
                                                                                                   cccAATTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCtat                                             >  1:1208843/1‑118 (MQ=255)
                                                                                                       aTTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  2:1343365/1‑51 (MQ=255)
                                                                                                       aTTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  1:1343365/51‑1 (MQ=255)
                                                                                                       aTTATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTa                    <  1:853571/139‑1 (MQ=255)
                                                                                                        ttATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAAtt                                                                                      <  2:883154/72‑1 (MQ=255)
                                                                                                        ttATCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAAtt                                                                                      >  1:883154/1‑72 (MQ=255)
                                                                                                           tCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  2:145061/47‑1 (MQ=255)
                                                                                                           tCTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  1:145061/1‑47 (MQ=255)
                                                                                                            cTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCattatt                                                                                           >  2:1576917/1‑63 (MQ=255)
                                                                                                            cTGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCattatt                                                                                           <  1:1576917/63‑1 (MQ=255)
                                                                                                             tGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCtata                                            <  1:1459713/109‑1 (MQ=255)
                                                                                                             tGGAATTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCtata                                            >  2:1459713/1‑109 (MQ=255)
                                                                                                                aaTTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCggg                                                                     <  2:543910/81‑1 (MQ=255)
                                                                                                                aaTTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCggg                                                                     >  1:543910/1‑81 (MQ=255)
                                                                                                                aaTTCCTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCtgtg           >  1:954927/1‑139 (MQ=255)
                                                                                                                     cTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCt                                               >  2:803309/1‑98 (MQ=255)
                                                                                                                     cTTATCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCt                                               <  1:803309/98‑1 (MQ=255)
                                                                                                                         tCCTGTTTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTga  <  1:977254/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                      
GTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑TTT‑TTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGA  >  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: ‑