Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,707,330:1 +A 10.8% intergenic (+106/+36) PP_1501 → / ← PP_1502 hypothetical protein/OmpA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,707,3301.A10.8% 128.4 / 11.7 37intergenic (+106/+36)PP_1501/PP_1502hypothetical protein/OmpA family protein
Reads supporting (aligned to +/- strand):  ref base . (17/16);  new base A (2/2);  total (19/18)
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.90e-01

GGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCGCTAC  >  NC_002947/1707182‑1707479
                                                                                                                                                     |                                                                                                                                                      
ggTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑t                                                                                                                                                       >  2:1153087/1‑150 (MQ=255)
                  ggAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTcc                                                                                                                                     <  1:178611/150‑1 (MQ=255)
                       cacCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCa                                                                                                                                <  2:163007/150‑1 (MQ=255)
                                       gTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAgc                                                                                                                >  2:774642/1‑150 (MQ=255)
                                          aGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCt                                                                                                             <  2:1016120/150‑1 (MQ=255)
                                                 ggCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTc                                                                                                                                      >  2:832326/1‑118 (MQ=255)
                                                 ggCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTc                                                                                                                                      <  1:832326/118‑1 (MQ=255)
                                                  gcgcAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGA‑CGCGGGGCtt                                                                                                                                         <  2:662772/114‑1 (MQ=255)
                                                  gcgcAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGA‑CGCGGGGCtt                                                                                                                                         >  1:662772/1‑114 (MQ=255)
                                                  gcgcAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                                                                                     >  1:1096435/1‑150 (MQ=255)
                                                      agcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCa                                                                                                 >  1:1201099/1‑150 (MQ=255)
                                                           ttttGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAg                                                                                            <  2:556005/150‑1 (MQ=255)
                                                                      gcgAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCAc                                                                                >  1:681020/1‑150 (MQ=255)
                                                                          gaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTa                                                                                                                            <  1:423035/103‑1 (MQ=255)
                                                                          gaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTa                                                                                                                            >  2:423035/1‑103 (MQ=255)
                                                                            tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTg                                                                          <  1:684111/150‑1 (MQ=255)
                                                                            tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGc                                                                                                                                           <  1:82735/86‑1 (MQ=255)
                                                                            tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCCAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGc                                                                                                                                           >  2:82735/1‑86 (MQ=37)
                                                                              gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTcc                                                                                                                                     >  2:1126117/1‑90 (MQ=255)
                                                                              gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTcc                                                                                                                                     <  1:1126117/90‑1 (MQ=255)
                                                                                     gcgcgGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTCGtt                                                                 <  1:774642/150‑1 (MQ=255)
                                                                                      cgcgGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTCGTTg                                                                <  1:678550/150‑1 (MQ=255)
                                                                                                 cTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGA‑CGCGGGGc                                                                                                                                           >  1:249153/1‑65 (MQ=37)
                                                                                                 cTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGA‑CGCGGGGc                                                                                                                                           <  2:249153/65‑1 (MQ=37)
                                                                                                    gATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTCGTTGGAGATGGTCACCTc                                                  >  1:264657/1‑150 (MQ=255)
                                                                                                    gATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTCGTTGGAGATGGTCACCTc                                                  >  1:112594/1‑150 (MQ=255)
                                                                                                       ctcACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCAc                                               >  2:501997/1‑150 (MQ=255)
                                                                                                        tcACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCttt                                                                                                                                        <  2:976578/61‑2 (MQ=255)
                                                                                                        tcACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCttt                                                                                                                                        >  1:976578/1‑60 (MQ=255)
                                                                                                          acaGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTc                                                                                                                                      >  2:1172707/1‑61 (MQ=255)
                                                                                                          acaGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTc                                                                                                                                      <  1:1172707/61‑1 (MQ=255)
                                                                                                                  cacaACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGc                                                                                           >  1:780814/1‑96 (MQ=255)
                                                                                                                  cacaACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGc                                                                                           <  2:780814/96‑1 (MQ=255)
                                                                                                                     aacaTTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGACCCCAGGATTGTAGCTGAT‑GTGGTGATCACCTCCACCCGACGGTTCTgcg                                 >  1:1172823/1‑150 (MQ=255)
                                                                                                                        aTTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCAc                              <  1:1153087/150‑1 (MQ=255)
                                                                                                                            tCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTcc                          >  2:1120259/1‑150 (MQ=255)
                                                                                                                                      aaGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTc                <  2:681020/150‑1 (MQ=255)
                                                                                                                                              ccccGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTc        >  1:1187528/1‑150 (MQ=255)
                                                                                                                                                 cgcg‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTc                                                                    >  1:86163/1‑87 (MQ=255)
                                                                                                                                                 cgcg‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTc                                                                    <  2:86163/87‑1 (MQ=255)
                                                                                                                                                    g‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGACTCACCGAACGTGATG‑CCACAAGATGGTTGTCGAGGGAGAGGGTGTCGTGCACCGTACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCGCTAc  >  2:1041157/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCGCTAC  >  NC_002947/1707182‑1707479

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: ‑