Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,707,335 Δ1 bp 9.5% intergenic (+111/+31) 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,3350T. 9.5% 126.6 / 7.9 42intergenic (+111/+31)PP_1501/PP_1502hypothetical protein/OmpA family protein
Reads supporting (aligned to +/- strand):  ref base T (19/19);  new base . (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 = 7.82e-01

GGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCGCTACCGG  >  NC_002947/1707200‑1707482
                                                                                                                                        |                                                                                                                                                    
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                                                                                                                                            >  1:662772/1‑114 (MQ=255)
                                gcgcAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGA‑CGCGGGGCtt                                                                                                                                            <  2:662772/114‑1 (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)
                                                                                                                                      gATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCGCTACCgg  <  2:112594/150‑1 (MQ=255)
                                                                                                                                      gATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCGCTACCgg  <  1:501997/150‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                    
GGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCG‑TCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTG‑CCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCGCTACCGG  >  NC_002947/1707200‑1707482

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