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.T93.3% 88.1 / ‑0.3 30intergenic (+173/+329)tyrB/yjbStyrosine aminotransferase, tyrosine‑repressible, PLP‑dependent/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base T (13/15);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTAAACTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGA  >  NC_000913/4268346‑4268613
                                                                                                                                       |                                                                                                                                      
aaTTAAACTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑tttt                                                                                                                                   >  2:452086/1‑139 (MQ=255)
       cTGTCGTCGCTTTCTCTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACCTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑ttttttttggg                                                                                                                            >  1:933002/1‑135 (MQ=255)
       cTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑TTTTTTTtggg                                                                                                                            <  1:1194538/139‑2 (MQ=255)
       cTGTCGTCCCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑ttttttttggg                                                                                                                            <  2:933002/139‑5 (MQ=255)
           cgtcgCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGa                                                                                                                         <  1:340673/139‑1 (MQ=255)
                  ttCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGtaataa                                                                                                                  >  1:924597/1‑139 (MQ=255)
                         tttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCtac                                                                                                           <  1:452086/139‑3 (MQ=255)
                               tataGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCt                                                                                                     <  2:213044/139‑1 (MQ=255)
                               tataGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCt                                                                                                     <  2:1392571/139‑1 (MQ=255)
                                     tgatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCat                                                                                               >  1:346050/1‑139 (MQ=255)
                                       atTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑ttttttttggg                                                                                                                            >  2:317068/1‑103 (MQ=255)
                                       atTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑TTTTTTTtggg                                                                                                                            <  1:317068/107‑2 (MQ=255)
                                              tGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAAt                                                                                       <  2:472138/139‑1 (MQ=255)
                                                              cGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCgg                                                                      >  2:110861/1‑139 (MQ=255)
                                                               gTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCggg                                                                     <  2:346050/139‑1 (MQ=255)
                                                                             aCGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTCTCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGa                                                       <  2:521779/139‑1 (MQ=255)
                                                                                  gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGc                                                                                                           >  1:163325/1‑82 (MQ=255)
                                                                                  gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGc                                                                                                           <  2:163325/82‑1 (MQ=255)
                                                                                           tttATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAata                                                    <  1:875237/128‑1 (MQ=255)
                                                                                           tttATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAata                                                    >  2:875237/1‑128 (MQ=255)
                                                                                            ttATTTTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCt                                        <  1:1459878/139‑1 (MQ=255)
                                                                                            ttATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCt                                        >  1:1462505/1‑139 (MQ=255)
                                                                                                         tATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  2:1168938/1‑58 (MQ=255)
                                                                                                         tATCCCAATTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  1:1168938/58‑1 (MQ=255)
                                                                                                               aaTTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataat                                                                               <  2:183633/81‑1 (MQ=255)
                                                                                                               aaTTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataat                                                                               >  1:183633/1‑81 (MQ=255)
                                                                                                                aTTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  2:437351/51‑1 (MQ=255)
                                                                                                                aTTATCTGGAATTCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  1:437351/1‑51 (MQ=255)
                                                                                                                           ttCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGt                                                                >  1:473401/1‑84 (MQ=255)
                                                                                                                           ttCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGt                                                                <  2:473401/84‑1 (MQ=255)
                                                                                                                           ttCCTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGAc         >  2:1114266/1‑139 (MQ=255)
                                                                                                                              cTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  2:156354/37‑1 (MQ=255)
                                                                                                                              cTTATCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  1:156354/1‑37 (MQ=255)
                                                                                                                                  tCCTGT‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTga  >  1:1320526/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                      
AATTAAACTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGA  >  NC_000913/4268346‑4268613

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