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.1% 161.1 / ‑5.6 52Q105* (CAA→TAA) ycjGL‑Ala‑D/L‑Glu epimerase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (24/27);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CGCAAATTATGAGCGTCGTGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGCAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAA  >  NC_000913/1389108‑1389367
                                                                                                                                      |                                                                                                                             
cgcAAATTATGAGCGTCGTGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATc                                                                                                                           >  2:1785798/1‑139 (MQ=255)
  cAAATTATGAGCGTCGTGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCgc                                                                                                                         >  1:364377/1‑139 (MQ=255)
  cAAATTATGAGCGTCGTGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCgc                                                                                                                         >  1:235805/1‑139 (MQ=255)
   aaaTTATGAGCGTCGTGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCgct                                                                                                                        >  2:231484/1‑139 (MQ=255)
        aTGAGCGTCGTGCAGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCgctggctg                                                                                                                   >  1:331133/1‑139 (MQ=255)
             cgtcgtGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTg                                                                                                              >  2:1709791/1‑139 (MQ=255)
              gtcgtGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGa                                                                                                             <  2:235805/139‑1 (MQ=255)
               tcgtGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                            <  2:167769/139‑1 (MQ=255)
               tcgtGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGAt                                                                                                            >  1:1706115/1‑139 (MQ=255)
                cgtGCCGCAACTAAGGAAAGGGCTGACCCGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCGGCACGTAATCCGCTGGATTGGGCATCGGGGGATCTGGCCGCGCGAAGCCAGCAGCAATCGCTGGCTGATTTGATc                                                                                                           >  1:280539/1‑139 (MQ=255)
                       cAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATaa                                                                                                    >  2:69125/1‑139 (MQ=255)
                         aCTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACg                                                                                                  <  2:1706115/139‑1 (MQ=255)
                             gagaAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTc                                                                                              >  1:1360143/1‑139 (MQ=255)
                              agaAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTcc                                                                                             <  1:1785798/139‑1 (MQ=255)
                                aaaGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCATGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCACTAATCGCTGGCTGATTTGATCTGAATAACGCTTCCCg                                                                                           <  1:1257393/139‑1 (MQ=255)
                                  aGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCgag                                                                                         <  2:55736/139‑1 (MQ=255)
                                          cacGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTAt                                                                                 <  1:1418256/139‑1 (MQ=255)
                                                aggagTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGc                                                                           <  2:331133/139‑1 (MQ=255)
                                                     ttGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGa                                                                      <  2:1232196/139‑1 (MQ=255)
                                                     ttGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGCAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGa                                                                      <  2:280539/139‑1 (MQ=255)
                                                                     cTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTAc                                                      <  2:1360143/139‑1 (MQ=255)
                                                                           gcgcAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCgctggctg                                                                                                                   >  2:1515073/1‑72 (MQ=255)
                                                                           gcgcAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCgctggctg                                                                                                                   <  1:1515073/72‑1 (MQ=255)
                                                                             gcagcaCGTAATGCGCTGGATTGTGCATTGTGGGATCTGGTCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATc                                              <  2:692017/139‑1 (MQ=255)
                                                                                 caCGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATc                                                                                                                           <  1:40586/58‑1 (MQ=255)
                                                                                 caCGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATc                                                                                                                           >  2:40586/1‑58 (MQ=255)
                                                                                 caCGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACCGTTGTCATCGGTACGCCTGATCAGAt                                          >  1:1073728/1‑139 (MQ=255)
                                                                                  aCGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATg                                         >  2:694962/1‑139 (MQ=255)
                                                                                       aTGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCaa                                    <  2:944817/139‑1 (MQ=255)
                                                                                          cgcTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAg                                 <  2:364377/139‑1 (MQ=255)
                                                                                           gcTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGcaca                                                                        <  2:1647128/99‑1 (MQ=255)
                                                                                           gcTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGcaca                                                                        >  1:1647128/1‑99 (MQ=255)
                                                                                            cTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCa                                                                                                       >  1:1840649/1‑67 (MQ=255)
                                                                                            cTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCa                                                                                                       <  2:1840649/67‑1 (MQ=255)
                                                                                                  tgtgCATTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCaa                         <  1:1709791/139‑1 (MQ=255)
                                                                                                       aTTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACActc                    >  1:936838/1‑139 (MQ=255)
                                                                                                       aTTGTGGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACActc                    <  1:231484/139‑1 (MQ=255)
                                                                                                         tgtgGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTc                                                              >  1:654780/1‑95 (MQ=255)
                                                                                                         tgtgGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTc                                                              <  2:654780/95‑1 (MQ=255)
                                                                                                         tgtgGGATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTg                  >  2:215347/1‑139 (MQ=255)
                                                                                                              gATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTAt                                                                                 <  2:1221247/71‑1 (MQ=255)
                                                                                                              gATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTAt                                                                                 >  1:1221247/1‑71 (MQ=255)
                                                                                                              gATCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTggcagg             <  2:914886/139‑1 (MQ=255)
                                                                                                               aTCTGGCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACCGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTggcaggc            <  2:1073728/139‑1 (MQ=255)
                                                                                                                   ggCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCg                                                                                           >  2:1771599/1‑56 (MQ=255)
                                                                                                                   ggCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCg                                                                                           <  1:1771599/56‑1 (MQ=255)
                                                                                                                   ggCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGcaca                                                                        >  1:1751583/1‑75 (MQ=255)
                                                                                                                   ggCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGcaca                                                                        <  2:1751583/75‑1 (MQ=255)
                                                                                                                    gCCGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGgcg       <  1:694962/139‑1 (MQ=255)
                                                                                                                     ccGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAg                                            <  1:881358/101‑1 (MQ=255)
                                                                                                                     ccGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAg                                            >  2:881358/1‑101 (MQ=255)
                                                                                                                     ccGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGgcgc      >  2:1455599/1‑139 (MQ=255)
                                                                                                                     ccGCGCGAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGCCCATAATTCAACCAAAACTCTGGAAGGCAGgcgc      >  2:1798990/1‑139 (MQ=255)
                                                                                                                         gcgAAGACAGCAGTAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGaaa  >  1:1582089/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                             
CGCAAATTATGAGCGTCGTGCCGCAACTAGAGAAAGGGCTGACACGGGAGGAGTTGCAAAAAATTCTCCCTGCCGGCGCAGCACGTAATGCGCTGGATTGTGCATTGTGGGATCTGGCCGCGCGAAGACAGCAGCAATCGCTGGCTGATTTGATCGGCATAACGCTTCCCGAGACAGTTATCACTGCACAGACGGTTGTCATCGGTACGCCTGATCAGATGGCCAATAGTGCATCAACACTCTGGCAGGCAGGCGCGAAA  >  NC_000913/1389108‑1389367

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