Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 873,807 G→A G277S (GGT→AGT)  yliE → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913873,8070GA96.3% 180.8 / ‑1.7 54G277S (GGT→AGT) yliEputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (27/25);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

TATAATATTATTGATTCTGTCGATACTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCGGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGTTTATTTAATGACCACTCGTGT  >  NC_000913/873711‑873901
                                                                                                |                                                                                              
tataATATTATTGATTCTGTCGATACTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTa                                                                                             >  2:861902/1‑100 (MQ=255)
    atatTATTGATTCTGTCGATACTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAAttt                                                                                         >  1:1182860/1‑100 (MQ=255)
    atatTATTGATTCTGTCGATACTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAAttt                                                                                         >  1:364634/1‑100 (MQ=255)
     tattatTGATTCTGTCGATACTCAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTa                                                                                        >  2:759660/1‑100 (MQ=255)
      attattGATTCTGTCGATACTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTAt                                                                                       <  2:246694/100‑1 (MQ=255)
           tGATTCTGTCGATACTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAAc                                                                                  >  1:797444/1‑100 (MQ=255)
            gATTCTGTCGATACTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACg                                                                                 <  1:948063/100‑1 (MQ=255)
                cTGTCGATACTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATc                                                                             >  2:1034405/1‑100 (MQ=255)
                 tGTCGATACTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCa                                                                            <  1:937709/100‑1 (MQ=255)
                 tGTCGATACTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCa                                                                            >  1:973664/1‑100 (MQ=255)
                 tGTCGATACTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCa                                                                            <  2:256451/100‑1 (MQ=255)
                   tCGATACTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGa                                                                          <  1:982304/100‑1 (MQ=255)
                     gATACTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATa                                                                        >  2:303864/1‑100 (MQ=255)
                      aTACTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATaa                                                                       <  1:14038/100‑1 (MQ=255)
                         cTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACt                                                                    >  1:456100/1‑100 (MQ=255)
                            aaTATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACtttt                                                                 <  1:1073475/100‑1 (MQ=255)
                              tatatCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACtttttt                                                               <  1:1036419/100‑1 (MQ=255)
                                   cGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTg                                                          >  1:133455/1‑100 (MQ=255)
                                   cGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTg                                                          >  1:537608/1‑100 (MQ=255)
                                        aaaCATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGtt                                                     <  1:26944/100‑1 (MQ=255)
                                            aTCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCg                                                 <  1:768539/100‑1 (MQ=255)
                                              cAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGtt                                               <  1:1034415/100‑1 (MQ=255)
                                                  ggTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCatt                                           <  1:759665/100‑1 (MQ=255)
                                                    tGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCattat                                         >  2:12362/1‑100 (MQ=255)
                                                    tGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCattat                                         >  2:971149/1‑100 (MQ=255)
                                                    tGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCattat                                         <  1:351450/100‑1 (MQ=255)
                                                     gCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTGATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCattatt                                        >  2:725624/1‑100 (MQ=255)
                                                     gCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCattatt                                        <  1:1133936/100‑1 (MQ=255)
                                                        aaaCAACGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATc                                     >  1:928506/1‑100 (MQ=255)
                                                             acGCATTATTCGAAATCTATTTTTATTATCCTGGCGGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTAtt                                >  2:44312/1‑100 (MQ=255)
                                                              cGCATTATTCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTg                               >  1:596128/1‑100 (MQ=255)
                                                                cATTATTCGAAATCTATTTTTATTATCCTGGCGGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCt                             >  1:623050/1‑100 (MQ=255)
                                                                  ttattCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTTTTTTTTTTTTTCCTTCATTAGTATCGTCTTGCt                             >  2:451051/1‑98 (MQ=255)
                                                                  ttattCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCtgg                           <  1:1006523/100‑1 (MQ=255)
                                                                  ttattCGAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCtgg                           <  1:574092/100‑1 (MQ=255)
                                                                        gAAATCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGttt                     <  2:420832/100‑1 (MQ=255)
                                                                          aaTCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTGTTGCTGGTGGtttat                   >  2:475695/1‑100 (MQ=255)
                                                                          aaTCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGtttat                   <  1:12362/100‑1 (MQ=255)
                                                                           aTCTATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGtttatt                  <  2:133453/100‑1 (MQ=255)
                                                                              tATTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGTTTATTTaa               >  2:401689/1‑100 (MQ=255)
                                                                               aTTTTTATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGTTTATTTAAt              <  2:747046/100‑1 (MQ=255)
                                                                                tttttATTATTCTGGCAGTAGTTTATTGAGGGCGACAGCTAAACTTTTTTATCTGCCGTTTGGGTTGATTATTCTCGTATTACGGGTGCgttattggatg             >  2:623613/1‑95 (MQ=255)
                                                                                tttttATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTTTTTTTTTTTTTCGTTCAATACTACCTTATTGCTGCTGGtttatttatt              >  2:903403/1‑97 (MQ=255)
                                                                                tttttATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGTTTATTTAATg             <  1:980117/100‑1 (MQ=255)
                                                                                tttttATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGTTTATTTAATg             >  2:1090580/1‑100 (MQ=255)
                                                                                tttttATTATCCTGGCAGTAATTTATTGAACGCATCAGATAAAATTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGTTTATTTAATg             <  2:832417/100‑1 (MQ=255)
                                                                                       tatCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTCTTTATCTGCCGTTTGCGTTCATTATTATAGTATTGCTGGTGGTTTATTTAATGACCACTc      >  2:847677/1‑100 (MQ=255)
                                                                                        atCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGTTTATTTAATGACCACTCg     <  2:417784/100‑1 (MQ=255)
                                                                                        atCCTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGCTTTATTTAATGACCACTCg     >  2:734931/1‑100 (MQ=255)
                                                                                          ccTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGTTTATTTAATGACCACTCgtg   <  1:44653/100‑1 (MQ=255)
                                                                                          ccTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGTTTATTTAATGACCACTCgtg   <  1:1083249/100‑1 (MQ=255)
                                                                                          ccTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGTTTATTTAATGACCACTCgtg   >  2:894355/1‑100 (MQ=255)
                                                                                          ccTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGTTTATTTAATGACCACTCgtg   >  1:211642/1‑100 (MQ=255)
                                                                                           cTGGCAGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGTTTATTTAATGACCACTCgtgt  >  2:304417/1‑100 (MQ=255)
                                                                                                |                                                                                              
TATAATATTATTGATTCTGTCGATACTGAATATATCGAAAAAACATCAGCGGTGCCAAACAACGCATTATTCGAAATCTATTTTTATTATCCTGGCGGTAATTTATTGAACGCATCAGATAAACTTTTTTATCTGCCGTTTGCGTTCATTATTATCGTATTGCTGGTGGTTTATTTAATGACCACTCGTGT  >  NC_000913/873711‑873901

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