Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,194,888 C→A 17.4% intergenic (+30/+23) 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,8880CA17.4% 93.7 / 11.0 44intergenic (+30/+23)PP_4574/PP_4575hypothetical protein/urea amidolyase‑like protein
Reads supporting (aligned to +/- strand):  ref base C (15/21);  new base A (5/3);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

AGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGC  >  NC_002947/5194746‑5195028
                                                                                                                                               |                                                                                                                                            
aGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCcgc                                                                                                                                        <  8:116479/149‑1 (MQ=255)
aGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCcgc                                                                                                                                        <  3:252966/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                                                                                                                                              >  5:136247/1‑133 (MQ=255)
          gTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTc                                                                                                                                              <  6:136247/133‑1 (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                                                                                           <  3:109412/149‑1 (MQ=255)
                                             tCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgct                                                                                           <  7:1016/149‑1 (MQ=255)
                                                aGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgctgc                                                                                         >  5:71022/1‑148 (MQ=255)
                                                       gATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgctgct                                                                                        <  7:136014/142‑1 (MQ=255)
                                                       gATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgctgct                                                                                        >  8:136014/1‑142 (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)
                                                                                                                                               |                                                                                                                                            
AGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTT‑CGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGC  >  NC_002947/5194746‑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: ‑