Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,991,565:1 +A 11.8% intergenic (‑43/+32) vgrG‑I ← / ← hcpC‑I type VI secretion system protein/hemolysin‑coregulated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,991,5651.A11.8% 137.8 / 20.4 51intergenic (‑43/+32)vgrG‑I/hcpC‑Itype VI secretion system protein/hemolysin‑coregulated protein
Reads supporting (aligned to +/- strand):  ref base . (24/21);  new base A (3/3);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

AAGGGTGCGCTGATGCCCTCGCGGAGTTCGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTGGGTGTAGTTTTCATTGTCT  >  NC_002947/2991419‑2991714
                                                                                                                                                   |                                                                                                                                                       
aaGGGTGCGCTGATGCCCTCGCGGAGTTCGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCcgggcggg                                                                                                                                                      >  1:653603/1‑150 (MQ=32)
aaGGGTGCGCTGATGCCCTCGCGGAGTTCGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCcgggcggg                                                                                                                                                      >  1:879013/1‑150 (MQ=32)
        gcTGATGCCCTCGCGGAGTTCGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCt                                                                                                                                              >  2:408122/1‑150 (MQ=32)
        gcTGATGCCCTCGCGGAGTTCGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCt                                                                                                                                              >  1:860102/1‑150 (MQ=32)
         cTGATGCCCTCGCGGAGTTCGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACCCAAAGCCCGGGCGGGCATGGCCTG                                                                                                                                             <  2:973523‑M1/150‑1 (MQ=255)
              gCCCTCGCGGAGTTCGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGG                                                                                                                                       <  2:1114954‑M1/150‑1 (MQ=255)
               cccTCGCGGAGTTCGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGG                                                                                                                                      <  1:1113313‑M1/150‑1 (MQ=255)
                ccTCGCGGAGTTCGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGACTCCAGCTAGACAAAGCCCGGGCGGGCATGGCCTGCCCgggg                                                                                                                                     >  2:820012/1‑150 (MQ=11)
                          ttCGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGcaca                                                                                                                           >  1:908264/1‑150 (MQ=33)
                             gAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGt                                                                                                                        <  1:768780/150‑1 (MQ=35)
                                      cacTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCAGGCCATGCCCGCCCGGGctttgtgttgattcagacggcgg                                                                                                               >  1:694122‑M1/1‑127 (MQ=255)
                                      cacTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGgcggcg                                                                                                               >  2:424926/1‑150 (MQ=35)
                                       acTTCGAATGCAGCTTTGCCAACGCGTGCCTCTAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCAGGCCATGCCCGCCCGGGctttgtgttgattcagacggcgga                                                                                                              <  2:694122‑M1/150‑25 (MQ=255)
                                       acTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGcggcgg                                                                                                              <  2:987943/150‑1 (MQ=35)
                                       acTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGcggcgg                                                                                                              <  1:180885/150‑1 (MQ=35)
                                       acTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGcggcgg                                                                                                              <  1:270425/150‑1 (MQ=35)
                                         ttCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGt                                                                                                            <  1:737137/150‑1 (MQ=35)
                                                 cAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCa                                                                                                    <  1:334265/150‑1 (MQ=32)
                                                   gCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAgt                                                                                                  >  1:668652/1‑150 (MQ=32)
                                                    cTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAgtc                                                                                                 >  2:441338/1‑150 (MQ=32)
                                                    cTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAgtc                                                                                                 >  2:969820/1‑150 (MQ=32)
                                                     tttGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAgtcg                                                                                                >  1:1109657/1‑150 (MQ=32)
                                                          cAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCAGGCCATGCCCGCCCGGGcttt                                                                                                                                  <  1:9099‑M1/111‑5 (MQ=255)
                                                          cAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCAGGCCATGCCCGCCCGGGcttt                                                                                                                                  >  2:9099‑M1/1‑107 (MQ=255)
                                                                gTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAgtc                                                                                                 >  2:363710/1‑138 (MQ=33)
                                                                gTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAgtc                                                                                                 <  1:363710/138‑1 (MQ=33)
                                                                  gCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCgg                                                                                   <  2:668652/150‑1 (MQ=32)
                                                                   ccTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGgt                                                                                  <  1:441338/150‑1 (MQ=32)
                                                                     tCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCAGGCCATGCCCGCCCGGGctttgtgttgat                                                                                                                          <  1:793880‑M1/108‑13 (MQ=255)
                                                                     tCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCAGGCCATGCCCGCCCGGGctttgtgttgat                                                                                                                          >  2:793880‑M1/1‑96 (MQ=255)
                                                                     tCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGcggcgg                                                                                                              >  1:736949/1‑120 (MQ=33)
                                                                     tCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGcggcgg                                                                                                              <  2:736949/120‑1 (MQ=33)
                                                                     tCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGgtgt                                                                                >  1:394590/1‑150 (MQ=34)
                                                                     tCAAAGCTGTAGCGCAAATCAGATTGGCGAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGgtgt                                                                                >  1:1051664/1‑150 (MQ=21)
                                                                       aaaGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTg                                                                              >  2:1073880/1‑150 (MQ=35)
                                                                       aaaGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTg                                                                              <  2:908264/150‑1 (MQ=35)
                                                                                gcgcAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCg                                                                     >  1:875405/1‑150 (MQ=35)
                                                                                               gCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACAt                                                                       >  1:479406/1‑133 (MQ=35)
                                                                                               gCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATc                                                     >  2:817888/1‑150 (MQ=35)
                                                                                                                  tGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTCCTCGCGGTGCCCACGTGATCTTTCGATACGCCAGCGATa                                  <  1:1112577/150‑1 (MQ=11)
                                                                                                                      cgccgTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTc                              >  2:498428/1‑150 (MQ=35)
                                                                                                                       gccgTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCa                             >  2:988955/1‑150 (MQ=35)
                                                                                                                               aaTCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTg                     >  2:687764/1‑150 (MQ=38)
                                                                                                                                aTCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGcatttgtgttgattcagacggcggacaaggatcttttatgcctagccaatc                                                                                   <  2:520993‑M1/89‑52 (MQ=255)
                                                                                                                                aTCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGcatttgtgttgattcagacggcggacaaggatcttttatgcctagccaatc                                                                                   >  1:520993‑M1/1‑38 (MQ=255)
                                                                                                                                aTCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCATATCCTTCGGTGTTGGTTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTgg                    <  2:879013/150‑1 (MQ=255)
                                                                                                                                aTCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTgg                    <  1:498428/150‑1 (MQ=38)
                                                                                                                                     cacaAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGNTCAGNCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTAGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAAACAGATAGgtgga               >  1:51259/1‑148 (MQ=11)
                                                                                                                                         aaGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTGGGGGTAGttt           >  1:414921/1‑150 (MQ=37)
                                                                                                                                           gCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTGGGTGTAGTTTTc         <  1:408122/150‑1 (MQ=38)
                                                                                                                                             ccGGGCGGGCATGGCCTGCCCGGGctttgtgttgattcagacggcggacaaggatcttttatgcctagccaatctgatttgcgctacagctttgaggcactcgttggcaaagctgaattcgaagtggtgtcgttcgaactccgcg             >  2:740317‑M1/1‑24 (MQ=255)
                                                                                                                                             ccGGGCGGGCATGGCCTGCCCGGGctttgtgttgattcagacggcggacaaggatcttttatgcctagccaatctgatttgcgctacagctttgaggcactcgttggcaaagctgaattcgaagtggtgtcgttcgaactccgcg             <  1:740317‑M1/145‑122 (MQ=255)
                                                                                                                                             ccGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTGGGTGTAGTTTTCAt       <  1:653152/150‑1 (MQ=38)
                                                                                                                                               gggcgggcATGGCCTGCCCGGGctttgtgttgattcagacggcggacaaggatcttttatgcctagccaatctgatttgcgctacagctttg                                                                >  2:176614‑M1/1‑22 (MQ=255)
                                                                                                                                               gggcgggcATGGCCTGCCCGGGctttgtgttgattcagacggcggacaaggatcttttatgcctagccaatctgatttgcgctacagctttg                                                                <  1:176614‑M1/92‑71 (MQ=255)
                                                                                                                                                ggcgggcATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTGGGTGTAGTTTTCATTGt    <  2:1109657/150‑1 (MQ=39)
                                                                                                                                                  cgggcATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTGGGTGTAGTTTTCATTGTCt  <  1:1073880/150‑1 (MQ=38)
                                                                                                                                                   |                                                                                                                                                       
AAGGGTGCGCTGATGCCCTCGCGGAGTTCGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTGGGTGTAGTTTTCATTGTCT  >  NC_002947/2991419‑2991714

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: