Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 2,921,460 G→A G3494R (GGG→AGG)  PP_2561 → hemolysin‑type calcium‑binding bacteriocin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,921,4600GA98.1% 178.6 / NA 53G3494R (GGG→AGG) PP_2561hemolysin‑type calcium‑binding bacteriocin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (33/19);  minor base T (0/1);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATC  >  NC_002947/2921316‑2921600
                                                                                                                                                |                                                                                                                                            
ggaacattGGGGGAGGCCGGTAACGACAACCTGACCGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCagg                                                                                                                                            <  6:280987/140‑1 (MQ=255)
   ttGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGATGGCGATGACAGGGTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGGAAGGACAGATTGGGCGGGGACGTGGGTACCGAGACCCTGATCagggggg                                                                                                                                        >  8:226988/1‑148 (MQ=255)
     gCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAAATGGACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGcc                                                                                                                                      >  1:274606/1‑148 (MQ=255)
       tGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGACTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGCCACCCTGATCAGGGGGGCggg                                                                                                                                    >  1:135244/1‑148 (MQ=255)
                  tgTAACGACAACCTGAACGGGGGCGCTGGAGATGAAAGCATGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGAGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTcc                                                                                                                         >  2:101570/2‑148 (MQ=255)
                  ggTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGAGGGCGGCAACGGCAATGACAGCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTc                                                                                                                          <  5:170701/147‑1 (MQ=255)
                   cgaaccccAACCTGAACGCGGGCGCTGGCGATGACATCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCt                                                                                                                        <  4:110690/141‑1 (MQ=255)
                       cGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAAc                                                                                                                    >  7:27782/1‑148 (MQ=255)
                              cTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCCGCGATGGTAACGACACATTGTTAGGCGACGTGGGTACTGACACCCTGATCAGGGGGTCCGGCACCGACTTCCTCAACGGTTCgg                                                                                                             >  8:155385/1‑148 (MQ=255)
                                       cacaatGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAAc                                                                                                     <  5:194505/142‑1 (MQ=255)
                                               cGATGACACCATGGATGGAGTCTACGGCTATGACACCCTGCAGGGCGGCGTTGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACggtg                                                                                            <  2:137451/148‑1 (MQ=255)
                                                gATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACggtgg                                                                                           <  3:175155/148‑1 (MQ=255)
                                                 aTGACAGCCTGGATGGCGGCAACGGCAATGTCACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACggt                                                                                             <  4:27535/145‑1 (MQ=255)
                                                 aTGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACggt                                                                                             >  3:27535/1‑145 (MQ=255)
                                                      aGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTggtgggtg                                                                                     >  8:273218/1‑148 (MQ=255)
                                                       gCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGgtgggtgg                                                                                    <  8:27782/148‑1 (MQ=255)
                                                         cTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTgggg                                                                                  >  3:188674/1‑148 (MQ=255)
                                                            ggtGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCCACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCg                                                                               >  3:215434/3‑148 (MQ=255)
                                                            gATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGcc                                                                                <  7:273218/147‑1 (MQ=255)
                                                                gctgcaacggcaaGGACACCCTGCAGTGCGGCGATGGTAACGACACATTGTTCTGCTACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTa                                                                            <  2:274606/144‑1 (MQ=255)
                                                                gctcccccggcaaTTTCACCCTGCACGGGGGCTCTGGTAACGACCCATTGTTCGGCGACGTGGGTACCGACCCCCTCATCTGGTTGGGCGGCAACGACATCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTaa                                                                           <  2:27631/141‑1 (MQ=255)
                                                                gcggcAACGGCAATGACCCCCTGCAGGGCGTCGATGGTAACGACACATTGTTCGGCGACGTGTGTACTGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGTTGGTGGGTGGGGCCGGTaa                                                                           <  7:155385/148‑1 (MQ=255)
                                                                        ggcaaTGACACCCTGCAGGGCGGCTATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACg                                                                                                    <  4:264182/115‑1 (MQ=255)
                                                                        ggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACg                                                                                                    >  3:264182/1‑115 (MQ=255)
                                                                         gcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCAt                                                                  >  7:266203/1‑148 (MQ=255)
                                                                            aTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGa                                                                >  6:310363/1‑147 (MQ=255)
                                                                               acacCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCg                                                             >  5:127477/1‑147 (MQ=255)
                                                                                cacCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGt                                                                                                                 <  3:199032/94‑1 (MQ=255)
                                                                                 acCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGcc                                                           >  6:214762/1‑147 (MQ=255)
                                                                                 acCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGCGTACCGACACCCTGATCAGGGCGGCCGGCAACGACTTCATCAACGGTg                                                                                                                >  4:199032/1‑94 (MQ=255)
                                                                                    cTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACCCCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCcccc                                                        >  4:39236/1‑147 (MQ=255)
                                                                                    cTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCAc                                                         >  5:112030/1‑146 (MQ=255)
                                                                                       cAGGGCGGCGAGGGTACCCACACATTGTTCTGTTACTTTTGTACCGCCACCCTCCTCATGGTTGCCGGCAAAGACTTCCTCAACGGTTCCGCAGGTTACGATTCGGTGGTGGGTGTGGCCGGTAACGACACCATGATCGCCACCGATg                                                    <  2:135244/148‑1 (MQ=255)
                                                                                           gcggcgATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGaaa                                                >  2:47810/1‑148 (MQ=255)
                                                                                             ggcgATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCgg                                                                              <  2:122271/116‑1 (MQ=255)
                                                                                             ggcgATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCgg                                                                              >  1:122271/1‑116 (MQ=255)
                                                                                              gcgATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACg                                              <  5:310363/147‑1 (MQ=255)
                                                                                                       aCGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTcc                                     >  6:200104/1‑147 (MQ=255)
                                                                                                       aCGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTc                                      >  2:135033/1‑146 (MQ=255)
                                                                                                           cacaTTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGtt                                 >  5:144769/1‑147 (MQ=255)
                                                                                                               ttGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAACCGATGTGTTCCAGTTCGCTg                            >  4:77141/1‑148 (MQ=255)
                                                                                                                 gTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGc                           <  8:266203/147‑1 (MQ=255)
                                                                                                                   tCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAg                         <  5:214762/147‑1 (MQ=255)
                                                                                                                       cGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCttt                    >  4:75328/1‑148 (MQ=255)
                                                                                                                              ggTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACCCCATGATCGCCCCCGATGGAAAGGAGGTGTTCCAGTTCGCTGCAGGCTTGGGCATCg             >  4:253425/1‑148 (MQ=255)
                                                                                                                              ggTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAAc              >  1:258271/1‑147 (MQ=255)
                                                                                                                                tACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGTTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGa            <  6:132059/147‑1 (MQ=255)
                                                                                                                                    gACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGa       >  1:329803/1‑148 (MQ=255)
                                                                                                                                      cacCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGat      >  2:62350/1‑147 (MQ=255)
                                                                                                                                      cacCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCGGCAGGCTTTGGCAACGACCTGat      >  1:328182/1‑147 (MQ=255)
                                                                                                                                      cacCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGAAACCATGATCGCCCGCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGat      >  2:202496/1‑147 (MQ=255)
                                                                                                                                        cccTGATGAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGTAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcat   <  4:188674/148‑1 (MQ=255)
                                                                                                                                          cTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCCCCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGataatc  >  7:146900/1‑147 (MQ=255)
                                                                                                                                          cTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcatc  >  2:141666/1‑147 (MQ=255)
                                                                                                                                          cTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcat   >  6:169063/1‑146 (MQ=255)
                                                                                                                                                |                                                                                                                                            
ACGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATC  >  NC_002947/2921316‑2921600

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