Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 92,906 G→A 6.9% R230C (CGT→TGT)  trpB ← tryptophan synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294792,9060GA 6.9% 163.9 / 3.0 58R230C (CGT→TGT) trpBtryptophan synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (25/29);  new base A (2/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

TGGCGGCATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGG  >  NC_002947/92762‑93051
                                                                                                                                                |                                                                                                                                                 
tGGCGGCATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCCGTCTGCCAGACAGcctg                                                                                                                                              >  2:855182/1‑149 (MQ=255)
 ggcggcATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCCGTCTGCCAGACAGcctgg                                                                                                                                             <  1:855182/150‑3 (MQ=255)
  gcggcATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGNAGCCGCNTCCGACGCAGGCNACCAGGCTGTCTGGCAGACGGccttcc                                                                                                                                            <  1:809211/150‑1 (MQ=255)
        tGCTTGTCGGTGTGTCAGCCGTGGCCACCGGCTTTAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAGTGCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTctct                                                                                                                                      <  1:222846/150‑1 (MQ=255)
             gTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCaa                                                                                                                                 <  1:777695/150‑1 (MQ=255)
                 gtgtgtACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGGGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTTTTGCAACTgg                                                                                                                             >  1:620760/1‑150 (MQ=255)
                 gcgtgtACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTgg                                                                                                                             <  2:481741/148‑1 (MQ=255)
                  tggctACGCCGTGGCCACCGTCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTggg                                                                                                                            <  2:620760/146‑1 (MQ=255)
                   gtgtACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGc                                                                                                                           <  2:752062/150‑1 (MQ=255)
                   gtgtACGCCGTGGCCACCGGCTTCAACGCCGATGAGCTGCACGCTCGGCTCTTCGAGGAACTAATGACAAAGGCCCATGGCATTGCAGCAACGGCGAAAACAAACAACTATCCGATCTGCCGTTCGACCTTCCTTCTCTTGCAACTGGGc                                                                                                                           >  2:1243612/1‑150 (MQ=255)
                      tACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTGCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCCGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACg                                                                                                                        <  1:1125787/150‑1 (MQ=255)
                       aCGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACgg                                                                                                                       >  2:438789/1‑150 (MQ=255)
                                ccACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTg                                                                                                                              <  1:425881/134‑1 (MQ=255)
                                ccACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTg                                                                                                                              >  2:425881/1‑134 (MQ=255)
                                 cACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGc                                                                                                             <  1:1243612/150‑1 (MQ=255)
                                            aCGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGAc                                                                                                  >  1:813565/1‑150 (MQ=255)
                                            aCGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGGATCTTGGCCGTAGATCGAc                                                                                                  >  2:582212/1‑150 (MQ=255)
                                              gCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTg                                                                                                <  2:86979/150‑1 (MQ=255)
                                                 gatgatCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTctct                                                                                                                                      <  2:1143607/109‑1 (MQ=255)
                                                 gatgatCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTctct                                                                                                                                      >  1:1143607/1‑109 (MQ=255)
                                                 gatgatCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGaa                                                                                             >  2:742665/1‑150 (MQ=255)
                                                   tgatCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGccttcct                                                                                                                                           <  2:923050/102‑1 (MQ=255)
                                                   tgatCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGccttcct                                                                                                                                           >  1:923050/1‑102 (MQ=255)
                                                       cTGCGCGCTCTGCTCCGCGAGGAGCTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTcgcg                                                                                       <  2:1036246/150‑1 (MQ=255)
                                                        tGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCCGGCTCTCTGGCAGACGGCCGTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCgcgc                                                                                      <  2:247231/150‑1 (MQ=255)
                                                        tGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCgcgc                                                                                      >  1:256852/1‑150 (MQ=255)
                                                        tGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCCGTCTGCCAGACAGcctg                                                                                                                                              <  1:57182/94‑2 (MQ=255)
                                                        tGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCCGTCTGCCAGACAGcctg                                                                                                                                              >  2:57182/1‑93 (MQ=255)
                                                          cACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCCTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCAc                                                                                    <  2:172995/150‑1 (MQ=255)
                                                                   ctctTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCggg                                                                           >  2:850982/1‑150 (MQ=255)
                                                                    tctTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCgggg                                                                          >  2:413020/1‑150 (MQ=255)
                                                                           ggAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTctct                                                                                                                 >  2:353009/1‑104 (MQ=255)
                                                                           ggAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTctct                                                                                                                 <  1:353009/104‑1 (MQ=255)
                                                                            gcaCTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGAACACCATGGCGGGGTACGGgtg                                                                  <  1:39133/148‑1 (MQ=255)
                                                                               cTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTgg                                                               <  1:413020/150‑1 (MQ=255)
                                                                               cTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTgg                                                               <  2:205836/150‑1 (MQ=255)
                                                                                tCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGgtgt                                                                 <  1:1112363/147‑1 (MQ=255)
                                                                                tCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGgtgt                                                                 >  2:1112363/1‑147 (MQ=255)
                                                                                         aGGCCCATGGCAGTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTTTCTTGCAACCGGGCAAGGTTCTCTTTGCCGATGATCCACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACg                                                     <  1:582212/150‑1 (MQ=255)
                                                                                             ccATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGc                                                 >  1:1157989/1‑150 (MQ=255)
                                                                                             ccATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGc                                                 <  2:256852/150‑1 (MQ=255)
                                                                                             ccATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGAATGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGc                                                 >  2:1129128/1‑150 (MQ=255)
                                                                                              cATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGcc                                                >  1:839401/1‑150 (MQ=255)
                                                                                              cATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGcc                                                >  2:46516/1‑150 (MQ=255)
                                                                                               aTGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCg                                               >  1:692033/1‑150 (MQ=255)
                                                                                                tGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGa                                              >  1:746587/1‑150 (MQ=255)
                                                                                                    aTTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAg                                          <  1:46516/150‑1 (MQ=255)
                                                                                                    aTTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAg                                          >  1:805602/1‑150 (MQ=255)
                                                                                                          gccgccgccGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGaa                                    <  1:193899/150‑1 (MQ=255)
                                                                                                               cgccgACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGgtgt                               >  1:272994/1‑150 (MQ=255)
                                                                                                                   gACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTg                                                  >  2:1156516/1‑127 (MQ=255)
                                                                                                                   gACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTg                                                  <  1:1156516/127‑1 (MQ=255)
                                                                                                                          gCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTg                    <  1:438789/150‑1 (MQ=255)
                                                                                                                          gCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTg                    <  2:839401/150‑1 (MQ=255)
                                                                                                                               catccTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAAc               >  1:887627/5‑150 (MQ=255)
                                                                                                                                          gCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACgg                                                                     >  1:948275/1‑85 (MQ=255)
                                                                                                                                          gCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACgg                                                                     <  2:948275/85‑1 (MQ=255)
                                                                                                                                            aGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAgg  <  1:31917/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TGGCGGCATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACCCAGTCACGCAGG  >  NC_002947/92762‑93051

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