Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 92,900 G→C 7.3% P232A (CCA→GCA)  trpB ← tryptophan synthase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294792,9000GC 7.3% 159.0 / 2.6 51P232A (CCA→GCA) trpBtryptophan synthase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (19/28);  new base C (2/2);  total (21/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 = 1.00e+00

TGGCGGCATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACC  >  NC_002947/92762‑93039
                                                                                                                                          |                                                                                                                                           
tGGCGGCATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCCGTCTGCCAGACAGcctg                                                                                                                                  <  1:15436/150‑2 (MQ=255)
  gcggcATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGccttcc                                                                                                                                <  2:349089/150‑1 (MQ=255)
     gcATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCcttccttc                                                                                                                             >  2:700972/1‑150 (MQ=255)
       aTGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTtctc                                                                                                                           >  1:404900/1‑150 (MQ=255)
          cTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGAGCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTCCTCTTg                                                                                                                        <  1:295128/150‑1 (MQ=255)
              tCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAAc                                                                                                                    <  2:256000/150‑1 (MQ=255)
                 gtgtgtACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACtag                                                                                                                 <  2:713651/150‑3 (MQ=255)
                 gtgtgtACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCATTCTCTTGCAACTgg                                                                                                                 <  2:1050030/150‑1 (MQ=255)
                             tGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTctct                                                                                                     >  1:1067544/1‑150 (MQ=255)
                                  aCCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGcc                                                                                                >  2:415430/1‑150 (MQ=255)
                                            aCGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGTCACGGGTCTCTTTGCCGATGATCGAc                                                                                      <  1:1146969/150‑1 (MQ=255)
                                            aCGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCCGTCTGCCAGACAGcctg                                                                                                                                  <  2:953473/106‑2 (MQ=255)
                                            aCGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCCGTCTGCCAGACAGcctg                                                                                                                                  >  1:953473/1‑105 (MQ=255)
                                             cGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACt                                                                                     <  2:404900/150‑1 (MQ=255)
                                                cGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGa                                                                                  >  1:683781/1‑150 (MQ=255)
                                                 gatgatCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGaa                                                                                 <  1:1142540/150‑1 (MQ=255)
                                                       cTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCCCCCCCCCCCCCCCCCCGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTcgcg                                                                           >  1:179295/1‑150 (MQ=255)
                                                        tGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCga                                                                                              >  1:596977/1‑130 (MQ=255)
                                                        tGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCga                                                                                              <  2:596977/130‑1 (MQ=255)
                                                        tGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACACAGGCAACCAGGCTGTCTGGCAGACGGcc                                                                                                                                    >  1:896200/1‑92 (MQ=255)
                                                        tGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACACAGGCAACCAGGCTGTCTGGCAGACGGcc                                                                                                                                    <  2:896200/92‑1 (MQ=255)
                                                         gCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACgg                                                                                                           <  2:755385/116‑1 (MQ=255)
                                                         gCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACgg                                                                                                           >  1:755385/1‑116 (MQ=255)
                                                          cACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCAc                                                                        <  2:690259/150‑1 (MQ=255)
                                                           aCGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCAcc                                                                       >  2:646871/1‑150 (MQ=255)
                                                              cTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATg                                                                    <  2:1086244/150‑1 (MQ=255)
                                                                cGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGc                                                                  <  2:899953/150‑1 (MQ=255)
                                                                         gAGGAACGCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACgg                                                         <  2:3496/150‑1 (MQ=255)
                                                                           ggAACTCATGAACCAGGCCCAAGGCATCGACGCTGTGGCCGACGCGGTCTAGCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGgt                                                       <  2:615672/150‑1 (MQ=255)
                                                                              aCTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGtgtg                                                    <  1:255248/150‑1 (MQ=255)
                                                                               cTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTgg                                                   <  1:363486/150‑1 (MQ=255)
                                                                                  aTGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACggg                                                        <  2:820052/142‑1 (MQ=255)
                                                                                  aTGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACggg                                                        >  1:820052/1‑142 (MQ=255)
                                                                                  aTGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGgcc                                                >  1:803049/1‑150 (MQ=255)
                                                                                           gCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTg                                                                                    <  2:1021943/105‑1 (MQ=255)
                                                                                           gCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTg                                                                                    >  1:1021943/1‑105 (MQ=255)
                                                                                                       ggAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTa                           <  2:803049/150‑1 (MQ=255)
                                                                                                        gAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGtgtg                                                    <  2:441223/124‑1 (MQ=255)
                                                                                                        gAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGtgtg                                                    >  1:441223/1‑124 (MQ=255)
                                                                                                               cgccgACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTgg                                                                                                                 >  2:1205084/1‑56 (MQ=255)
                                                                                                               cgccgACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTgg                                                                                                                 <  1:1205084/56‑1 (MQ=255)
                                                                                                                   gACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTc               <  1:700972/150‑1 (MQ=255)
                                                                                                                      gCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCa                               >  1:880401/1‑131 (MQ=255)
                                                                                                                      gCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCa                               <  2:880401/131‑1 (MQ=255)
                                                                                                                        aGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGt          >  2:882284/1‑150 (MQ=255)
                                                                                                                         ggCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGtt         >  2:1090742/1‑150 (MQ=255)
                                                                                                                            aaCCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGc                                                                                                 >  2:106442/1‑59 (MQ=255)
                                                                                                                            aaCCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGc                                                                                                 <  1:106442/59‑1 (MQ=255)
                                                                                                                              ccAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTaa    <  1:415430/150‑1 (MQ=255)
                                                                                                                               cAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAAc   >  2:160469/1‑150 (MQ=255)
                                                                                                                                aGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAAcc  <  2:179295/150‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                           
TGGCGGCATGCTTGTCGGTGTGTACGCCGTGGCCACCGGCTTCAACGCCGATGATCTGCACGCTCGGCTCTTCGAGGAACTCATGAAACAGGCCCATGGCATTGGAGCCGCCGCCGACGCAGGCAACCAGGCTGTCTGGCAGACGGCCTTCCTTCTCTTGCAACTGGGCACGGGTCTCTTTGCCGATGATCGACTGGAAGTCGCGCACCATGGCGGGGTACGGGTGTGGGCCAGCCACGGTGCCGATCAGGTAGAAGGTGTCTTCGACGTTGGTAACC  >  NC_002947/92762‑93039

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