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.T88.5% 147.4 / 12.0 52intergenic (+173/+329)tyrB/yjbStyrosine aminotransferase, tyrosine‑repressible, PLP‑dependent/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base T (18/28);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

AAACTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTA  >  NC_000913/4268350‑4268595
                                                                                                                                   |                                                                                                                   
aaaCTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGTTGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑tttttttt                                                                                                             >  1:1273202/1‑138 (MQ=255)
       cgtcgCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGa                                                                                                       <  1:386714/139‑1 (MQ=255)
                    ctttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATCAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTTGAGTAATAATTTCTg                                                                                          <  1:1199206/139‑1 (MQ=255)
                    ctttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          <  2:1006657/139‑1 (MQ=255)
                    ctttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          <  2:139930/139‑1 (MQ=255)
                      ttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca                                                                                        <  2:391510/139‑1 (MQ=255)
                      ttctttATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca                                                                                        <  1:736648/139‑1 (MQ=255)
                           tataGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCt                                                                                   <  2:1273202/139‑1 (MQ=255)
                           tataGATGATTTTTTTGATGCCATAGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          >  1:1236751/1‑132 (MQ=255)
                           tataGATGATTTTTTTGATGCCATAGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          <  2:1236751/132‑1 (MQ=255)
                               gatgatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGa                                                                                                       >  1:1517484/1‑115 (MQ=255)
                               gatgatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGa                                                                                                       <  2:1517484/115‑1 (MQ=255)
                               gatgatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑TTTTTTTtggg                                                                                                          >  2:458404/1‑110 (MQ=255)
                               gatgatTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑TTTTTTTtggg                                                                                                          <  1:458404/111‑2 (MQ=255)
                                          tGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAAtt                                                                    <  1:1634007/139‑1 (MQ=255)
                                          tGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAAtt                                                                    <  2:1732420/139‑1 (MQ=255)
                                                cATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataa                                                              <  2:1548724/139‑1 (MQ=255)
                                                 aTCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          >  1:454804/1‑110 (MQ=255)
                                                 aTCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          <  2:454804/110‑1 (MQ=255)
                                                  tCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTtaataata                                                            >  2:1568104/1‑139 (MQ=255)
                                                       cTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          <  1:340212/104‑1 (MQ=255)
                                                       cTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          >  2:340212/1‑104 (MQ=255)
                                                              agagaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCTAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGTACATCTAACCATTATTAAATTTATTAATAATAATCGGGttt                                                >  1:239645/1‑139 (MQ=255)
                                                                agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          <  1:644027/95‑1 (MQ=255)
                                                                agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          >  2:644027/1‑95 (MQ=255)
                                                                agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGca                                                                                        <  2:1223697/97‑1 (MQ=255)
                                                                agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAAc                                                                                >  1:981497/1‑105 (MQ=255)
                                                                agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAAc                                                                                <  2:981497/105‑1 (MQ=255)
                                                                agaTAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCGGTTTTTTTTGGGTGGAGTAATAATTTCTGca                                                                                        >  1:1223697/1‑97 (MQ=255)
                                                                  ataataAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAAcc                                                                               >  1:1375421/1‑104 (MQ=255)
                                                                  ataataAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAAcc                                                                               <  2:1375421/104‑1 (MQ=255)
                                                                             tgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          >  2:1566298/1‑82 (MQ=255)
                                                                             tgttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTCATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          <  1:1566298/82‑1 (MQ=255)
                                                                              gttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATAtt                                >  1:1536484/1‑139 (MQ=255)
                                                                               ttAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATAttt                               >  1:77688/1‑139 (MQ=255)
                                                                                 aGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTAt                                                       <  1:767402/113‑1 (MQ=255)
                                                                                 aGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTAt                                                       >  2:767402/1‑113 (MQ=255)
                                                                                 aGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCt                             <  1:151454/139‑1 (MQ=255)
                                                                                    tCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCtata                          <  2:77688/139‑1 (MQ=255)
                                                                                              ttAAGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGaa                <  1:1568104/139‑1 (MQ=255)
                                                                                                aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          <  2:615265/63‑1 (MQ=255)
                                                                                                aaGCTTATCCCAATTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          >  1:615265/1‑63 (MQ=255)
                                                                                                        cccAATTATCTGGAATTCCTTATCCTG‑TTTTTTTtggg                                                                                                          >  1:191722/1‑37 (MQ=255)
                                                                                                        cccAATTATCTGGAATTCCTTATCCTG‑TTTTTTTtggg                                                                                                          <  2:191722/38‑2 (MQ=255)
                                                                                                            aTTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCatta                                                                           >  2:1335575/1‑66 (MQ=255)
                                                                                                            aTTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCatta                                                                           <  1:1335575/66‑1 (MQ=255)
                                                                                                            aTTATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTa  <  2:239645/139‑1 (MQ=255)
                                                                                                             ttATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCt                      <  2:1506927/118‑1 (MQ=255)
                                                                                                             ttATCTGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTAACGGGTTTGTTACGATAGAATATTTCTATATCCt                      >  1:1506927/1‑118 (MQ=255)
                                                                                                             ttATCTGGAATTCCTTATCCTG‑TTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          >  1:1579634/1‑49 (MQ=255)
                                                                                                             ttATCTGGAATTCCTTATCCTG‑TTTTTTTGGGTGGAGTAATAATTTCTg                                                                                          <  2:1579634/49‑1 (MQ=255)
                                                                                                                  tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                  >  1:1285367/1‑37 (MQ=255)
                                                                                                                  tGGAATTCCTTATCCTGTTTTTTTTGGGTGGAGtaat                                                                                                  <  2:1285367/37‑1 (MQ=255)
                                                                                                                                   |                                                                                                                   
AAACTGTCGTCGCTTTCACTCTTTCTTTATAGATGATTTTTTTGATGCCATCGTTCTACGTGAGAGATAATAAACGTTGTTAGTTCTTTTATTGTTAAGCTTATCCCAATTATCTGGAATTCCTTATCCTG‑TTTTTTTGGGTGGAGTAATAATTTCTGCACATCTAACCATTATTAAATTTAATAATAAGTATCGGGTTTGTTACGATAGAATATTTCTATATCCTCAAGAATATGCTTTGCCTTA  >  NC_000913/4268350‑4268595

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