Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,288,746	C→T	9.2%	G126E (GGG→GAG)	PP_3760 ←	methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,288,746	0	C	T	9.2%	118.4 / 2.2	43	G126E (GGG→GAG)	PP_3760	methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (18/21);  new base T (2/2);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGGCGAAGTCACGGCGGCCACCGGCCAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAAC  >  NC_002947/4288601‑4288889                                                                                                                                                    |                                                                                                                                                cTCGGCGAAGTCACGGCGGCCACCGGCCAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCt                                                                                                                                               <  2:298849/148‑1 (MQ=255)    ggCGAAGTCACGGCGGCCACCGGCCAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTcc                                                                                                                                            <  7:225529/148‑1 (MQ=255)     gCGAAGTCACGGCGGCCACCGGCCAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCt                                                                                                                                           >  2:88711/1‑148 (MQ=255)       gAAGTCACGGCGGCCACCGGCCAGGCGGTAGTTTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGGGGCGTAGATGATGGTGCGTTCA‑AGCAGCTCTTCCTCa                                                                                                                                         >  1:201974/1‑149 (MQ=255)       gAAGTCACGGCGGCCACCGGCCAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCTCTTCCTCAc                                                                                                                                        >  1:31271/1‑149 (MQ=255)        aaGTCACGGCGGCCACCGGCCAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACg                                                                                                                                       <  5:252565/149‑1 (MQ=255)              cggcggCCACCGGCCAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAg                                                                                                                                 >  8:53553/1‑149 (MQ=255)                gcggcCACCGGCCAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCTCTTCCTCACGCAACAGGa                                                                                                                               <  2:201974/149‑1 (MQ=255)                gcggcCACCGGCCAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCTCTTCCTCACGCAACAGGa                                                                                                                               <  2:31271/149‑1 (MQ=255)                   gcCACCGGCCAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATgg                                                                                                                            <  4:291987/149‑1 (MQ=255)                    ccACCGGCCAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGc                                                                                                                           <  7:108239/149‑1 (MQ=255)                     cACCGGCCAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGc                                                                                                                           >  8:83123/1‑148 (MQ=255)                      aCCGGCCAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGAGGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCa                                                                                                                         >  7:293255/1‑149 (MQ=255)                          gCCAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATgg                                                                                                                      >  2:38036/1‑148 (MQ=255)                            cAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCGCGCAACAGGATGGCCATGGaa                                                                                                                    >  5:21797/1‑148 (MQ=255)                                 ggTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACAcc                                                                                                              >  7:278162/1‑149 (MQ=255)                                            cTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTccccc                                                                                                     <  4:284798/148‑1 (MQ=255)                                               gTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGtgc                                                                                                 <  7:53553/149‑1 (MQ=255)                                                  cTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCAGAGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGtgctg                                                                                               <  1:88711/149‑1 (MQ=255)                                                     gcgcATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAg                                                                                            <  8:278162/148‑1 (MQ=255)                                                          tGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCcgc                                                                                       >  7:59646/1‑148 (MQ=255)                                                                 cATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACAcc                                                                                                              >  3:15888/1‑118 (MQ=255)                                                                 cATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACAcc                                                                                                              <  4:15888/118‑1 (MQ=255)                                                                  aTCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGAc                                                                              >  7:235176/1‑149 (MQ=255)                                                                     gAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACttt                                                                           >  2:230824/1‑149 (MQ=255)                                                                       gTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCGCGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGa                                                                         <  6:21797/149‑1 (MQ=255)                                                                             gCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACg                                                                   >  1:159394/1‑149 (MQ=255)                                                                                      ggCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACGGCCAGGt                                                            <  1:38036/147‑1 (MQ=255)                                                                                      ggCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACGGCCAGGt                                                            <  7:83123/147‑1 (MQ=255)                                                                                       gCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGtgc                                                                                                 <  1:68263/109‑1 (MQ=255)                                                                                       gCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGtgc                                                                                                 >  2:68263/1‑109 (MQ=255)                                                                                               cAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACgg                                                                  <  5:74978/132‑1 (MQ=255)                                                                                               cAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACgg                                                                  >  6:74978/1‑132 (MQ=255)                                                                                                     aTGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACGGCCAGGTGCGCAGCAGCGGcacca                                           >  1:254702/1‑149 (MQ=255)                                                                                                            ttGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCAGCGATCCAGACTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGtt                                    <  8:59646/149‑1 (MQ=255)                                                                                                              gATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTg                                   <  2:159394/148‑1 (MQ=255)                                                                                                                    ggTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACgg                            <  1:230824/149‑1 (MQ=255)                                                                                                                         cGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAgg                       <  5:327689/149‑1 (MQ=255)                                                                                                                                       cGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTc         <  2:254702/149‑1 (MQ=255)                                                                                                                                        gTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCa        >  4:212050/1‑149 (MQ=255)                                                                                                                                           cA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACg      >  4:139070/1‑148 (MQ=255)                                                                                                                                            a‑aGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACgg     >  3:123374/1‑148 (MQ=255)                                                                                                                                               gcagcCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAAc  <  8:235176/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                CTCGGCGAAGTCACGGCGGCCACCGGCCAGGCGGTAG‑TTTTCCTCGTACTCGCGCATGCTCTGCATCGAGTAGATGCCCTGCTTGGCCCTGTCCAGCGAATGCGGGTTGATGTCGGTGGCGTAGATGATGGTGCGTTCA‑AGCAGCCCTTCCTCACGCAACAGGATGGCCATGGAATACACCTCCTCCCCCGTGCTGCAGCCCGCGATCCAGACTTTGATCGACGGCCAGGTGCGCAGCAGCGGCACCACTTCGTTGCGCACGGCAAGGAAGTGGTCAGGGTCACGGAAC  >  NC_002947/4288601‑4288889

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