Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,890,154 C→T 9.2% intergenic (+10/‑81) tnaA → / → tnaB tryptophanase/L‑cysteine desulfhydrase, PLP‑dependent/tryptophan transporter of low affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,890,1540CT 9.2% 136.5 / 3.7 53intergenic (+10/‑81)tnaA/tnaBtryptophanase/L‑cysteine desulfhydrase, PLP‑dependent/tryptophan transporter of low affinity
Reads supporting (aligned to +/- strand):  ref base C (31/17);  new base T (3/2);  total (34/19)
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.78e-01

TCAAACACATATGGACTTCATTATTGAAGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCACTCTGCATTTTGGGGTGTTATGGTTATAGCAGGTACA  >  NC_000913/3890012‑3890297
                                                                                                                                              |                                                                                                                                               
tCAAACACATATGGACTTCATTATTGAAGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTg                                                                                                                                           >  1:270962/1‑149 (MQ=255)
tCAAACACATATGGACTTCATTATTGAAGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAAAGTg                                                                                                                                           >  1:461712/1‑149 (MQ=255)
  aaaCACATATGGACTTCATTATTGAAGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTAACTTCACCGCCAAACTTAAAGAAGTTTAATTAATACTACAAAgggga                                                                                                                                         >  2:423924/1‑148 (MQ=255)
              aCTTCATTATTGAAGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTATAGAGTGGCTATAAGGATGtt                                                                                                                             >  2:739666/1‑149 (MQ=255)
              aCTTCATTATTGAAGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTATAGAGTGGCTATAAGGATGtt                                                                                                                             >  2:732072/1‑149 (MQ=255)
              aCTTCATTATTGAAGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGAGTGtt                                                                                                                             >  1:393174/1‑149 (MQ=255)
               cTTCATTATTGAAGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGttc                                                                                                                            <  1:545077/149‑2 (MQ=255)
                ttCATTATTGAAGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTATAGAGTGGCTATAAGGATGTTAg                                                                                                                           >  1:520822/1‑149 (MQ=255)
                ttCATTATTGAAGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGCCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGAAGGTAg                                                                                                                           >  2:280379/1‑149 (MQ=255)
                  cATTATTGAAGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGcc                                                                                                                         >  2:770077/1‑149 (MQ=255)
                   attattGAAGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGATATAAGGATGTTAGCCa                                                                                                                        >  2:217067/1‑149 (MQ=255)
                    ttattGAAGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCAc                                                                                                                       <  1:747645/149‑1 (MQ=255)
                    ttattGAAGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACCTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCAc                                                                                                                       <  1:27358/149‑1 (MQ=255)
                          aaGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAATGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCtt                                                                                                                 >  1:471789/1‑149 (MQ=255)
                          aaGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCtt                                                                                                                 >  1:381503/1‑149 (MQ=255)
                          aaGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCtt                                                                                                                 >  1:620447/1‑149 (MQ=255)
                           aGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGCGGCTATAAGGATGTTAGCCACTCTCTTa                                                                                                                <  1:440921/149‑1 (MQ=255)
                               tttAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTCCCCt                                                                                                            >  1:81908/1‑149 (MQ=255)
                                  aaaCATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCACTTACCCTaaa                                                                                                         >  1:66159/1‑147 (MQ=255)
                                  aaaCATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGGCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATCCTACAGAGTGGCTCTACGGATGTTAGCCACCTTCTTACCCTACa                                                                                                         >  1:175185/1‑149 (MQ=255)
                                   aaCATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGCACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACAt                                                                                                        <  1:607249/149‑1 (MQ=255)
                                    aCATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGAGGCTATAAGGATATTAGCCACTCGCTTACCCTACATc                                                                                                       >  2:220636/1‑149 (MQ=255)
                                          gAAAGATAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAAt                                                                                                 <  1:169736/149‑1 (MQ=255)
                                           aaaGAGAACGCGGCGAATATTAAAGGATTAACCATTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCGTCAATa                                                                                                >  2:617539/1‑149 (MQ=255)
                                              gagaACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGCAGTTTAATTACTACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATaaaa                                                                                             >  2:76380/1‑147 (MQ=255)
                                                        cgAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCt                                                                                   >  1:621774/1‑149 (MQ=255)
                                                          aaTATTAAAGGGTTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTc                                                                                 >  2:482586/1‑149 (MQ=255)
                                                             atTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTAGCCTACATCCTCAATAACAAAAATACCCTTCctc                                                                              >  1:52293/1‑149 (MQ=255)
                                                             atTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCctc                                                                              >  1:642646/1‑149 (MQ=255)
                                                                    gagtAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGt                                                                       >  2:11096/4‑149 (MQ=255)
                                                                    gATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGGTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAGGGt                                                                       >  1:620742/1‑149 (MQ=255)
                                                                      ttAACCTTAACGAACAAACAGAAAGAATTGCGTCACTTCACCGGAAAACTTAAAGAAGTTCAATTAATACTACAGAGTGACTATAAAGATTTTAGCCACCCTCTTACCCTATAAAAAAAAAAAAAAAAATAGCCTTCCTCTAAAGGTgg                                                                     <  2:17801/149‑1 (MQ=255)
                                                                         aCCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcat                                                                  >  2:704604/1‑149 (MQ=255)
                                                                         aCCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcat                                                                  >  2:378816/1‑149 (MQ=255)
                                                                           cTTTACGTACGAACCGAAAGTATCGCGTAACTTCACCGCAAAACTTAAAGAAGATTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGcatca                                                                >  1:431438/1‑149 (MQ=255)
                                                                              tacgtacgACCCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGGGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGa                                                             <  1:619235/149‑1 (MQ=255)
                                                                                  tccgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGa                                                         <  2:277058/147‑1 (MQ=255)
                                                                                  tacgAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATCCTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCGCAATAACAAAAAGAGCCTTCCTCTAAAGGTGGCATCATGACTGa                                                         <  2:740670/149‑1 (MQ=255)
                                                                                     gAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCa                                                      <  1:144007/149‑1 (MQ=255)
                                                                                      aaCCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCaa                                                     >  2:690234/1‑149 (MQ=255)
                                                                                              gtgttGCGTCGCTTCACCGGCAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCGACATCCTCAATAACGAAAATAGCATTTCTCTGAAGGTGGCATCAGGACTGATCAAGCTGaaaa                                             <  2:134300/146‑1 (MQ=255)
                                                                                              gTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTATAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGaaaa                                             <  1:739666/149‑1 (MQ=255)
                                                                                              gTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTATAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGaaaa                                             <  1:732072/149‑1 (MQ=255)
                                                                                              gTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCCCTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGaaaa                                             <  2:393174/149‑1 (MQ=255)
                                                                                                 ttGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCCCTAAAGGTGGCATCATGACTGATCAAGCTGaaaaaaa                                          <  1:710231/149‑1 (MQ=255)
                                                                                                    cGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCACCACGACTGACCAAGCTGAAAAAAAGCa                                       >  2:644570/1‑149 (MQ=255)
                                                                                                       cACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCActc                                    >  2:279265/1‑149 (MQ=255)
                                                                                                        aCTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCGCTCTCTTACCCTACATCCTCAATAACAAAAATTGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCActct                                   >  2:433012/1‑149 (MQ=255)
                                                                                                        aCTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGATAGCCCCTCTCTTACCCTGCACACTTCCAAACACATATTCTCTTGCTTGAAATTGTGCATCATGACTGATCAAGCTGAAAAAAAGCGctcg                                   >  2:552552/1‑148 (MQ=255)
                                                                                                                  aaaaCTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCACTCTGCATTTTggg                         <  1:445202/149‑1 (MQ=255)
                                                                                                                      cTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCACTCTGCATTTTGGGgtgt                     >  1:719926/1‑149 (MQ=255)
                                                                                                                           agaagTTTAATTAATACTACAGTGGTGCTATAAGGATGTTAGCCACTCTCTTACCCTTCATCCTCAAAAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCACTCTGCATTTTGGGGTGTTATgg                <  2:13485/149‑1 (MQ=255)
                                                                                                                              agTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCACTCTGCATTTTGGGGTGTTATGGTta             <  1:455723/149‑1 (MQ=255)
                                                                                                                                         tactacAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAATGTGGCATCATGACTGATCAAGCTGAAAAAAAGCACTCTGCATTTTGGGGTGTTATGGTTATAGCAGGTACa  <  1:561581/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
TCAAACACATATGGACTTCATTATTGAAGCCTTTAAACATGTGAAAGAGAACGCGGCGAATATTAAAGGATTAACCTTTACGTACGAACCGAAAGTATTGCGTCACTTCACCGCAAAACTTAAAGAAGTTTAATTAATACTACAGAGTGGCTATAAGGATGTTAGCCACTCTCTTACCCTACATCCTCAATAACAAAAATAGCCTTCCTCTAAAGGTGGCATCATGACTGATCAAGCTGAAAAAAAGCACTCTGCATTTTGGGGTGTTATGGTTATAGCAGGTACA  >  NC_000913/3890012‑3890297

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑