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.T97.2% 112.0 / ‑3.4 36intergenic (+173/+329)tyrB/yjbStyrosine aminotransferase, tyrosine‑repressible, PLP‑dependent/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (18/17);  total (19/17)
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.86e-01

TGAAATTAAACTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGA  >  NC_000913/4268343‑4268613
                                                                                                                                          |                                                                                                                                     
tGAAATTAAACTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGt                                                                                                                                      <  2:137080/139‑1 (MQ=255)
          cTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGttttttttggg                                                                                                                            <  2:598022/139‑5 (MQ=255)
               gtcgCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAg                                                                                                                        >  2:150221/1‑139 (MQ=255)
                       cACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                                    >  1:941709/1‑135 (MQ=255)
                       cACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                                    <  2:941709/135‑1 (MQ=255)
                             ttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca                                                                                                          <  2:398629/139‑1 (MQ=255)
                                  tataGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca                                                                                                          >  1:1205662/1‑134 (MQ=255)
                                  tataGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca                                                                                                          <  2:1205662/134‑1 (MQ=255)
                                    taGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGttttttttggg                                                                                                                            <  1:1192697/113‑5 (MQ=255)
                                    taGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGttttttttggg                                                                                                                            >  2:1192697/1‑109 (MQ=255)
                                         gatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  2:376284/125‑1 (MQ=255)
                                         gatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  1:376284/1‑125 (MQ=255)
                                         gatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCatt                                                                                              >  1:377272/1‑139 (MQ=255)
                                                       cATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataa                                                                                <  1:152514/139‑1 (MQ=255)
                                                                       agaTAATAAACGTTGTTAGTTCTTTTATTGTTATGCTTAGCCCAATTATCTGGAATTCCTTATCCTGttttttttgg                                                                                                                             <  2:914606/77‑4 (MQ=255)
                                                                       agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGc                                                                                                           <  1:382195/96‑1 (MQ=255)
                                                                       agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGc                                                                                                           >  2:382195/1‑96 (MQ=255)
                                                                       agaTAATAAACGTTGCTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTtgg                                                                                                                             >  1:914606/1‑77 (MQ=255)
                                                                           aataaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTTTCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGt                                                                <  2:459389/135‑1 (MQ=255)
                                                                           aataaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTTTCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGt                                                                >  1:459389/1‑135 (MQ=255)
                                                                           aataaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaat                                                                                  <  1:921895/117‑1 (MQ=255)
                                                                           aataaACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaat                                                                                  >  2:921895/1‑117 (MQ=255)
                                                                                 cGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                                    >  1:827049/1‑77 (MQ=255)
                                                                                 cGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                                    <  2:827049/77‑1 (MQ=255)
                                                                                            cTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCa                                                                                                >  1:827778/1‑86 (MQ=255)
                                                                                            cTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCa                                                                                                <  2:827778/86‑1 (MQ=255)
                                                                                              tttATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATcc                                         >  1:576255/1‑139 (MQ=255)
                                                                                               ttATTGTTCAGTTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCt                                        <  1:831385/139‑1 (MQ=255)
                                                                                               ttATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCt                                        <  2:576255/139‑1 (MQ=255)
                                                                                                       aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAata                                >  1:503997/1‑139 (MQ=255)
                                                                                                             aTCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataa                                                                                >  2:418446/1‑85 (MQ=255)
                                                                                                             aTCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataa                                                                                <  1:418446/85‑1 (MQ=255)
                                                                                                                   aTTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  2:313942/51‑1 (MQ=255)
                                                                                                                   aTTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  1:313942/1‑51 (MQ=255)
                                                                                                                    ttATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTaataat                                                                                                                 >  2:1329087/1‑45 (MQ=255)
                                                                                                                    ttATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTaataat                                                                                                                 <  1:1329087/45‑1 (MQ=255)
                                                                                                                         tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCt                                                                                                     >  2:301496/1‑52 (MQ=39)
                                                                                                                         tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCt                                                                                                     <  1:301496/52‑1 (MQ=39)
                                                                                                                                 cTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            <  2:177035/37‑1 (MQ=255)
                                                                                                                                 cTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                            >  1:177035/1‑37 (MQ=255)
                                                                                                                                     tCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTga  >  1:197515/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                     
TGAAATTAAACTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGA  >  NC_000913/4268343‑4268613

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