Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,999,680 (G)7→8 coding (212/756 nt) ygeR ← LysM domain‑containing M23 family putative peptidase, septation lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,999,6731.G95.9% 164.7 / ‑0.5 49coding (219/756 nt)ygeRLysM domain‑containing M23 family putative peptidase, septation lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (21/26);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTC  >  NC_000913/2999548‑2999788
                                                                                                                               |                                                                                                                    
ttCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATg                                                                                                         >  1:983946/1‑139 (MQ=255)
         gCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGt                                                                                                <  2:85140/139‑1 (MQ=255)
          cTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGGAAATGCCGTTCAGt                                                                                                <  1:655252/139‑1 (MQ=255)
          cTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTc                                                                                               >  1:819383/1‑139 (MQ=255)
          cTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTc                                                                                               >  1:811953/1‑139 (MQ=255)
             atgaCGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTcgcg                                                                                            >  2:604629/1‑139 (MQ=255)
                    tgtAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTc                                                                                     >  2:174451/1‑139 (MQ=255)
                       aaCCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCttt                                                                                  >  1:43021/1‑139 (MQ=255)
                         ccGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTa                                                                                >  1:1005022/1‑139 (MQ=255)
                           gATGCGGTTTTGGTCGTTGATTTGGCGGTTGCTTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACg                                                                              <  2:390776/139‑1 (MQ=255)
                            aTGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGc                                                                             >  2:377441/1‑139 (MQ=255)
                              gCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTg                                                                           >  1:317734/1‑139 (MQ=255)
                               cGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTgg                                                                          <  1:551326/139‑1 (MQ=255)
                                ggTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGt                                                                         >  1:918598/1‑139 (MQ=255)
                                    ttGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCg                                                                    <  2:183415/139‑1 (MQ=255)
                                       gTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgtg                                                                  >  1:427125/1‑139 (MQ=255)
                                         cGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAA‑‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgtggt                                                                <  2:817938/139‑1 (MQ=255)
                                            tGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTgcg                                                             <  1:558724/139‑1 (MQ=255)
                                            tGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTgcg                                                             <  1:604629/139‑1 (MQ=255)
                                              aTTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGcgcg                                                           <  1:174451/139‑1 (MQ=255)
                                              aTTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGcgcg                                                           <  2:1005022/139‑1 (MQ=255)
                                                       ttGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGat                                                  <  2:333335/139‑1 (MQ=255)
                                                          aTTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGAtata                                               <  1:377441/139‑1 (MQ=255)
                                                                    aaTACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATcccc                                     >  2:684925/1‑139 (MQ=255)
                                                                         tGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGcgcg                                                           >  1:92279/1‑112 (MQ=255)
                                                                         tGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGcgcg                                                           <  2:92279/112‑1 (MQ=255)
                                                                           cTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGaaa                                                        >  2:697139/1‑113 (MQ=255)
                                                                           cTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGaaa                                                        <  1:697139/113‑1 (MQ=255)
                                                                           cTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTAt                                         <  1:668766/128‑1 (MQ=255)
                                                                           cTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTc                              >  2:595197/1‑139 (MQ=255)
                                                                           cTACTTTTCGCCCCACCCAGTGTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTAGGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTAt                                         >  2:668766/1‑128 (MQ=255)
                                                                             aCTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCAc                            <  1:595197/139‑1 (MQ=255)
                                                                                   cGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTcgcg                                                                                            <  2:472946/69‑1 (MQ=255)
                                                                                   cGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTcgcg                                                                                            >  1:472946/1‑69 (MQ=255)
                                                                                   cGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTa                      >  1:1022071/1‑139 (MQ=255)
                                                                                       ccacccaGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACg                  <  2:317734/139‑1 (MQ=255)
                                                                                        cacccaGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACgg                 <  2:427125/139‑1 (MQ=255)
                                                                                           ccaGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGc              <  2:43021/139‑1 (MQ=255)
                                                                                           ccaGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGc              <  2:819383/139‑1 (MQ=255)
                                                                                           ccaGTTTTAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGc              <  2:983946/139‑1 (MQ=255)
                                                                                               ttttAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGa          <  1:684925/139‑1 (MQ=255)
                                                                                                  tAGTTTCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCTAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATa       <  2:212528/139‑1 (MQ=255)
                                                                                                       tCTGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTc  <  2:1022071/139‑1 (MQ=255)
                                                                                                         tGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgtggtg                                                               <  1:847623/76‑1 (MQ=255)
                                                                                                         tGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCgtggtg                                                               >  2:847623/1‑76 (MQ=255)
                                                                                                         tGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGc              >  2:198212/1‑125 (MQ=255)
                                                                                                         tGACCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGc              <  1:198212/125‑1 (MQ=255)
                                                                                                           aCCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCgcgc                                                                                           >  2:1076946/1‑46 (MQ=255)
                                                                                                           aCCAACTTCAATGGTGTAAGGGGGGGGAAATGCCGTTCAGTCgcgc                                                                                           <  1:1076946/46‑1 (MQ=255)
                                                                                                                               |                                                                                                                    
TTCGGTACCGCTGATGACGGTGTAACCGATGCGGTTTTGGTCGTTGATTTGGCGGTTGATTTACGTGTAATACTGCTACTTTTCGCCCCACCCAGTTTTAGTTTCTGACCAACTTCAATGGTGTAA‑‑GGGGGGGAAATGCCGTTCAGTCGCGCCAGTTCTTTTACGCTGGTTCCCGTGGTGCGCGAAATACGATATAGCGTATCCCCCCGTTTCACGGTGTAAACGGAGCCGGAATACGTTC  >  NC_000913/2999548‑2999788

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