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.T95.1% 187.4 / 2.0 61intergenic (+173/+329)tyrB/yjbStyrosine aminotransferase, tyrosine‑repressible, PLP‑dependent/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base T (28/30);  total (29/32)
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.20e-01

CTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTT  >  NC_000913/4268353‑4268609
                                                                                                                                |                                                                                                                                  
cTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑ttttttttggg                                                                                                                        <  2:864881/139‑5 (MQ=255)
cTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑TTTTTTTtggg                                                                                                                        <  1:371662/139‑2 (MQ=255)
   tcgtcgCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTgg                                                                                                                      >  1:865502/1‑139 (MQ=255)
    cgtcgCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑tttttttttgggtgg                                                                                                                    >  2:1458922/1‑131 (MQ=255)
        gCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaa                                                                                                                 >  1:941419/1‑139 (MQ=255)
            tCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTaataat                                                                                                             <  2:1460237/139‑1 (MQ=255)
                 ctttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                        <  2:865502/139‑1 (MQ=255)
                 ctttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                        <  1:1538332/139‑1 (MQ=255)
                 ctttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                        <  1:1078705/139‑1 (MQ=255)
                  tttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGc                                                                                                       >  1:330753/1‑139 (MQ=255)
                          taGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                        >  1:1375839/1‑130 (MQ=255)
                          taGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                        <  2:1375839/130‑1 (MQ=255)
                           aGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAAc                                                                                              <  2:1518282/139‑1 (MQ=255)
                            gatgatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑ttttttttggg                                                                                                                        >  1:1472652/1‑107 (MQ=255)
                            gatgatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑ttttttttggg                                                                                                                        <  2:1472652/111‑5 (MQ=255)
                               gatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑TTTTTTTtggg                                                                                                                        >  1:1382681/1‑107 (MQ=255)
                               gatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑TTTTTTTtggg                                                                                                                        <  2:1382681/108‑2 (MQ=255)
                                       tGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCCGTTTTTTTTGGGTg                                                                                                                       >  1:579455/1‑102 (MQ=255)
                                       tGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCCGTTTTTTTTGGGTg                                                                                                                       <  2:579455/102‑1 (MQ=255)
                                             cATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataa                                                                            <  2:330753/139‑1 (MQ=255)
                                                           agagaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                        <  2:30821/97‑1 (MQ=255)
                                                           agagaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                        >  1:30821/1‑97 (MQ=255)
                                                           agagaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCggg                                                                 <  2:129519/136‑1 (MQ=255)
                                                           agagaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCggg                                                                 >  1:129519/1‑136 (MQ=255)
                                                            gagaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCATTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATGAGTATCGGGTTTg                                                             >  1:1063560/1‑139 (MQ=255)
                                                             agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                        >  2:703228/1‑95 (MQ=255)
                                                             agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                        <  1:703228/95‑1 (MQ=255)
                                                                  ataaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGa                                                       <  2:1063560/139‑1 (MQ=255)
                                                                   taaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCTCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGAt                                                      <  2:53333/139‑1 (MQ=255)
                                                                   taaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGa                                                       <  1:1458922/139‑1 (MQ=255)
                                                                          tgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAatat                                               >  1:217791/1‑139 (MQ=255)
                                                                           gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATAtt                                              >  1:608727/1‑139 (MQ=255)
                                                                              aGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGt                                                            <  2:350805/122‑1 (MQ=255)
                                                                              aGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGt                                                            >  1:350805/1‑122 (MQ=255)
                                                                                  cTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGc                                                                                                       <  1:231087/75‑1 (MQ=255)
                                                                                  cTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGc                                                                                                       >  2:231087/1‑75 (MQ=255)
                                                                                    tttATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGc                                                                                                       <  1:106545/73‑1 (MQ=255)
                                                                                    tttATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGc                                                                                                       >  2:106545/1‑73 (MQ=255)
                                                                                    tttATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCattatt                                                                                       <  2:131728/89‑1 (MQ=255)
                                                                                    tttATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCattatt                                                                                       >  1:131728/1‑89 (MQ=255)
                                                                                             aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca                                                                                                      >  1:58068/1‑65 (MQ=255)
                                                                                             aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca                                                                                                      <  2:58068/65‑1 (MQ=255)
                                                                                             aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCt                                    >  1:821239/1‑131 (MQ=255)
                                                                                             aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCt                                    <  2:821239/131‑1 (MQ=255)
                                                                                                cTTATCCCAATTATCTTGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAAt                                                                                   <  1:1184616/81‑1 (MQ=255)
                                                                                                cTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAAt                                                                                   >  2:1184616/1‑81 (MQ=255)
                                                                                                   aTCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataat                                                                           <  1:245504/86‑1 (MQ=255)
                                                                                                   aTCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataat                                                                           >  2:245504/1‑86 (MQ=255)
                                                                                                   aTCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataat                                                                           >  1:476151/1‑86 (MQ=255)
                                                                                                   aTCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataat                                                                           <  2:476151/86‑1 (MQ=255)
                                                                                                               tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                        >  2:55403/1‑45 (MQ=255)
                                                                                                               tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                        >  1:1550332/1‑45 (MQ=255)
                                                                                                               tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                        <  2:1550332/45‑1 (MQ=255)
                                                                                                               tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                        <  1:55403/45‑1 (MQ=255)
                                                                                                               tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGaa                                                  >  2:659151/1‑99 (MQ=255)
                                                                                                               tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGaa                                                  <  1:659151/99‑1 (MQ=255)
                                                                                                                ggAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataat                                                                           <  2:436204/73‑1 (MQ=255)
                                                                                                                ggAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataat                                                                           >  1:436204/1‑73 (MQ=255)
                                                                                                                  aaTTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCtata                                        <  2:742432/106‑1 (MQ=255)
                                                                                                                  aaTTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCtata                                        >  1:742432/1‑106 (MQ=255)
                                                                                                                       cTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                        >  2:529036/1‑37 (MQ=255)
                                                                                                                       cTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                        <  1:529036/37‑1 (MQ=255)
                                                                                                                       cTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGt                                                            <  2:1203642/81‑1 (MQ=255)
                                                                                                                       cTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGt                                                            >  1:1203642/1‑81 (MQ=255)
                                                                                                                       cTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTATGCCTTATTGCCTGTGACttt  >  2:85807/1‑139 (MQ=255)
                                                                                                                                |                                                                                                                                  
CTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTT  >  NC_000913/4268353‑4268609

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: