Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,389,242 C→T Q105* (CAA→TAA)  ycjG → L‑Ala‑D/L‑Glu epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,389,2420CT98.3% 191.4 / ‑5.0 60Q105* (CAA→TAA) ycjGL‑Ala‑D/L‑Glu epimerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/31);  minor base A (1/0);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AAATTATGAGCGTCGTGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGCAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGA  >  NC_000913/1389111‑1389365
                                                                                                                                   |                                                                                                                           
aaaTTATGAGCGTCGTGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCgct                                                                                                                      >  2:1437725/1‑139 (MQ=255)
 aaTTATGAGCGTCGTGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCgctg                                                                                                                     >  2:890321/1‑139 (MQ=255)
            tcgtGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                          <  2:22039/139‑1 (MQ=255)
            tcgtGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                          >  2:1378457/1‑139 (MQ=255)
            tcgtGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                          >  2:21802/1‑139 (MQ=255)
              gtGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCg                                                                                                        >  1:609792/1‑139 (MQ=255)
                      aCTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACg                                                                                                <  1:1243647/139‑1 (MQ=255)
                               aGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCgag                                                                                       >  1:108374/1‑139 (MQ=255)
                               aGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCgag                                                                                       <  2:108374/139‑1 (MQ=255)
                               aGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCgag                                                                                       >  1:226057/1‑139 (MQ=255)
                                  gCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCgag                                                                                       <  1:1584436/136‑1 (MQ=255)
                                  gCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAAAGCTTCCCgag                                                                                       >  2:1584436/1‑136 (MQ=255)
                                    tGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCATGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCgctgg                                                                                                                    >  2:809489/1‑105 (MQ=255)
                                    tGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCgctgg                                                                                                                    <  1:809489/105‑1 (MQ=255)
                                     gACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGtt                                                                                 >  1:672695/1‑139 (MQ=255)
                                       cacGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTAt                                                                               <  1:1437725/139‑1 (MQ=255)
                                          gggAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCAc                                                                            <  2:1545587/139‑1 (MQ=255)
                                           ggaggaGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGTGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCgag                                                                                       <  1:933649/127‑1 (MQ=255)
                                           ggaggaGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCgag                                                                                       >  2:933649/1‑127 (MQ=255)
                                           ggaggaGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACt                                                                           >  1:351342/1‑139 (MQ=255)
                                                 gTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAg                                                                     >  2:963113/1‑139 (MQ=255)
                                                  ttGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGa                                                                    >  2:1191874/1‑139 (MQ=255)
                                                     cAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACgg                                                                 <  2:672695/139‑1 (MQ=255)
                                                     cAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACgg                                                                 <  1:1191874/139‑1 (MQ=255)
                                                     cAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACgg                                                                 <  1:890321/139‑1 (MQ=255)
                                                      aaaaaaTTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGt                                                                >  2:924583/1‑139 (MQ=255)
                                                          aaTTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGcaca                                                                      <  2:264076/129‑1 (MQ=255)
                                                          aaTTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGcaca                                                                      >  1:264076/1‑129 (MQ=255)
                                                            ttCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCAt                                                          <  2:351342/139‑1 (MQ=255)
                                                               tcCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCgg                                                       <  2:15551/139‑1 (MQ=255)
                                                                  cTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTAc                                                    <  1:924583/139‑1 (MQ=255)
                                                                    gCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                          <  1:172720/83‑1 (MQ=255)
                                                                    gCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                          >  2:172720/1‑83 (MQ=255)
                                                                            agcaCGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAg                                          <  1:1378457/139‑1 (MQ=255)
                                                                            agcaCGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAg                                          <  2:609792/139‑1 (MQ=255)
                                                                            agcaCGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAg                                          <  2:226057/139‑1 (MQ=255)
                                                                               aCGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAg                                          >  1:884983/1‑136 (MQ=255)
                                                                               aCGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAg                                          <  2:884983/136‑1 (MQ=255)
                                                                                   aaTGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                          <  1:649777/68‑1 (MQ=255)
                                                                                   aaTGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                          >  2:649777/1‑68 (MQ=255)
                                                                                       cgcTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAg                               <  1:870876/139‑1 (MQ=255)
                                                                                         cTGGATTGTGCATTGTGCGACCTGGGCGCGCGAAGACTGCAGTAATCGCTGCCTGATTTGATCGGCATAACGCCTCCCGAGACAGTTATCACTGCCCAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATACTg                             <  1:21802/139‑1 (MQ=255)
                                                                                            gATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                          <  1:1585722/59‑1 (MQ=255)
                                                                                            gATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                          >  2:1585722/1‑59 (MQ=255)
                                                                                                 tgCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTcc                                                                                           >  2:420466/1‑69 (MQ=255)
                                                                                                 tgCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTcc                                                                                           <  1:420466/69‑1 (MQ=255)
                                                                                                      tgtgGGATCTGGCCGCGCGAAGACAGCAGTAATCg                                                                                                                        >  1:1381845/1‑35 (MQ=255)
                                                                                                      tgtgGGATCTGGCCGCGCGAAGACAGCAGTAATCg                                                                                                                        <  2:1381845/35‑1 (MQ=255)
                                                                                                      tgtgGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGc                                                                                               <  1:1242476/60‑1 (MQ=255)
                                                                                                      tgtgGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGc                                                                                               >  2:1242476/1‑60 (MQ=255)
                                                                                                      tgtgGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTg                >  2:1558660/1‑139 (MQ=255)
                                                                                                      tgtgGGATCTGGCCGCGCGAAGACAGCAGAAATCGCTGGCTGATTGGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTg                >  1:931997/1‑139 (MQ=255)
                                                                                                           gATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                          >  1:1073873/1‑44 (MQ=255)
                                                                                                           gATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                          <  2:1073873/44‑1 (MQ=255)
                                                                                                               tGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGt                              >  2:634220/1‑116 (MQ=255)
                                                                                                               tGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGt                              <  1:634220/116‑1 (MQ=255)
                                                                                                                 gCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGgcg     >  1:541988/1‑139 (MQ=255)
                                                                                                                    gcgcgAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGa  <  1:1558660/139‑1 (MQ=255)
                                                                                                                       cgAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCAt                                                                                                    <  2:67777/38‑1 (MQ=255)
                                                                                                                       cgAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCAt                                                                                                    >  1:67777/1‑38 (MQ=255)
                                                                                                                                   |                                                                                                                           
AAATTATGAGCGTCGTGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGCAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGA  >  NC_000913/1389111‑1389365

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: