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

CGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGAC  >  NC_000913/4268357‑4268606
                                                                                                                            |                                                                                                                               
cgtcgCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGa                                                                                                                  >  2:1708755/1‑139 (MQ=255)
 gtcgCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAg                                                                                                                 >  1:1365025/1‑139 (MQ=255)
     cTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                             <  2:356716/139‑1 (MQ=255)
         cACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAAtt                                                                                                         >  2:1471882/1‑139 (MQ=255)
             ctttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                     <  2:774173/139‑1 (MQ=255)
             ctttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                     <  2:1365025/139‑1 (MQ=255)
             ctttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTACTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                     <  1:1471882/139‑1 (MQ=255)
             ctttctttATAGATGATTTTTTTAATGCCATCCTTCTACGTCGGAAATTATAAATGATGTTAGTTCTTTTATTGTTAAGCTTATCACAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                     <  1:1445407/139‑1 (MQ=255)
                tctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                     <  2:1658615/136‑1 (MQ=255)
                tctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                     >  1:1658615/1‑136 (MQ=255)
                tctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑ttttttttggg                                                                                                                     >  1:1001348/1‑115 (MQ=255)
                tctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑ttttttttggg                                                                                                                     <  2:1001348/119‑5 (MQ=255)
                       aGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTaa                                                                                            <  1:1621760/139‑1 (MQ=255)
                           gatTTTTTTGATGCCATTGTTATACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑TTTTTTTtggg                                                                                                                     <  1:863700/108‑2 (MQ=255)
                           gatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑TTTTTTTtggg                                                                                                                     >  2:863700/1‑107 (MQ=255)
                            atTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCatta                                                                                      >  1:1949799/1‑139 (MQ=255)
                              ttttttGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCattatt                                                                                    <  1:1708755/139‑1 (MQ=255)
                                          aTCGTTCTACGTGAGAGATTATAAATGTTGTTCGTTCTTTTATTGTTAAGCTTATCACCATTATCTGGACTTCCTTCCCCTGTTTTTTTTGGATGGAGTAATAATTTCTGCACATCTACCCATTATTAAATTtaataat                                                                        <  1:108047/139‑1 (MQ=255)
                                                       agagaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGttt                                                           >  2:1657523/1‑139 (MQ=255)
                                                        gagaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTT‑‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGt                                                         <  2:1452344/139‑1 (MQ=255)
                                                         agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca                                                                                                   <  1:1084980/97‑1 (MQ=255)
                                                         agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca                                                                                                   >  2:1084980/1‑97 (MQ=255)
                                                                     ttgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑ttttttttgggt                                                                                                                    >  1:589182/1‑62 (MQ=255)
                                                                     ttgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑TTTTTTTtgggt                                                                                                                    <  2:589182/67‑3 (MQ=255)
                                                                      tgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑ttttttttg                                                                                                                       >  1:2024381/1‑61 (MQ=255)
                                                                      tgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑ttttttttg                                                                                                                       <  2:2024381/63‑3 (MQ=255)
                                                                       gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGaa                                               <  2:1577945/135‑1 (MQ=255)
                                                                       gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGaa                                               >  1:1577945/1‑135 (MQ=255)
                                                                       gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATAtt                                           >  2:820627/1‑139 (MQ=255)
                                                                          aGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCt                                        <  2:1949799/139‑1 (MQ=255)
                                                                                tttATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATcc                                  >  1:725256/1‑139 (MQ=255)
                                                                                               aTCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTc                                >  2:1468061/1‑126 (MQ=255)
                                                                                               aTCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTc                                <  1:1468061/126‑1 (MQ=255)
                                                                                               aTCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCttt                   >  1:483855/1‑139 (MQ=255)
                                                                                                tCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                             >  1:173808/1‑48 (MQ=255)
                                                                                                tCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                             <  2:173808/48‑1 (MQ=255)
                                                                                                     aTTATCTGGAATTCCTTATTTTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTGCTATATCCTCAAGAATATGCTTTGCCTTa             <  1:820627/139‑1 (MQ=255)
                                                                                                      ttATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAatat                        >  1:29877/1‑127 (MQ=255)
                                                                                                      ttATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAatat                        <  2:29877/127‑1 (MQ=255)
                                                                                                      ttATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTAt            >  1:1816215/1‑139 (MQ=255)
                                                                                                           tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                     >  2:638483/1‑45 (MQ=255)
                                                                                                           tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                                     <  1:638483/45‑1 (MQ=255)
                                                                                                             gAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCtgt     <  2:1816215/139‑1 (MQ=255)
                                                                                                                ttCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGAc  >  1:1042888/1‑139 (MQ=255)
                                                                                                                            |                                                                                                                               
CGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGAC  >  NC_000913/4268357‑4268606

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