Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,497,336:1	+G	8.5%	coding (261/765 nt)	PP_5562 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,497,336	1	.	G	8.5%	176.0 / 11.5	47	coding (261/765 nt)	PP_5562	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (20/23);  new base G (2/2);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

AAGGCGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGATCAGCGGTGCCACGGTGG  >  NC_002947/3497193‑3497484                                                                                                                                                 |                                                                                                                                                     aaGGCGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCTGTAACAGTTCTGGGGTGGTTCTGATCATGTTCAGCGGTAGTTCCAGCTCGGTGACC‑AGCTg                                                                                                                                                 <  4:290269/149‑1 (MQ=255)     cGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAAc                                                                                                                                             >  6:79756/1‑149 (MQ=255)      gTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACg                                                                                                                                            <  5:79756/149‑1 (MQ=255)         gAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCCGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGgcg                                                                                                                                         >  7:331532/1‑149 (MQ=255)              ccaccaCCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgcc                                                                                                                                    <  4:271248/149‑1 (MQ=255)                               gggCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTg                                                                                                                                                 <  6:105409/119‑1 (MQ=255)                               gggCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTg                                                                                                                                                 >  5:105409/1‑119 (MQ=255)                                                  ccTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTg                                                                                                 >  1:331440/1‑148 (MQ=255)                                                   cTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgccgagccga                                                                                                                             <  6:234539/119‑1 (MQ=255)                                                   cTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGgccgagccga                                                                                                                             >  5:234539/1‑119 (MQ=255)                                                        cGCTGTGTTCCAGGGTCTTGTCGTATTCCTGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTTGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTgg                                                                                          <  1:250802/149‑1 (MQ=255)                                                         gCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGc                                                                                         <  5:251415/149‑1 (MQ=255)                                                            gtgtTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTgg                                                                                                <  1:91027/139‑1 (MQ=255)                                                            gtgtTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTgg                                                                                                >  2:91027/1‑139 (MQ=255)                                                             tgtTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTg                                                                                     <  5:137036/149‑1 (MQ=255)                                                               ttCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTgg                                                                                    <  2:331440/148‑1 (MQ=255)                                                                     ggTCTTGTCGTATTCCACCGGGCACCGTTCGGGGGTGGGGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTcgc                                                                             <  7:218207/149‑1 (MQ=255)                                                                         ttGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGccc                                                                         >  8:68372/1‑149 (MQ=255)                                                                                    cAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCt                                                              >  4:150827/1‑149 (MQ=255)                                                                                        gggAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTgc                                                          <  3:150827/149‑1 (MQ=255)                                                                                         ggAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTgcg                                                         >  3:252768/1‑149 (MQ=255)                                                                                           aaCAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCgcgc                                                                           <  1:162709/129‑1 (MQ=255)                                                                                           aaCAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCgcgc                                                                           >  2:162709/1‑129 (MQ=255)                                                                                              aGTTCGGGGATGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTg                                                    >  3:253311/1‑149 (MQ=255)                                                                                                ttCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCa                                                  <  2:221902/149‑1 (MQ=255)                                                                                                 tCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAg                                                 <  7:68372/149‑1 (MQ=255)                                                                                                    ggggTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTa                                              >  4:213504/1‑149 (MQ=255)                                                                                                         ggtgCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGc                                         >  8:22106/1‑149 (MQ=255)                                                                                                         ggtgCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGc                                         <  3:213504/149‑1 (MQ=255)                                                                                                                     gTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTg                                                                                           <  4:143716/87‑1 (MQ=255)                                                                                                                     gTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTg                                                                                           >  3:143716/1‑87 (MQ=255)                                                                                                                          gCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcgga                        <  7:22106/149‑1 (MQ=255)                                                                                                                           cGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggac                       <  4:44647/149‑1 (MQ=255)                                                                                                                             gTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTc                                           >  1:208639/1‑127 (MQ=255)                                                                                                                             gTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTc                                           <  2:208639/127‑1 (MQ=255)                                                                                                                             gTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggatgg                     <  4:252768/149‑1 (MQ=255)                                                                                                                              tACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggacgga                    >  6:291051/1‑149 (MQ=255)                                                                                                                              tACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGc                           >  6:281826/1‑142 (MQ=255)                                                                                                                              tACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGc                           <  5:281826/142‑1 (MQ=255)                                                                                                                                gttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggacgg                     <  4:139333/145‑1 (MQ=255)                                                                                                                                gttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGcggacgg                     >  3:139333/2‑146 (MQ=255)                                                                                                                                 ttCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGAt                   >  6:286445/1‑147 (MQ=255)                                                                                                                                 ttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCt                                                                                                  <  4:196757/68‑1 (MQ=37)                                                                                                                                 ttCCAGCT‑GGTCACCGAGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCt                                                                                                  >  3:196757/1‑68 (MQ=37)                                                                                                                                          ggTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGATCAGCgg             <  5:217701/144‑1 (MQ=255)                                                                                                                                          ggTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGATCAGCgg             >  6:217701/1‑144 (MQ=255)                                                                                                                                             gACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGATCAGCGGTGCCACgg     >  5:299109/1‑149 (MQ=255)                                                                                                                                                c‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGATCAGCGGTGCCACGGTgg  >  6:272232/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     AAGGCGTGGAAGCCCACCACCCAGCGGGTCGGGCTGGTGATGCGGATGACCTGGCCGCTGTGTTCCAGGGTCTTGTCGTATTCCAGCGGGAACAGTTCGGGGGTGGTGCTGATCAGGTTCAGCGGTACTTCCAGCTCGGTGACC‑AGCTGGAACGGCGCGGCCGAGCCGATACGTTCGTAGAGTTCGACCAGGGCCTGGAAGTGGCCCTGGGTTTCGCGCCCGCTGCTGCGCTGTGCGCCTTGCAGGTATTCGCCAAACAGGGTTTGCGGACGGATCAGCGGTGCCACGGTGG  >  NC_002947/3497193‑3497484

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