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.T86.8% 107.7 / 8.8 38intergenic (+173/+329)tyrB/yjbStyrosine aminotransferase, tyrosine‑repressible, PLP‑dependent/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base T (17/16);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GTGAAATTAAACTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTG  >  NC_000913/4268342‑4268602
                                                                                                                                           |                                                                                                                           
gTGAAATTAAACTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTg                                                                                                                              <  2:1017997/139‑1 (MQ=255)
           cTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTtggg                                                                                                                 <  2:1159177/139‑2 (MQ=255)
                            ctttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                 <  1:776910/139‑1 (MQ=255)
                                ctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                         <  1:284609/127‑1 (MQ=255)
                                ctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                         >  2:284609/1‑127 (MQ=255)
                                      aGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAAc                                                                                       <  1:47090/139‑1 (MQ=255)
                                        atgatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGa                                                                                                              >  2:160486/1‑114 (MQ=255)
                                        atgatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGa                                                                                                              <  1:160486/114‑1 (MQ=255)
                                        atgatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                         >  2:1639724/1‑119 (MQ=255)
                                        atgatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                         <  1:1639724/119‑1 (MQ=255)
                                           atTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCatta                                                                                  >  1:1870071/1‑139 (MQ=255)
                                           atTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGttttttttggg                                                                                                                 >  2:510310/1‑103 (MQ=255)
                                           atTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGttttttttggg                                                                                                                 <  1:510310/107‑5 (MQ=255)
                                             ttttttGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCattatt                                                                                >  2:1113854/1‑139 (MQ=255)
                                                    aTGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTta                                                                         <  1:1113854/139‑1 (MQ=255)
                                                         aTCGTTCTACGTGAGAGATAATAAACGTTGTTATTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataat                                                                    <  2:1870071/139‑1 (MQ=255)
                                                                       gagaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                 >  1:1641204/1‑96 (MQ=255)
                                                                       gagaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                 <  2:1641204/96‑1 (MQ=255)
                                                                        agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGt                                                     >  1:294611/1‑139 (MQ=255)
                                                                            aataaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACg                                                 >  2:669160/1‑139 (MQ=255)
                                                                                  cGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTaataat                                                                                                      >  2:1043265/1‑80 (MQ=255)
                                                                                  cGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTCTTTGGGTGGAGTaataat                                                                                                      <  1:1043265/80‑1 (MQ=255)
                                                                                     tgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGttttttttggg                                                                                                                 <  1:347376/65‑5 (MQ=255)
                                                                                     tgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGttttttttggg                                                                                                                 >  2:347376/1‑61 (MQ=255)
                                                                                      gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGttttttttggg                                                                                                                 >  1:445170/1‑60 (MQ=255)
                                                                                      gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGttttttttggg                                                                                                                 <  2:445170/64‑5 (MQ=255)
                                                                                                 tATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTc                            >  1:1620834/1‑139 (MQ=255)
                                                                                                     gTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGa                        >  1:889035/1‑139 (MQ=255)
                                                                                                      ttAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGaa                       <  2:889035/139‑1 (MQ=255)
                                                                                                      ttAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGaa                       <  2:294611/139‑1 (MQ=255)
                                                                                                      ttAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGaa                       <  2:650256/139‑1 (MQ=255)
                                                                                                      ttAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTtggg                                                                                                                 >  2:1351987/1‑47 (MQ=255)
                                                                                                      ttAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTtggg                                                                                                                 <  1:1351987/48‑2 (MQ=255)
                                                                                                        aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCattatt                                                                                >  2:163588/1‑80 (MQ=255)
                                                                                                        aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCattatt                                                                                <  1:163588/80‑1 (MQ=255)
                                                                                                        aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGttttttttggg                                                                                                                 >  2:187102/1‑42 (MQ=255)
                                                                                                        aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTtggg                                                                                                                 <  1:187102/46‑2 (MQ=255)
                                                                                                              aTCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGa                                                <  2:203121/106‑1 (MQ=255)
                                                                                                              aTCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGa                                                >  1:203121/1‑106 (MQ=255)
                                                                                                                    aTTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAAtttatg                                                                                                 >  2:1563927/1‑48 (MQ=255)
                                                                                                                    aTTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAAtttatg                                                                                                 <  1:1563927/51‑4 (MQ=255)
                                                                                                                     ttATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGt                                                     <  2:1670103/94‑1 (MQ=255)
                                                                                                                     ttATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGt                                                     >  1:1670103/1‑94 (MQ=255)
                                                                                                                        tCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                 <  2:1099560/47‑1 (MQ=255)
                                                                                                                        tCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                 >  1:1099560/1‑47 (MQ=255)
                                                                                                                           ggAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCtg  >  2:1819409/1‑139 (MQ=255)
                                                                                                                                           |                                                                                                                           
GTGAAATTAAACTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTG  >  NC_000913/4268342‑4268602

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