Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,998,849	G→C	9.1%	Q47E (CAG→GAG) ‡	PP_2622 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,998,849	0	G	C	9.1%	156.3 / 2.7	51	Q47E (CAG→GAG) ‡	PP_2622	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/26);  new base C (3/2);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCACTGCGTTGAGCGTGCGTACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCAC  >  NC_002947/2998703‑2998983                                                                                                                                                    |                                                                                                                                         gggCACTGCGTTGAGCGTGCGTACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCTGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTgg                                                                                                                                         <  1:182639/148‑1 (MQ=255) gggCACTGCGTTGAGCGTGCGTACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTgg                                                                                                                                         <  7:174832/148‑1 (MQ=255)         cGTTGAGCGTGCGTACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCAcc                                                                                                                                  <  7:342107/147‑1 (MQ=255)           ttGAGCGTGCGTACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTg                                                                                                                              >  7:26175/1‑149 (MQ=255)              aGCGTGCGTACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑Gc                                                                                                                           <  2:331763/148‑1 (MQ=255)               gcgtgcgtACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTCAAGGGCTGCTCTGGGAGCACCCCTG‑GC‑g                                                                                                                         <  8:75477/149‑1 (MQ=255)                   gcgtACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCAt                                                                                                                    >  1:348275/1‑149 (MQ=255)                       aCCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTcag                                                                                                                >  1:294228/1‑149 (MQ=255)                         cTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTcag                                                                                                                >  2:97995/1‑147 (MQ=255)                         cTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTcag                                                                                                                >  2:173994/1‑147 (MQ=255)                                 gcgGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTc                                                                                                       >  1:358338/1‑148 (MQ=255)                                  cgGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGa                                                                                                     <  2:63699/149‑1 (MQ=255)                                   ggCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAg                                                                                                    >  5:119150/1‑149 (MQ=255)                                     cTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGt                                                                                                  <  1:9641/149‑1 (MQ=255)                                     cTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGt                                                                                                  <  1:97995/149‑1 (MQ=255)                                        gAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTgg                                                                                                <  2:294228/148‑1 (MQ=255)                                          aGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGgcgc                                                                                             <  5:3030/149‑1 (MQ=255)                                               tGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGAt                                                                                        >  3:183190/1‑149 (MQ=255)                                                 gAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATgg                                                                                      >  5:59070/1‑149 (MQ=255)                                                   gggTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATggcg                                                                                    <  1:177446/149‑1 (MQ=255)                                                     gTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATggcggc                                                                                  >  5:143074/1‑149 (MQ=255)                                                        cggcggATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGaa                                                                               <  6:308936/149‑1 (MQ=255)                                                        cggcggATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGaa                                                                               <  6:119150/149‑1 (MQ=255)                                                           cggATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGAGGGCGGCGAACtt                                                                            <  2:348275/149‑1 (MQ=255)                                                             gATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTtgt                                                                          <  6:220718/149‑1 (MQ=255)                                                                        tCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGTGCACCCCTG‑GC‑CGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGct                                                                >  7:121713/1‑148 (MQ=255)                                                                          aCCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCtctc                                                             <  6:59070/149‑1 (MQ=255)                                                                          aCCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCtctc                                                             <  5:365686/149‑1 (MQ=255)                                                                           ccTCCACAAACAAGTCGCCGTTGCTCAAGTCATGCCCATTGAAGTCGCGCAAGCCCAGCTC‑AGGGCTCCTCTCCCAGCACCCCTC‑GC‑GGGCATTCAGCAGCGTATCGAGGTGGCTCTAGATGGCGGCGAACTTGTGGCAACCGctctac                                                            >  4:274808/1‑147 (MQ=255)                                                                            cTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGACCTTGTGGGAAGCGCTCTccc                                                           >  6:44227/1‑149 (MQ=255)                                                                                acaAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCg                                                       <  1:173994/149‑1 (MQ=255)                                                                                  aagcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCggg                                                     <  6:143074/149‑1 (MQ=255)                                                                                   agcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGa                                                    <  5:354342/149‑1 (MQ=255)                                                                                    gcaagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAg                                                   <  7:251931/149‑1 (MQ=255)                                                                                      aagTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGc                                                 <  2:358338/149‑1 (MQ=255)                                                                                       agTCGCCGGTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCGGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGACGGGAgagc                                                <  7:286870/149‑4 (MQ=255)                                                                                           gCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCAtt                                                                                                                   <  2:58463/78‑1 (MQ=255)                                                                                           gCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCAtt                                                                                                                   >  1:58463/1‑78 (MQ=255)                                                                                                          aTGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGt                             <  2:302079/149‑1 (MQ=255)                                                                                                            gCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTc                            >  7:276932/1‑148 (MQ=255)                                                                                                                    aaGTCGCGCAGGCCCAGCGC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCa                    <  4:93861/148‑1 (MQ=255)                                                                                                                     aGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTcca                          >  4:183592/1‑141 (MQ=255)                                                                                                                     aGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTcca                          <  3:183592/141‑1 (MQ=255)                                                                                                                       tCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGc                 >  2:341147/1‑148 (MQ=255)                                                                                                                         gcgcAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGAt              >  2:356682/1‑149 (MQ=255)                                                                                                                               gCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCaaa          >  2:175563/1‑147 (MQ=255)                                                                                                                                   ccgcc‑AGGGGTGCTCCCAGAGCAGCCCTGAGCTGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGa                                                    >  2:271675/5‑103 (MQ=21)                                                                                                                                   ccgcc‑AGGGGTGCTCCCAGAGCAGCCCTGAGCTGGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGa                                                    <  1:271675/99‑1 (MQ=21)                                                                                                                                     cTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCAc  >  3:89845/1‑149 (MQ=255)                                                                                                                                         aGGGGTGCTCCCAGAGCAGCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCaaaaa        >  1:102287/1‑140 (MQ=255)                                                                                                                                         aGGGGTGCTCCCAGAGCAGCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCaaaaa        <  2:102287/140‑1 (MQ=255)                                                                                                                                           ggCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTtgt                                                                          >  2:234705/1‑72 (MQ=255)                                                                                                                                           ggCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTtgt                                                                          <  1:234705/72‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                         GGGCACTGCGTTGAGCGTGCGTACCTGTATGCGCGGCTCGAAGCGCTGGAGGGTGCGGCGGATATCCGCGCTCACCTGCACAAGCAAGTCGCCGTTGCTCACGTCATGCCCATTGAAGTCGCGCAGGCCCAGCTC‑AGGGCTGCTCTGGGAGCACCCCTG‑GC‑GGGCATTCAGCAGCGTTTCGAGGTGGCGCTTGATGGCGGCGAACTTGTGGGAAGCGCTCTCCCGTCGGGAGGCCTCTCGGCGGGTGCACTGGTCCACCATCAGGCGATCAAAAAGACCAC  >  NC_002947/2998703‑2998983

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