Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 92,066 G→T 10.0% P105H (CCT→CAT)  trpA ← tryptophan synthase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294792,0660GT10.0% 113.3 / 2.9 40P105H (CCT→CAT) trpAtryptophan synthase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (20/16);  new base T (2/2);  total (22/18)
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.19e-01

GTCAGGCGGATGAAGTCGATGCCTGCGGCCTGGGCCGGGGTAGCCAGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCAT  >  NC_002947/91919‑92208
                                                                                                                                                   |                                                                                                                                              
gtCAGGCGGATGAAGTCGATGCCTGCGGCCTGGGCCGGGGTAGCCAGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAggg                                                                                                                                              >  1:534087/1‑150 (MQ=255)
   aGGCGGATGAAGTCGATGCCTGCGGCCTGGGCCGGGGTAGCCAGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTg                                                                                                                                           <  2:498082/150‑1 (MQ=255)
     gCGGATGAAGTCGATGCCTGCGGCCTGGGCCGGGGTAGCCAGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGCTCAGGCCCTAGACGCCCTCTTCCTTTTCTTGAGGCACTAACGTGTCCACGCCAAAGCGGTGGATAGGGTTgta                                                                                                                                         <  1:14265/150‑1 (MQ=255)
     gCGGATGAAGTCGATGCCTGCGGCCTGGGCCGGGGTAGCCAGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTgta                                                                                                                                         <  2:199664/150‑1 (MQ=255)
        gATGAAGTCGATGCCTGCGGCCTGGGCCGGGGTAGCCAGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTgtagta                                                                                                                                      <  2:534087/150‑1 (MQ=255)
               tCGATGCCTGCGGCCTGGGCCGGGGTAGCCAGTTCCGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTAACGAACGTGTCAACCCCAAAGCGGTGGCTAGGGTTGTAGTAGCCCATc                                                                                                                               >  1:176291/1‑150 (MQ=255)
                      cTGCGGCCTGGGCCGGGGTAGCCAGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGa                                                                                                                        >  2:548277/1‑150 (MQ=255)
                        gtggCCTGGGCCGCGGTCCCCACACCACCGTCGTGCGATGGCAACAGGTCGACGATGATCACCCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGACCACGCCAAAGCGGTGGATAGGGTTCTAGTAGCCCATCAGCACGATc                                                                                                                      <  2:39198/148‑1 (MQ=255)
                          ggCCTGGGCAGGCGTCGCCACTTCGGCGTCGACATTGCGCGCGAGGACCGCCAGGTTCAGGCCATTGACGCCCGCTTCCTTGGCTGGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCgg                                                                                                                    >  2:7912/1‑150 (MQ=255)
                              tGGGCCGGGGTAGCCAGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTg                                                                                                                >  2:494647/1‑150 (MQ=255)
                              tGGGCCGGGGTAGCCAGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTg                                                                                                                >  2:494580/1‑150 (MQ=255)
                                    ggggTCGCCAGTTCGGCTTCGGGCTCGGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTgg                                                                                                          <  2:233161/150‑1 (MQ=255)
                                      ggTAGCCAGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATATGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGtt                                                                                                        >  1:110492/1‑150 (MQ=255)
                                         aGCCAGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATc                                                                                                     <  2:241225/150‑1 (MQ=255)
                                          gCCAGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATATGGTTGTAGTAGCCCATCa                                                                                                                              <  2:243310/124‑1 (MQ=255)
                                          gCCAGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATATGGTTGTAGTAGCCCATCa                                                                                                                              >  1:243310/1‑124 (MQ=255)
                                           ccAGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCca                                                                                                   <  1:276158/150‑1 (MQ=255)
                                             aGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCcaca                                                                                                 >  1:636329/1‑150 (MQ=255)
                                                   gCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAAt                                                                                           >  1:309664/1‑150 (MQ=255)
                                                        gtGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTcgcg                                                                                      >  2:597711/1‑150 (MQ=255)
                                                            tctGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCg                                                                                                                     >  1:303812/1‑115 (MQ=255)
                                                            tctGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCg                                                                                                                     <  2:303812/115‑1 (MQ=255)
                                                              tGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCAt                                                                                >  1:720802/1‑150 (MQ=255)
                                                                gCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCt                                                                              <  1:51231/150‑1 (MQ=255)
                                                                  ggAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGTTCGCGCGGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGc                                                                            >  1:610435/1‑150 (MQ=255)
                                                                          gacgaTGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTGATCCACACGGAATTCGCGAACCATCTGCAGGGTCtt                                                                    >  2:550949/1‑150 (MQ=255)
                                                                             gatgatCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTgtagta                                                                                                                                      <  2:245778/81‑1 (MQ=255)
                                                                             gatgatCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTgtagta                                                                                                                                      >  1:245778/1‑81 (MQ=255)
                                                                                 atCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGCAGGGTCTTGACCAAg                                                             <  2:636329/150‑1 (MQ=255)
                                                                                        cATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGc                                                      >  2:619853/1‑150 (MQ=255)
                                                                                            gACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAggg                                                                                                                                              >  2:155242/1‑58 (MQ=255)
                                                                                            gACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAggg                                                                                                                                              <  1:155242/58‑1 (MQ=255)
                                                                                                 ccGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCaa                                             >  2:50176/1‑150 (MQ=255)
                                                                                                     ttCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGTGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGCAg                                                                          >  1:524122/1‑117 (MQ=255)
                                                                                                     ttCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGTGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGCAg                                                                          <  2:524122/117‑1 (MQ=255)
                                                                                                       ccTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGgcgc                                       >  1:720955/1‑150 (MQ=255)
                                                                                                              ttGAGTCACGAACGTGTCTACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGCAGGGTCTTGCCCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGt                                >  2:276776/1‑150 (MQ=255)
                                                                                                                                cACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCAc              <  2:176291/150‑1 (MQ=255)
                                                                                                                                cACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCAc              >  2:143960/1‑150 (MQ=255)
                                                                                                                                     cAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGc                                                  <  1:167365/109‑1 (MQ=255)
                                                                                                                                     cAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGc                                                  >  2:167365/1‑109 (MQ=255)
                                                                                                                                      aaaGCGGTGGATATGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCTACgccatc        <  2:110492/150‑1 (MQ=255)
                                                                                                                                            gTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACgccatcggccat  <  1:597711/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GTCAGGCGGATGAAGTCGATGCCTGCGGCCTGGGCCGGGGTAGCCAGTTCGGCGTCGTGCTCTGGCGGAAGGTCGACGATGATCAGGCCATCGACGCCCGCTTCCTTGGCTTGAGTCACGAACGTGTCCACGCCAAAGCGGTGGATAGGGTTGTAGTAGCCCATCAGCACGATCGGCGTGGTCTGGTTATCCACACGGAATTCGCGAACCATCTGCAGGGTCTTGACCAAGGTCTGGCCAGCTTCAAGGGCGCGCAGGGTGGCCAGCTGGATGGCCACGCCATCGGCCAT  >  NC_002947/91919‑92208

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: