Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,194,887 T→G 19.5% intergenic (+29/+24) PP_4574 → / ← PP_4575 hypothetical protein/urea amidolyase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,194,8870TG19.5% 82.3 / 17.6 46intergenic (+29/+24)PP_4574/PP_4575hypothetical protein/urea amidolyase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (15/21);  new base G (5/4);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CCCAAGCAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGC  >  NC_002947/5194739‑5195028
                                                                                                                                                     |                                                                                                                                             
cccAAGCAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCt                                                                                                                                               >  4:109412/1‑149 (MQ=255)
       aGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCcgc                                                                                                                                        <  3:252966/149‑1 (MQ=255)
       aGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCcgc                                                                                                                                        <  8:116479/149‑1 (MQ=255)
           cTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGaaa                                                                                                                                    <  1:112824/149‑1 (MQ=255)
           cTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCGAGCCCGCGaaa                                                                                                                                    <  5:192592/149‑1 (MQ=255)
              ggAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAgg                                                                                                                                  <  6:224985/148‑1 (MQ=255)
              ggAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAg                                                                                                                                   <  5:184746/147‑1 (MQ=255)
              ggAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCGAGCCCGCGAAAg                                                                                                                                   <  3:291888/147‑1 (MQ=255)
               gAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGc                                                                                                                                >  1:60057/1‑149 (MQ=255)
                 gTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTc                                                                                                                                              <  6:136247/133‑1 (MQ=255)
                 gTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTc                                                                                                                                              >  5:136247/1‑133 (MQ=255)
                   acacCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCGAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCcag                                                                                                                             >  7:85902/1‑148 (MQ=255)
                             gACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTAc                                                                                                                    <  1:299505/147‑1 (MQ=255)
                                        tGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGcccc                                                                                                        <  2:60057/148‑1 (MQ=255)
                                         gcgcGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAg                                                                                                      >  3:90368/1‑149 (MQ=255)
                                          cgcgACCTGATCGAGCACAAGATGTAGCCGCTGACCCTGCATGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGGTGCCTGTCCCGGCCCTTT‑CGCGGGCTCGCCCGCTAATGGGCCAGCAGTCTTACGACCAACGCCCCAgg                                                                                                     <  2:147504/149‑1 (MQ=255)
                                           gcgACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAgg                                                                                                     >  8:1016/1‑148 (MQ=255)
                                               ccTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAg                                                                                                <  4:187280/149‑1 (MQ=255)
                                                  gATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctg                                                                                             >  1:255288/1‑149 (MQ=255)
                                                    tCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgct                                                                                           <  7:1016/149‑1 (MQ=255)
                                                    tCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgct                                                                                           <  3:109412/149‑1 (MQ=255)
                                                       aGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgctgc                                                                                         >  5:71022/1‑148 (MQ=255)
                                                              gATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgctgct                                                                                        >  8:136014/1‑142 (MQ=255)
                                                              gATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgctgct                                                                                        <  7:136014/142‑1 (MQ=255)
                                                              gATGTATACGCTGACCCTGAAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTTCCTGTACCGGCCCTTT‑CGAGGGCACGCCCGAGAAAGGGCAAGCAGTATTAAAACCAACTCCCCAGTTAAAGCTGATGATCTCgaa                                                                                  >  4:65788/1‑146 (MQ=255)
                                                                           caccTGCAAGGGCTGGAGGACGGGATCCGCTACCAGTCCAAGTGAGGTTCCCTGAACCGGCCCTTT‑CCCGGGCTCGCCCGCGAAAGGGCCAGCAGTCATACGACCAACGCCCCAGCTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGc                                                                    <  8:85902/147‑1 (MQ=255)
                                                                           aCCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTAATACGACCAACGCCCCAGGTAAAGCTGGTGCTATCGCCAGGCCGCGTccag                                                                     >  6:183470/1‑146 (MQ=255)
                                                                                gCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCCGCAGTCTTTCGCCCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCCGGCCTCGTCCTGcacccc                                                               >  3:165463/1‑147 (MQ=255)
                                                                                    gggTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGcc                                                           <  7:38071/149‑1 (MQ=255)
                                                                                                 ggATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGcacca                                              >  8:168056/1‑149 (MQ=255)
                                                                                                      cgcTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTc                                         <  6:71022/149‑1 (MQ=255)
                                                                                                            cTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGc                                   <  2:255288/149‑1 (MQ=255)
                                                                                                                 cAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCa                              >  1:200632/1‑149 (MQ=255)
                                                                                                                         gTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTgcgc                        >  7:13707/1‑147 (MQ=255)
                                                                                                                              cTGTACCGGCGCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACGGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgc                 <  5:183470/149‑1 (MQ=255)
                                                                                                                               actgctggCCCTTT‑CGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcg                >  6:121616/7‑149 (MQ=255)
                                                                                                                                 gtccGGCCCTTTGCGGGGGCGCCCCCGCTAACGGGCCGGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgc                 <  6:202117/145‑1 (MQ=255)
                                                                                                                                   ctggCCCTTT‑CGCGGGCGAGCCCGCGAAAGGGCCGGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgc                 >  5:202117/3‑144 (MQ=255)
                                                                                                                                     ggCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccagcg          >  3:246243/1‑149 (MQ=255)
                                                                                                                                      gCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccagcgc         >  6:176597/1‑149 (MQ=255)
                                                                                                                                        ccTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccagcgcca       <  4:90368/149‑1 (MQ=255)
                                                                                                                                        ccTTT‑CGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGt                    >  6:24923/1‑136 (MQ=255)
                                                                                                                                        ccTTT‑CGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGt                    <  5:24923/136‑1 (MQ=255)
                                                                                                                                        ccTTT‑CGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCACCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccagcgc         >  7:258072/1‑147 (MQ=255)
                                                                                                                                          ttt‑CGCGGGCGAGCCCGCGAAAGGGCCGGCAGTCTTACGACCAACGCCCCCGGTAAAGCTGCTGCACTCGCCCGGCCTCGTCCTGCACCACCGCCCTGAACCGCCCCACCGTCCCCGGCATGCACTGCGCCAGTTGCgccagcgccagc     >  1:231301/1‑149 (MQ=255)
                                                                                                                                            t‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCgccagcgccagc     <  2:200632/147‑1 (MQ=255)
                                                                                                                                               gcgGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAgcggc  <  4:165463/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
CCCAAGCAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGC  >  NC_002947/5194739‑5195028

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