Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 511,638 C→T 7.0% S204S (AGC→AGT)  trpC → indole‑3‑glycerol phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947511,6380CT 7.0% 168.4 / 2.6 57S204S (AGC→AGT) trpCindole‑3‑glycerol phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base T (2/2);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

GCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATCAACGAGGTTTACTCGTTCCTGGTGGGTGAGGCGTTCATGCGCGCC  >  NC_002947/511492‑511785
                                                                                                                                                  |                                                                                                                                                   
gCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCt                                                                                                                                                   >  5:233695/1‑149 (MQ=255)
gCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAAACCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCt                                                                                                                                                   >  8:247663/1‑149 (MQ=255)
     aCTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGaaa                                                                                                                                              <  1:7775/149‑1 (MQ=255)
                cTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGAc                                                                                                                                   <  4:159972/149‑1 (MQ=255)
                   ccAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCt                                                                                                                                       >  5:277645/1‑142 (MQ=255)
                   ccAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCt                                                                                                                                       <  6:277645/142‑1 (MQ=255)
                     aaGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGTCTGGAAACCACCCTTGACctgct                                                                                                                              >  6:13177/1‑149 (MQ=255)
                     aaGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGTCTGGAAACCACCCTTGACctgct                                                                                                                              <  5:13177/149‑1 (MQ=255)
                     aaGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGTCTGGAAACCACCCTTGACctgct                                                                                                                              <  5:18261/149‑1 (MQ=255)
                     aaGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGTCTGGAAACCACCCTTGACctg                                                                                                                                >  6:18261/1‑147 (MQ=255)
                                ccTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTc                                                                                                                   >  5:162590/1‑149 (MQ=255)
                                 cTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTcc                                                                                                                  >  2:21526/1‑149 (MQ=255)
                                          cTGGTGGAAGTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGAcc                                                                                                          <  4:277913/148‑1 (MQ=255)
                                                   gTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCCCACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGAt                                                                                                >  6:319289/1‑149 (MQ=255)
                                                   gTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGAt                                                                                                <  1:21526/149‑1 (MQ=255)
                                                    tGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGCATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGAtt                                                                                               <  2:104533/149‑1 (MQ=255)
                                                              cGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCg                                                                                     <  6:162590/149‑1 (MQ=255)
                                                              cGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCg                                                                                     <  3:53853/149‑1 (MQ=255)
                                                               gATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCgg                                                                                    >  1:346118/1‑149 (MQ=255)
                                                                     cTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCg                                                                                     <  7:124106/142‑1 (MQ=255)
                                                                     cTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATtct                                                                              <  7:116962/149‑1 (MQ=255)
                                                                     cTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCAGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCg                                                                                     >  8:124106/1‑142 (MQ=255)
                                                                      tGGAGCGCGCGCTGAAAACCCTGGATGCGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTgg                                                                                                    >  5:219568/1‑126 (MQ=255)
                                                                      tGGAGCGCGCGCTGAAAACCCTGGATGCGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTgg                                                                                                    <  6:219568/126‑1 (MQ=255)
                                                                        gAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATtctc                                                                             <  1:117775/147‑1 (MQ=255)
                                                                        gAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCaa                                                                           <  2:346118/149‑1 (MQ=255)
                                                                            gcgcgcTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCgg                                                                       <  7:247663/149‑1 (MQ=255)
                                                                                 cTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGa                                                                  <  4:170744/149‑1 (MQ=255)
                                                                                 cTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGa                                                                  >  3:170744/1‑149 (MQ=255)
                                                                                          cTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCt                                                         >  5:277243/1‑149 (MQ=255)
                                                                                               tACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGCTCAGCCTCGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTCGAGCTGATgg                                                    >  1:172453/1‑149 (MQ=255)
                                                                                                    cgcTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATc                                               <  6:233695/149‑1 (MQ=255)
                                                                                                       tggtggGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGa                                                                                                 >  8:47640/1‑96 (MQ=255)
                                                                                                         gtggGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCACATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCaa                                                 <  1:325630/142‑1 (MQ=255)
                                                                                                         gtggGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCACATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCaa                                                 >  2:325630/1‑142 (MQ=255)
                                                                                                              ggTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGa                                                                  <  1:92725/120‑1 (MQ=255)
                                                                                                              ggTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGa                                                                  >  2:92725/1‑120 (MQ=255)
                                                                                                                  aacaacCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTATCAACCGGGCCGATGTGTAGCTGATGGCAATCAAAGAGGTTTaga                                  >  8:227224/1‑146 (MQ=255)
                                                                                                                  aacaacCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATGCCGCGTGACCGTCTGGCGATTACCGAAAGCGTTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATCAACGAGGtt                                      <  7:227224/144‑1 (MQ=255)
                                                                                                                   acaacCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGa                                                            <  2:86288/121‑1 (MQ=255)
                                                                                                                     aacCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTgg                                                             >  5:51667/1‑118 (MQ=255)
                                                                                                                     aacCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTgg                                                             <  6:51667/118‑1 (MQ=255)
                                                                                                                     aacCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATCAACGAGGTTTACTCGtt                              >  8:27048/1‑149 (MQ=255)
                                                                                                                             ccTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATCAACGAGGTTTACTCGTTCCtggtgg                      >  1:98701/1‑149 (MQ=255)
                                                                                                                                gCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGACTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATCAACGAGGTTTACTCGTTCCTggtgggt                    >  5:310853/1‑148 (MQ=255)
                                                                                                                                    acCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATCAACGAGGTTTACTCGTTCCTGGTGGGTGa                  <  4:17056/146‑1 (MQ=255)
                                                                                                                                    acCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATCAACGAGGTTTACTCGTTCCTGGTGGGTGa                  >  3:17056/1‑146 (MQ=255)
                                                                                                                                       ttCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCt                                                         >  3:20310/1‑104 (MQ=255)
                                                                                                                                       ttCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCt                                                         <  4:20310/104‑1 (MQ=255)
                                                                                                                                       ttCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATCAACGAGGTTTACTCGTTCCtggtg                       >  4:230108/1‑138 (MQ=255)
                                                                                                                                       ttCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATCAACGAGGTTTACTCGTTCCtggtg                       <  3:230108/138‑1 (MQ=255)
                                                                                                                                         cGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGAGCGTCTGGCGTTTACCGAAAGCGGTATTCTCAACCGGGCCGTTGTGGAGCTGATGGCAATCAACGAGGTTTACTCGTTCCTGGTGGGTGAGGCGtt            <  6:110616/147‑1 (MQ=255)
                                                                                                                                          gAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATCAACGAGGTTTACTCGTTCCTGGTGGGTGAGGCGTTCa          <  6:310853/148‑1 (MQ=255)
                                                                                                                                          gAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATCAACGAGGTTTACTCGTTCCTGGTGGGTGAGGCGTTCAt         >  4:118465/1‑149 (MQ=255)
                                                                                                                                          gAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCAGGCCGATGTGGAGCTGATGGCAATCAACGAGGTTTACTCGTTCCTGGTGGGTGAGGCGTTCAt         <  5:319289/149‑1 (MQ=255)
                                                                                                                                             gTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATCAACGa                                          >  2:215350/1‑113 (MQ=255)
                                                                                                                                             gTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATCAACGa                                          <  1:215350/113‑1 (MQ=255)
                                                                                                                                                  ccTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATCAACGAGGTTTACTCGTTCCTGGTGGGTGAGGCGTTCATGCGCGcc  >  4:344620/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAACTGGAGCGCGCGCTGAAAACCCTGGATACGCCGCTGGTGGGGGTGAACAACCGCAACCTGCACACCTTCGAGGTCAGCCTGGAAACCACCCTTGACCTGCTGCCGCGCATTCCGCGTGACCGTCTGGCGATTACCGAAAGCGGTATTCTCAACCGGGCCGATGTGGAGCTGATGGCAATCAACGAGGTTTACTCGTTCCTGGTGGGTGAGGCGTTCATGCGCGCC  >  NC_002947/511492‑511785

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