Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,262,902 G→A 8.2% L21L (TTG→TTA)  pspG → phage shock protein G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,262,9020GA 8.2% 151.8 / 4.0 49L21L (TTG→TTA) pspGphage shock protein G
Reads supporting (aligned to +/- strand):  ref base G (28/17);  new base A (2/2);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCGGTGGTGTGGGTTATTAAGGCGATTAAAGCACCAAAAGTGCCGA  >  NC_000913/4262754‑4263050
                                                                                                                                                    |                                                                                                                                                    
cgGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTg                                                                                                                                                      <  1:181373/149‑1 (MQ=255)
                      aaTAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTgg                                                                                                                                >  4:51721/1‑149 (MQ=255)
                      aaTAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTgg                                                                                                                                >  3:205616/1‑149 (MQ=255)
                           attattTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGt                                                                                                                               <  4:66435/145‑1 (MQ=255)
                           attattTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGt                                                                                                                               >  3:66435/1‑145 (MQ=255)
                           attattTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTtgtg                                                                                                                           >  3:155167/1‑149 (MQ=255)
                            ttattTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTgg                                                                                                                          >  3:477776/1‑149 (MQ=255)
                                      cATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTACTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCAt                                                                                                                >  1:183267/1‑149 (MQ=255)
                                      cATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTACTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCAt                                                                                                                >  1:92444/1‑149 (MQ=255)
                                                       cAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTa                                                                                               <  1:303341/149‑1 (MQ=255)
                                                       cAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTa                                                                                               >  3:510879/1‑149 (MQ=255)
                                                        aGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTAt                                                                                              >  3:287600/1‑149 (MQ=255)
                                                         gCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATg                                                                                             >  2:539654/1‑149 (MQ=255)
                                                           aaGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCt                                                                                           >  4:238135/1‑149 (MQ=255)
                                                            agagATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTg                                                                                          >  2:556487/1‑149 (MQ=255)
                                                            agagATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTg                                                                                          >  4:15734/1‑149 (MQ=255)
                                                                aTTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCAt                                                                                      <  1:539654/149‑1 (MQ=255)
                                                                    cATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATtgat                                                                                  <  4:477776/149‑1 (MQ=255)
                                                                     atttttGGGAGAGCATCACGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTAGTTGTGGCGGCGGCCATTATGTTCCTCGGCGGTATGCTGGCATtgatg                                                                                 >  4:47384/4‑149 (MQ=255)
                                                                              agagCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTtg                                                                        >  2:417007/1‑149 (MQ=255)
                                                                                  catcatGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTtgctgc                                                                    >  4:363411/1‑149 (MQ=255)
                                                                                   atcatGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGcc                                                                   >  2:618785/1‑149 (MQ=255)
                                                                                     catGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGt                                                                 <  2:63997/149‑1 (MQ=255)
                                                                                         cTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGtt                                                             <  4:258153/149‑1 (MQ=255)
                                                                                           ggAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTtac                                                           >  1:295488/1‑149 (MQ=255)
                                                                                              actactTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTc                                                                                                         >  3:262797/1‑100 (MQ=255)
                                                                                              actactTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTc                                                                                                         <  4:262797/100‑1 (MQ=255)
                                                                                              actactTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTtactac                                                        >  1:302853/1‑149 (MQ=255)
                                                                                               ctactTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTactact                                                       >  3:135834/1‑149 (MQ=255)
                                                                                                tactTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTg                                                      >  2:89946/1‑149 (MQ=255)
                                                                                                 actTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCGTTGATGATTAAGTTGCTGCCGTGGTTactactcg                                                     >  3:506849/1‑147 (MQ=255)
                                                                                                  ctTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCAtt                                                                                                               >  1:204732/1‑90 (MQ=255)
                                                                                                  ctTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCAtt                                                                                                               <  2:204732/90‑1 (MQ=255)
                                                                                                     tttGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTACTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGAt                                                 <  2:92444/149‑1 (MQ=255)
                                                                                                     tttGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTACTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGAt                                                 <  2:183267/149‑1 (MQ=255)
                                                                                                          gATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCgg                                            <  3:122738/149‑1 (MQ=255)
                                                                                                          gATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCgg                                            >  2:585066/1‑149 (MQ=255)
                                                                                                           aTTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCggt                                           >  1:349375/1‑149 (MQ=255)
                                                                                                           aTTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCggt                                           <  2:124407/149‑1 (MQ=255)
                                                                                                                      tGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCGGTGGTGTGGGtta                                <  4:287600/149‑1 (MQ=255)
                                                                                                                        tCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCGGTGGTGTGGGttatt                              >  3:361246/1‑149 (MQ=255)
                                                                                                                           tGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCGGTGGTGTGGGTTATTAAg                           <  3:363411/149‑1 (MQ=255)
                                                                                                                            gCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCGGTGGTGTGGGTTATTAAgg                          <  1:585066/149‑1 (MQ=255)
                                                                                                                             cTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCGGTGGTGTGGGTTATTAAGGc                         <  1:556487/149‑1 (MQ=255)
                                                                                                                                aTGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCGGTGGTGTGGGTTATTAAGGCGAt                      <  1:247652/149‑1 (MQ=255)
                                                                                                                                      aCCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCGGTGGTGTGGGTTATTAAGGCGATTAAAGc                >  3:204570/1‑149 (MQ=255)
                                                                                                                                      aCCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCGGTGGTGTGGGTTATTAAGGCGATTAAAGc                >  4:617819/1‑149 (MQ=255)
                                                                                                                                         ggCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCGGTGGTGTGGGTTATTAAGGCGATTAAAGCAcc             <  3:466547/149‑1 (MQ=255)
                                                                                                                                             tttCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCgg                                            >  4:409309/1‑114 (MQ=255)
                                                                                                                                             tttCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCgg                                            <  3:409309/114‑1 (MQ=255)
                                                                                                                                                    gCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCGGTGGTGTGGGTTATTAAGGCGATTAAAGCACCAAAAGTGCCGa  >  1:233336/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGGGGCGTTTTGCTGTTAAATCAATAGATTATTTTTGGCATGATTCTTGTAATGCCAGCAAGAGATTTCATATTTGGGAGAGCATCATGCTGGAACTACTTTTTGTGATTGGCTTTTTTGTCATGCTGATGGTCACCGGCGTTTCGTTGCTGGGCATTATCGCCGCGCTGGTTGTGGCGACGGCCATTATGTTCCTCGGCGGTATGCTGGCATTGATGATTAAGTTGCTGCCGTGGTTACTACTGGCGATTGCGGTGGTGTGGGTTATTAAGGCGATTAAAGCACCAAAAGTGCCGA  >  NC_000913/4262754‑4263050

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