Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,991,528 T→A 7.7% intergenic (‑49/+26) 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,5280TA 7.7% 152.9 / 3.3 52intergenic (‑49/+26)vgrG‑I/hcpC‑Itype VI secretion system protein/hemolysin‑coregulated protein
Reads supporting (aligned to +/- strand):  ref base T (26/22);  new base A (2/2);  total (28/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.26e-01

GCGCTGATGCCCTCGCGGAGTTCGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTGGGTGTAGTTTTCATTGTCTTT  >  NC_002947/2991382‑2991673
                                                                                                                                                  |                                                                                                                                                 
gcgcTGATGCCCTCGCGGAGTTCGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATgg                                                                                                                                                 >  6:280702/1‑149 (MQ=32)
          catCGCGGAGTTCGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGATGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGACCGGGCGGGCATGGCCTGCCAggg                                                                                                                                       >  6:54274/3‑149 (MQ=12)
               cgGAGTTCGAACGACACCACTTCGAATGCAGCCTTGCCAACGCGCGCCTCAAAGCTGTATCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCAtt                                                                                                                                  <  7:344451/149‑1 (MQ=255)
                      cGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAg                                                                                                                           >  5:95729/1‑149 (MQ=33)
                                  cTTCGAATGCGGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGcggcgg                                                                                                               >  4:135450/1‑149 (MQ=21)
                                  cTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGcggcgg                                                                                                               <  4:151123/149‑1 (MQ=34)
                                         tGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGt                                                                                                             >  1:199397/1‑144 (MQ=33)
                                         tGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGt                                                                                                             <  2:199397/144‑1 (MQ=33)
                                                      aCGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCagag                                                                                           <  4:28280/149‑1 (MQ=32)
                                                              cTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGgt                                                                                   >  2:343013/1‑149 (MQ=32)
                                                                 aaaGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGtt                                                                                >  8:262887/1‑149 (MQ=34)
                                                                          gcgcAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATc                                                                       <  3:135450/149‑1 (MQ=34)
                                                                              aaaTCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGa                                                                   <  5:280702/149‑1 (MQ=34)
                                                                              aaaTCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGa                                                                   <  1:169094/149‑1 (MQ=34)
                                                                                         gCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATc                                                                       <  3:100882/134‑1 (MQ=34)
                                                                                         gCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATc                                                                       >  4:100882/1‑134 (MQ=34)
                                                                                         gCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGAt                                                        >  6:388339/1‑149 (MQ=34)
                                                                                               cATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCg                                                  <  5:388339/149‑1 (MQ=34)
                                                                                                aTAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGAt                                                        >  2:45089/1‑142 (MQ=34)
                                                                                                aTAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGAt                                                        <  1:45089/142‑1 (MQ=34)
                                                                                                aTAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGa                                                 <  2:382222/149‑1 (MQ=34)
                                                                                                 tAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCg                                                                                     <  7:25072/112‑1 (MQ=32)
                                                                                                 tAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCg                                                                                     >  8:25072/1‑112 (MQ=32)
                                                                                                 tAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGa                                                 >  2:24822/1‑148 (MQ=34)
                                                                                                  aaaaGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATa                                               >  4:136973/1‑149 (MQ=34)
                                                                                                      gATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGcc                                           <  6:95729/149‑1 (MQ=34)
                                                                                                      gATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGc                                            >  4:101639/1‑148 (MQ=34)
                                                                                                           ttGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGt                                                                     >  6:91135/1‑118 (MQ=35)
                                                                                                           ttGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGt                                                                     <  5:91135/118‑1 (MQ=35)
                                                                                                           ttGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGa                                                 >  2:215378/1‑138 (MQ=34)
                                                                                                           ttGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGa                                                 <  1:215378/138‑1 (MQ=34)
                                                                                                                cgccgTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCCTTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATGTTTCGATACGCCAGCGCTACt                                  >  3:43095/1‑148 (MQ=255)
                                                                                                                cgccgTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACtt                                 >  5:286622/1‑149 (MQ=34)
                                                                                                                    gTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGACGTCAGAGCCTTCGGTGtt                                                                                <  8:10687/98‑1 (MQ=37)
                                                                                                                    gTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGACGTCAGAGCCTTCGGTGtt                                                                                >  7:10687/1‑98 (MQ=37)
                                                                                                                        gAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGt                         >  5:410203/1‑149 (MQ=35)
                                                                                                                         aaTCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCCGGCGGCGGGTGCGCGCCAGGCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCa                           <  4:256114/146‑1 (MQ=21)
                                                                                                                         aaTCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAg                          <  6:286622/147‑1 (MQ=34)
                                                                                                                           tCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCTATACGCCAGCGATACTTCAACCAGTTgg                      <  1:24822/149‑1 (MQ=37)
                                                                                                                           tCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCGAGCGATACTTAAACCAGtttg                      >  6:289325/1‑147 (MQ=255)
                                                                                                                             aaCACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGANCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTGGgt                    >  1:130187/1‑149 (MQ=37)
                                                                                                                             aaCACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTGGgt                    >  1:125415/1‑149 (MQ=37)
                                                                                                                              acacAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGa                                      >  6:339408/1‑130 (MQ=35)
                                                                                                                              acacAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGa                                      <  5:339408/130‑1 (MQ=35)
                                                                                                                               cacaAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTGGgtg                   >  3:352288/1‑148 (MQ=37)
                                                                                                                                     gCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTGGGTGTAGtttt            <  4:352288/149‑1 (MQ=37)
                                                                                                                                     ccccGGGCAGGCCATGCCCGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGtt                        >  1:8645/2‑137 (MQ=11)
                                                                                                                                     ccccGGGCAGGCCATGCCCGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGtt                        <  2:8645/136‑1 (MQ=11)
                                                                                                                                      cccGGGCAGGCCATGCCCGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCAc                                                             >  4:134102/1‑99 (MQ=11)
                                                                                                                                      cccGGGCAGGCCATGCCCGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCAc                                                             <  3:134102/99‑1 (MQ=11)
                                                                                                                                       ccGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTGGGTGTAGTTTTCa          <  7:262887/149‑1 (MQ=37)
                                                                                                                                               gcATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTTAACCAGTTGGGTGTAGTTTTCATTGTCttt  >  3:38779/1‑149 (MQ=37)
                                                                                                                                                  |                                                                                                                                                 
GCGCTGATGCCCTCGCGGAGTTCGAACGACACCACTTCGAATGCAGCTTTGCCAACGAGTGCCTCAAAGCTGTAGCGCAAATCAGATTGGCTAGGCATAAAAGATCCTTGTCCGCCGTCTGAATCAACACAAAGCCCGGGCGGGCATGGCCTGCCCGGGGCATTTGCACAGGTCAGGCGGCGGGTGCGCGCCAGTCGTCAGAGCCTTCGGTGTTGGCTACATCGTGAGCCCACGTGATCTTTCGATACGCCAGCGATACTTCAACCAGTTGGGTGTAGTTTTCATTGTCTTT  >  NC_002947/2991382‑2991673

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