Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,564,769 A→G 9.9% intergenic (‑865/‑33) vgrG‑III ← / → glgA type VI secretion system protein/glycogen synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,564,7690AG 9.9% 164.9 / 7.7 61intergenic (‑865/‑33)vgrG‑III/glgAtype VI secretion system protein/glycogen synthase
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base G (3/3);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

CAGCCACGGCCACCAACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCGGCGCCCAGCGCCAGCACAA  >  NC_002947/4564621‑4564917
                                                                                                                                                    |                                                                                                                                                    
cAGCCACGGCCACCAACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAg                                                                                                                                                     <  1:1330846/150‑1 (MQ=255)
           accaACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCTGGCACCCgtggc                                                                                                                                           <  2:1230697/149‑5 (MQ=255)
           accaACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCTGGCACCCgtggc                                                                                                                                           >  1:1230697/1‑145 (MQ=255)
           accaACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCaa                                                                                                                                          >  2:1220607/1‑150 (MQ=255)
                   aacAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAggg                                                                                                                                  >  2:821555/1‑150 (MQ=255)
                             tCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGaa                                                                                                                        <  1:1158768/150‑1 (MQ=255)
                             tCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGaa                                                                                                                        <  1:463213/150‑1 (MQ=255)
                             tCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGaa                                                                                                                        <  1:1220607/150‑1 (MQ=255)
                               gCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGt                                                                                                                      <  1:821555/150‑1 (MQ=255)
                                        cGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGAACCCGACCAAGCATAGGGACGGAGAGAAGTATGATAAGt                                                                                                             >  1:75687/1‑150 (MQ=255)
                                           aaCGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGcc                                                                                                          <  1:710498/150‑1 (MQ=255)
                                              gCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATc                                                                                                                >  1:1047017/1‑141 (MQ=255)
                                              gCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATc                                                                                                                <  2:1047017/141‑1 (MQ=255)
                                                cGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGt                                                                                                     <  2:1217495/150‑1 (MQ=255)
                                                cGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGt                                                                                                     >  1:1189181/1‑150 (MQ=255)
                                                cGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGt                                                                                                     <  1:927518/150‑1 (MQ=255)
                                                cGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGt                                                                                                     <  1:40887/150‑1 (MQ=255)
                                                     gcgcCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGc                                                                                                <  2:1189181/150‑1 (MQ=255)
                                                        tcgtgcgcgCACGCACGGCACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACGTCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGc                                                                                                <  1:25561/146‑1 (MQ=255)
                                                        tcgtgcgcgCACGCACGGCACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACGTCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGc                                                                                                >  2:25561/2‑147 (MQ=255)
                                                        ccGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTc                                                                                             <  1:5632/150‑1 (MQ=255)
                                                                gcgcgcACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGAc                                                                                     <  2:714553/150‑1 (MQ=255)
                                                                    gcACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCCGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATc                                                                                                                >  1:1381028/1‑119 (MQ=255)
                                                                    gcACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATc                                                                                                                <  2:1381028/119‑1 (MQ=255)
                                                                      aCGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCa                                                                                                                                                      >  1:327360/1‑79 (MQ=255)
                                                                      aCGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCa                                                                                                                                                      <  2:327360/79‑1 (MQ=255)
                                                                           acaGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAg                                                                                                                                    <  1:1328333/92‑1 (MQ=255)
                                                                           acaGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAg                                                                                                                                    >  2:1328333/1‑92 (MQ=255)
                                                                           acaGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAAcc                                                                          <  2:240974/150‑1 (MQ=255)
                                                                                   aTGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCttt                                                                                <  1:1173236/136‑1 (MQ=255)
                                                                                   aTGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCttt                                                                                >  2:1173236/1‑136 (MQ=255)
                                                                                     ggAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACcgc                                                                <  2:62447/150‑1 (MQ=255)
                                                                                        aggCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAAGAGCGCCTCGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAg                                                             >  2:764506/3‑150 (MQ=255)
                                                                                           cACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACcg                                                                 <  1:1251892/143‑1 (MQ=255)
                                                                                           cACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACcg                                                                 >  2:1251892/1‑143 (MQ=255)
                                                                                           cACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACAAAGCAAAGGGACGGGGAGAAGTATGATCAGGACGGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCg                                                          >  2:691854/1‑150 (MQ=255)
                                                                                                      aTTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCg                                                                         <  1:240930/124‑1 (MQ=255)
                                                                                                      aTTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCg                                                                         >  2:240930/1‑124 (MQ=255)
                                                                                                       ttCAGCCAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGACCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGAt                                              >  2:1384244/1‑150 (MQ=255)
                                                                                                             aagcaGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGcc                                        >  1:1088588/1‑150 (MQ=255)
                                                                                                               gcaGAGAGAACGAGCACGAAACATCGCCAAGGGTGCCAGCACCCGGCCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGACTCACGTAGACGCTTTCAAAGCGGATAACTGCGAGCCGCTCACCCCGGGTTTCGCCAc                                      >  2:202207/1‑150 (MQ=255)
                                                                                                                   agagaACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACa                                  >  1:113201/1‑150 (MQ=255)
                                                                                                                   agagaACCAGCACGAAACATCGCCACGGGGGCCAGAACCGGAGGATGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCAGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACa                                  >  2:275766/1‑150 (MQ=255)
                                                                                                                         ccAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCaaa                                                                            >  1:691769/1‑102 (MQ=255)
                                                                                                                         ccAGCACGAAACATCGCCAAGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCaaa                                                                            <  2:691769/102‑1 (MQ=255)
                                                                                                                          cAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGt                                                                                                             >  1:191558/1‑68 (MQ=255)
                                                                                                                          cAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGt                                                                                                             <  2:191558/68‑1 (MQ=255)
                                                                                                                           aGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGc                                                            <  1:265783/116‑1 (MQ=255)
                                                                                                                           aGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGc                                                            >  2:265783/1‑116 (MQ=255)
                                                                                                                             cACGAAACATCGCCACGGGTGCCAGCACCCGACCaa                                                                                                                                          <  2:1369874/36‑1 (MQ=255)
                                                                                                                             cACGAAACATCGCCACGGGTGCCAGCACCCGACCaa                                                                                                                                          >  1:1369874/1‑36 (MQ=255)
                                                                                                                               cGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGc                                         >  1:623386/1‑131 (MQ=255)
                                                                                                                               cGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGc                                         <  2:623386/131‑1 (MQ=255)
                                                                                                                                gAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCg                     >  1:51808/1‑150 (MQ=255)
                                                                                                                                gAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCg                     >  1:306093/1‑150 (MQ=255)
                                                                                                                                 aaaCATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCgg                    >  2:508647/1‑150 (MQ=255)
                                                                                                                                         ccACGGGTGCCAGCACCCGACCAAGCATAGGGACGGTGTGACGAAACTTGCGTTAGCCAGGCGATCCTCACTGATACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCGGCGCCCAgc            >  1:739121/1‑150 (MQ=255)
                                                                                                                                          cACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCGGCGCCCAgcg           >  2:617140/1‑150 (MQ=255)
                                                                                                                                           tcGGGTGCTGGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGa                                                                       <  1:1018340/88‑1 (MQ=255)
                                                                                                                                           tcGGGTGCTGGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGa                                                                       >  2:1018340/2‑89 (MQ=255)
                                                                                                                                           tcGGGTGCTGGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCg                                                 <  2:1269886/110‑1 (MQ=255)
                                                                                                                                           tcGGGTGCTGGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCg                                                 >  1:1269886/2‑111 (MQ=255)
                                                                                                                                                   cAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGGAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCGGCGCCCAGCGCCAGCACaa  >  1:116956/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGCCACGGCCACCAACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCGGCGCCCAGCGCCAGCACAA  >  NC_002947/4564621‑4564917

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