Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	290,240	C→T	10.2%	A53A (GCC→GCT)	PP_0235 →	peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	290,240	0	C	T	10.2%	103.2 / 2.1	39	A53A (GCC→GCT)	PP_0235	peroxidase
Reads supporting (aligned to +/- strand):  ref base C (17/18);  new base T (2/2);  total (19/20)
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

CTCGGCGATATCGCCCCCGATTTCGAACAGGATTCCAGCGAAGGCAAGATCCGCTTCCACGAGTGGCTGGGCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACATCAACGAAACACAGAATACCGTGGTCAACTTCCCGATCATTGCTGACGCTGACC  >  NC_002947/290094‑290385                                                                                                                                                   |                                                                                                                                                  cTCGGCGATATCGCCCCCGATTTCGAACAGGATTCCAGCGAAGGCAAGATCCGCTTCCACGAGTGGCTGGGCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCaa                                                                                                                                                 <  7:199298/149‑1 (MQ=255)         tatCGCCCCCGATTTCGAACAGGATTCCAGCGAAGGCAAGATCCGCTTCCACGAGTGGCTGTGCAACAGCTGGGGTGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAgg                                                                                                                                         <  1:315484/149‑1 (MQ=255)             gCCCCCGATTTCGAACAGGATTCCAGCGAAGGCAAGATCCGCTTCCACGAGTGGCTGGGCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGa                                                                                                                                     <  3:252403/149‑1 (MQ=255)              cccccGATTTCGAACAGGATTCCAGCGAAGGCAAGATCCGCTTCCACGAGTGGCTGGGCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAg                                                                                                                                    >  8:136087/1‑149 (MQ=255)               ccccGAATACGAACAGGATTCCAGCGAAGGCAAGATACGCTTCCACGAGTGCCTGGGCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACACCGGTGTGAACCACCGAGTAGGGACTGACCGCCAAGCTGAAGGACGAGt                                                                                                                                   >  4:120330/1‑149 (MQ=255)                  cGATTTCGAACAGGATTCCAGCGAAGGCAAGATCCGCTTCCACGAGTGGCTGGGCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTg                                                                                                                                >  2:34996/1‑149 (MQ=255)                          aaCAGGATTCCAGCGAAGGCAAGCTCCGCTTCCACGAGTGGCTGGGCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCCCCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCTAAGCTGAAGGACGAGTTTGCCAAgcgc                                                                                                                        >  5:258250/1‑149 (MQ=255)                              ggATTCCAGCGAAGGCAAGATCCGCTTCCACGAGTGGCTGGGCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCgggg                                                                                                                    <  5:260923/149‑1 (MQ=255)                                          aGGCAAGATCCGCTTCCACGAGTGGCTGGGCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGcc                                                                                                                              >  4:184059/1‑127 (MQ=255)                                          aGGCAAGATCCGCTTCCACGAGTGGCTGGGCAACAGCTGGGGGGTGTTGTTATCCCACCCCGCCGACTTCACCCCGGTGTGCACCCCCGAGTTTGGCCTGACCGCCAAGCGGAAGGCGCAGATTGCCAAGCGCGGGGTCAAGGCCATCg                                                                                                        >  6:240510/1‑149 (MQ=255)                                                             gAGTGGCTGGGCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGt                                                                                     <  3:57826/149‑1 (MQ=255)                                                              aGTGGCTGGGCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTc                                                                                    <  1:317756/149‑1 (MQ=255)                                                                  gCTGGGCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCTAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGAc                                                                                 <  6:258250/148‑1 (MQ=255)                                                                    tGGGCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGc                                                                                                       <  6:216397/124‑1 (MQ=255)                                                                    tGGGCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGc                                                                                                       >  5:216397/1‑124 (MQ=255)                                                                      ggCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCa                                                                            <  3:140382/149‑1 (MQ=255)                                                                       gCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCTAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCAt                                                                           >  8:55881/1‑149 (MQ=255)                                                                       gCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCAt                                                                           >  5:123239/1‑149 (MQ=255)                                                                        cAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATc                                                                          <  7:136087/149‑1 (MQ=255)                                                                                      ggttgttCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGATTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGAt                                                            <  5:240510/148‑1 (MQ=255)                                                                                        ttgttCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGa                                                          >  7:177703/1‑149 (MQ=255)                                                                                              tcCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACATCaa                                                    >  3:171319/1‑149 (MQ=255)                                                                                                ccACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCTATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACATCAAc                                                   <  1:34791/148‑1 (MQ=255)                                                                                                ccACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCTATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACATCAAc                                                   >  2:34791/1‑148 (MQ=255)                                                                                                    ccGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGAt                                                            >  3:231614/1‑135 (MQ=255)                                                                                                    ccGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGAt                                                            <  4:231614/135‑1 (MQ=255)                                                                                                          gACTTCACCGCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACATCAACGAAACACAGaa                                        <  3:120330/149‑1 (MQ=255)                                                                                                             ttaaCCCCGCTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAACGACGAGTTTGCAAACCGCGGGGTCAAGCCAATCGCCCTGTCGCTGGACCCGGACGACTCGCATCACAAGAAGAAACATGACATCAACCAAACACAGAATAc                                     >  4:59321/4‑149 (MQ=255)                                                                                                               cACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACATCAACGAAACACAGAATACCg                                   >  3:148210/1‑149 (MQ=255)                                                                                                                  cccGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACATCAACGAAACACAGAATACCGTgg                                >  6:45986/1‑149 (MQ=255)                                                                                                                    cGGTGTGCACCACCGAGTTGGGCCTGACCGCTAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACATCAACGAAACACAGAATACCGTGGTc                              <  7:55881/149‑1 (MQ=255)                                                                                                                           caccacCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACTTCAACGAAACACAGAATACCGTGGTCAACTTc                        <  8:177703/148‑1 (MQ=255)                                                                                                                           caccacCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACATCAACGAAACACAGAATACCGTGGTCAACTTcc                       <  4:171319/149‑1 (MQ=255)                                                                                                                            accaccGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACATCAACGAAACACAGAATACCGTGGTCAACTTccc                      >  7:5290/1‑149 (MQ=255)                                                                                                                              caccGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACATCAACGAAACACAGAATACCGTGGTCAACTTCCCGa                    <  8:5290/149‑1 (MQ=255)                                                                                                                               accGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACATCAACGAAACACAGAATACCGTGGTCAACTTCCCGa                    >  7:240591/1‑148 (MQ=255)                                                                                                                                       gggCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACATCAACGAAACACAGAATACCGTGGTCAACTTGCCGATCATTgct            <  8:240591/148‑1 (MQ=255)                                                                                                                                               ccGCCAAGCTGTAGTACGAGTTTGCCTTGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGAGCCGGTCGACTCGCATCAAAAGTGGATCGATGAGATCAACGAAACACAGAATACCGTTGTCAACTTCCCGATCATTgcttacgctgac   <  3:59321/149‑1 (MQ=255)                                                                                                                                                cGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACATCAACGAAACACAGAATACCGTGGTCAACTTCCCGATCATTGATGACGCTGAcc  >  3:189539/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CTCGGCGATATCGCCCCCGATTTCGAACAGGATTCCAGCGAAGGCAAGATCCGCTTCCACGAGTGGCTGGGCAACAGCTGGGGGGTGTTGTTCTCCCACCCGGCCGACTTCACCCCGGTGTGCACCACCGAGTTGGGCCTGACCGCCAAGCTGAAGGACGAGTTTGCCAAGCGCGGGGTCAAGGCCATCGCCCTGTCGGTGGACCCGGTCGACTCGCATCACAAGTGGATCGATGACATCAACGAAACACAGAATACCGTGGTCAACTTCCCGATCATTGCTGACGCTGACC  >  NC_002947/290094‑290385

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