Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,032,761 C→G 14.0% R94G (CGA→GGA)  PP_4435 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,032,7610CG14.0% 110.9 / 11.0 43R94G (CGA→GGA) PP_4435hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/17);  new base G (2/4);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

TTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCG  >  NC_002947/5032621‑5032890
                                                                                                                                              |                                                                                                                                 
ttGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCGCGCAGGAGCCtga                                                                                                                           <  4:241281/149‑2 (MQ=255)
  gACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGGTGGGTCG‑CGGAAGGTCAGGCTCGCGCAGGAGCCtga                                                                                                                           <  1:231468/148‑2 (MQ=255)
   aCAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGtt                                                                                                                        >  2:272606/1‑149 (MQ=255)
     aGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGtt                                                                                                                        >  2:235993/1‑147 (MQ=255)
       gCTTCTTAAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAgg                                                                                                                    >  2:203477/1‑149 (MQ=255)
           cttGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCGACGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAggc                                                                                                                   <  1:235993/147‑1 (MQ=255)
            ttGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCGCGCAGGAGCCtga                                                                                                                           >  4:198591/1‑136 (MQ=255)
            ttGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCGCGCAGGAGCCtga                                                                                                                           <  3:198591/137‑2 (MQ=255)
                               caAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTcaagccaa                                                                                                  >  4:55205/1‑143 (MQ=255)
                               caAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTcaagccaa                                                                                                  <  3:55205/143‑1 (MQ=255)
                               caAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGc                                                                                            >  1:291844/1‑149 (MQ=255)
                                      aaatTCCTCAATGTGTCTATCTATCTGATCGACATCACCACCGTGTTCGAGGCGATGAACGAGGTTTTGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGATCCTTCGCGAGCCTGTGTTGATGCGGCTATGTTCAAGCCAAGCCTGCTCGTTga                                                                                     <  7:151134/149‑1 (MQ=255)
                                      aaaaTCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGtt                                                                                       >  3:333508/1‑147 (MQ=255)
                                            ctgagattGACTATCTATCTTAACGACTTCTCGACCGACTTCGCTGCCATGAACGAGGTTTGTGAAGA‑GTGTGTCG‑CGGAAGGTCATGGTCCTGCGCGAGCGTGTGTTGAGGCGGCTATGTGCAAGCCAAGCCTGCTCGTTGAGATGt                                                                                <  2:146895/141‑1 (MQ=255)
                                            cTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAg                                                                                                                     <  8:53589/111‑1 (MQ=255)
                                            cTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAg                                                                                                                     >  7:53589/1‑111 (MQ=255)
                                            cTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTga                                                                                     <  5:348041/143‑1 (MQ=255)
                                            cTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTga                                                                                     >  6:348041/1‑143 (MQ=255)
                                                      ctatctatCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTagca                                                                     <  2:291844/149‑1 (MQ=255)
                                                            atctGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTga                                                                                     >  3:301139/1‑127 (MQ=255)
                                                            atctGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTga                                                                                     <  4:301139/127‑1 (MQ=255)
                                                             tctGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCg                                                              >  5:34355/1‑149 (MQ=255)
                                                               tGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGcc                                                            >  3:97849/1‑149 (MQ=255)
                                                                gAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCt                                                           >  3:330493/1‑149 (MQ=255)
                                                                 aaCGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCGCGCAGGAGCCtga                                                                                                                           >  2:251824/1‑83 (MQ=255)
                                                                 aaCGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCGCGCAGGAGCCtga                                                                                                                           <  1:251824/84‑2 (MQ=255)
                                                                 aaCGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGc                                                             <  6:34355/146‑1 (MQ=255)
                                                                 aaCGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCt                                                           <  4:330493/148‑1 (MQ=255)
                                                                           ccaccGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTaaaa                                                <  1:203477/149‑1 (MQ=255)
                                                                              ccGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAgcg                                             <  4:97849/149‑1 (MQ=255)
                                                                                   ttCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCAccc                                        >  8:111374/1‑149 (MQ=255)
                                                                                          ccATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCAGCTATGTTCAAGCCAAGCCTGCTAGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGgc                                 >  6:255953/1‑149 (MQ=255)
                                                                                              gAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGGCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGgcggc                              >  7:216482/1‑146 (MQ=255)
                                                                                                  gAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATc                          >  7:41568/1‑148 (MQ=255)
                                                                                                     gTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTc                                                                               <  5:142245/92‑1 (MQ=255)
                                                                                                     gTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTc                                                                               >  6:142245/1‑92 (MQ=255)
                                                                                                          ggAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTgg                  <  1:272606/148‑1 (MQ=255)
                                                                                                                       cg‑cgGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAg            >  2:320602/1‑142 (MQ=255)
                                                                                                                       cg‑cgGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAg            <  1:320602/142‑1 (MQ=255)
                                                                                                                          cgGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGgc   >  4:326792/1‑149 (MQ=255)
                                                                                                                           ggAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGgcg  <  4:333508/149‑1 (MQ=255)
                                                                                                                                       tCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCagag     <  2:3918/134‑1 (MQ=255)
                                                                                                                                       tCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCagag     >  1:3918/1‑134 (MQ=255)
                                                                                                                                           gcgcGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTgg                  >  2:158040/1‑117 (MQ=255)
                                                                                                                                           gcgcGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTgg                  <  1:158040/117‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                 
TTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGA‑GTGGGTCG‑CGGAAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCG  >  NC_002947/5032621‑5032890

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