Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,610,345 C→G 9.5% A262A (GCC→GCG PP_4081 → type IV/VI secretion system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,610,3450CG 9.5% 122.4 / 5.1 42A262A (GCC→GCGPP_4081type IV/VI secretion system protein
Reads supporting (aligned to +/- strand):  ref base C (19/19);  new base G (2/2);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

ATTCACCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGTTTTTGCCTCGCT  >  NC_002947/4610199‑4610491
                                                                                                                                                  |                                                                                                                                                  
aTTCACCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATg                                                                                                                                                 >  2:846470/1‑150 (MQ=255)
     ccAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCATACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGt                                                                                                                                            >  1:19699/1‑150 (MQ=255)
      cAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCCCCGCCATCGCCATGGCGGTg                                                                                                                                           <  2:565405/150‑1 (MQ=255)
                     atcCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTAc                                                                                                                            >  2:959443/1‑150 (MQ=255)
                     atcCGCGAACGGCGCGCCCCCACCCCCGAAGAAGTCTGCGCCCCCTACACCAACGTCGCCCCCAGAAAACTCGCTGTGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTAc                                                                                                                            <  1:351447/150‑1 (MQ=18)
                             aCTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATTCTCTTTTGGTGGAGCCCGTTGGTCATCTCCGCCATAGCCATGCGGCTGGCCCACGGGATTTACAGCTACCg                                                                                                                    >  1:192158/1‑150 (MQ=18)
                              cTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGc                                                                                                                   <  1:846470/150‑1 (MQ=255)
                                 cgcgGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTGCGCACCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTg                                                                                                                >  1:68430/1‑150 (MQ=255)
                                  gcgGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGc                                                                                                               >  1:368708/1‑150 (MQ=255)
                                    ggCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTCGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCAc                                                                                                             <  2:575792/150‑1 (MQ=255)
                                        ccACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGa                                                                                                         <  1:259982/150‑1 (MQ=255)
                                            ccccGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATcac                                                                                                     >  1:661059/1‑150 (MQ=255)
                                               cGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAAc                                                                                                  >  1:769864/1‑150 (MQ=255)
                                                   gaagTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAg                                                                                              <  2:985030/150‑1 (MQ=255)
                                                         tGCGACCCATACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTg                                                                                        <  2:19699/150‑1 (MQ=255)
                                                           cGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGa                                                                                      <  2:919233/150‑1 (MQ=255)
                                                                ccTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGc                                                                                 >  1:773140/1‑150 (MQ=255)
                                                                ccTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGc                                                                                 <  1:200904/150‑1 (MQ=255)
                                                                                ccccAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTgca                                                                 <  1:1022213/150‑1 (MQ=255)
                                                                                  ccAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCt                                                                                                                 <  2:145669/100‑1 (MQ=255)
                                                                                  ccAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCt                                                                                                                 >  1:145669/1‑100 (MQ=255)
                                                                                   cAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTgcagca                                                              >  1:148189/1‑150 (MQ=255)
                                                                                          ttCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCCTCCCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCa                                                       >  2:427845/1‑150 (MQ=255)
                                                                                                               ccGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAAc                                                                            >  2:470434/1‑108 (MQ=255)
                                                                                                               ccGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAAc                                                                            <  1:470434/108‑1 (MQ=255)
                                                                                                                cGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGGCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAg                                 >  1:979175/1‑150 (MQ=255)
                                                                                                                 gtggtgGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGt                                <  2:368708/150‑1 (MQ=255)
                                                                                                                    gtggAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCa                             >  2:220908/1‑150 (MQ=255)
                                                                                                                     tggAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCaa                            <  1:533696/150‑1 (MQ=255)
                                                                                                                      ggAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAAc                           <  1:696655/150‑1 (MQ=255)
                                                                                                                         gCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAAc                           <  1:229134/147‑1 (MQ=40)
                                                                                                                         gCCCGTTGGTCATCTCCGCCATAGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAAc                           >  2:229134/1‑147 (MQ=255)
                                                                                                                                ggTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACTGTGTCACAGCACCACTGAGTTCAACGAATTCGCGa                 <  2:595503/150‑1 (MQ=38)
                                                                                                                                  tCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAgg               >  1:626589/1‑150 (MQ=37)
                                                                                                                                     tctcCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCAt                                                                       <  1:807581/91‑1 (MQ=37)
                                                                                                                                     tctcCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCAt                                                                       >  2:807581/1‑91 (MQ=37)
                                                                                                                                     tctcCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGttt            >  1:459402/1‑150 (MQ=34)
                                                                                                                                     ccaccGCCATGGCGATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGttt            <  1:392189/147‑1 (MQ=12)
                                                                                                                                      caccGCCATGGCGATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAAt                                                            <  2:433014/99‑1 (MQ=255)
                                                                                                                                      caccGCCATGGCGATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAAt                                                            >  1:433014/3‑101 (MQ=255)
                                                                                                                                      caccGCCATGGCGATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGtttt           >  2:428798/3‑150 (MQ=12)
                                                                                                                                       tcCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGCAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGttttt          <  2:769864/150‑1 (MQ=17)
                                                                                                                                               cgccatGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGTTTTTGCCTCGCt  <  2:661059/150‑1 (MQ=32)
                                                                                                                                                  |                                                                                                                                                  
ATTCACCAACTGCACGGGATCATCCGCGAACTGCGCGGCCCCACCCCCGAAGAAGTCTGCGACCCCTACACCAACGTCGCCCCCAGAAACTTCCGCATGAGCCGGATATGGCCGTGGTGGAGCCCGTTGGTCATCTCCGCCATCGCCATGGCGGTGGCCTACGGGATTTACAGCTACCGCCTGCACCTGATCACAACCGAGGTGCTGGAGTCGCTGAACGGCATTTTGCAGCAATAGCCAACCGTGTCACAGCACCACTGAGTTCAACGAATTCGCGAGGTTTTTGCCTCGCT  >  NC_002947/4610199‑4610491

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: