Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,797,359	G→A	9.1%	E138K (GAA→AAA)	PP_2449 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,797,359	0	G	A	9.1%	121.4 / 2.2	44	E138K (GAA→AAA)	PP_2449	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/21);  new base A (2/2);  total (21/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 = 4.10e-01

CTGCGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATCCAGCACG  >  NC_002947/2797212‑2797503                                                                                                                                                        |                                                                                                                                                 cTGCGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCg                                                                                                                                                  >  5:92786/1‑148 (MQ=255)      tATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAgg                                                                                                                                            <  3:190469/149‑1 (MQ=255)      tATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAgg                                                                                                                                            <  3:291332/149‑1 (MQ=255)               gCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCct                                                                                                                                   >  2:220523/1‑149 (MQ=255)                   tGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCAAAAAGGGCATCCCCTCTTc                                                                                                                               >  1:192037/1‑149 (MQ=255)                   tGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCAAAAAGGGCATCCCCTCTTc                                                                                                                               >  1:196958/1‑149 (MQ=255)                           aGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGt                                                                                                                       >  7:156850/1‑149 (MQ=255)                                     cAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACgg                                                                                                             >  7:38260/1‑149 (MQ=255)                                            gTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCg                                                                                                      <  8:156850/149‑1 (MQ=255)                                             tGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCTACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCgg                                                                                                     >  5:171063/1‑149 (MQ=255)                                              gCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCggt                                                                                                    >  3:28322/1‑149 (MQ=255)                                                    gcgcTGGCTGACAAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGAACTCAGCTATATGACGGCGCGCTTGATCGAAAAGGGCATCCCCCCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGtgc                                                                                               <  1:220523/148‑1 (MQ=255)                                                      gcTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGtgctgc                                                                                            >  3:81429/1‑149 (MQ=255)                                                       cTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCggtg                                                                                                   <  1:129745/141‑1 (MQ=255)                                                        tGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGctgctg                                                                                          <  6:92786/149‑1 (MQ=255)                                                                 aGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTc                                                                                 <  4:28322/149‑1 (MQ=255)                                                                   ccGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTccc                                                                               <  8:38260/149‑1 (MQ=255)                                                                   ccGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCAAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTccc                                                                               <  2:192037/149‑1 (MQ=255)                                                                   ccGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCAAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTccc                                                                               <  2:196958/149‑1 (MQ=255)                                                                      gcgcCCCATGTGTGGGATGTGGAATTTGTAAATGGCGCCAACACTCGCGCTCGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTcc                                                                                <  1:274595/149‑1 (MQ=255)                                                                          cccAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAggg                                                                                                                                           <  3:307548/82‑1 (MQ=255)                                                                          cccAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAggg                                                                                                                                           >  4:307548/1‑82 (MQ=255)                                                                          cccAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCAccc                                                                          >  6:20733/1‑147 (MQ=255)                                                                          cccAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCAccc                                                                          <  5:20733/147‑1 (MQ=255)                                                                            cAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTg                                                                       >  3:68/1‑148 (MQ=255)                                                                               gggTGGGATGTGGAAGTGGTACAGGGCGCCTACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGTGGCCGTTCCCCACCCGTgaaga                                                                   <  6:171063/149‑1 (MQ=255)                                                                                       tgtgGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAgg                                                           >  3:154060/1‑149 (MQ=255)                                                                                          ggAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGgcgc                                                        >  3:100062/1‑149 (MQ=255)                                                                                           gAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTc                                                                                 >  3:77134/1‑123 (MQ=255)                                                                                           gAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTc                                                                                 <  4:77134/123‑1 (MQ=255)                                                                                                 gTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAggg                                                  <  3:203455/148‑1 (MQ=255)                                                                                                   aCAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTc                                                                                 <  1:302997/115‑1 (MQ=255)                                                                                                   aCAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTc                                                                                 >  2:302997/1‑115 (MQ=255)                                                                                                     aGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGaa                                             <  4:154060/149‑1 (MQ=255)                                                                                                       ggCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTAGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGAAGGCCCGCAGGGACGaact                                           >  4:153469/1‑147 (MQ=255)                                                                                                         cgcCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGa                                         <  4:81429/149‑1 (MQ=255)                                                                                                           ccAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATg                                       <  5:43708/149‑1 (MQ=255)                                                                                                                a‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGcc                                    >  2:106353/1‑147 (MQ=255)                                                                                                                  t‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAg                                 >  5:35302/1‑149 (MQ=255)                                                                                                                      g‑‑cgcACTAAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATACCCTATTCGGTGGAGTACCGTGACGGAAAGCCGGTCGTGATGCTGGGGCAGTTCACCACCCGAGAAGAACCGCAGGCGCGCAGGCACGAAGTGATGGCCAAGATg                              >  2:113706/1‑149 (MQ=255)                                                                                                                             cTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTAc                         <  6:35302/149‑1 (MQ=255)                                                                                                                              tCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTAcc                        >  2:14214/1‑149 (MQ=255)                                                                                                                                aGCGATATGAGGGGGCGGCTGATCGAATAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGGTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGGGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGa                      <  3:153469/149‑1 (MQ=255)                                                                                                                                       tGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCggtg               <  4:68/149‑1 (MQ=255)                                                                                                                                                    gATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATCCAGCACg  >  4:262689/1‑149 (MQ=255)                                                                                                                                                        |                                                                                                                                                 CTGCGTATCCGTGAGCTCTGCGGTGGAGCCGACCTGCAGCGCCGTGCCGAGGCGCTGGCAGACCAGCCGGCGCCCCAGGGGTGGGATGTGGAAGTGGTACAGGGCGCCAACA‑T‑GGG‑‑CGCACTCAGCGATATGACGGCGCGCCTGATCGAAAAGGGCATCCCCTCTTCGGTGGAGTACCGTGACGGCAAGCCGGTGGTGCTGCTGGGGCCGTTCCCCACCCGTGAAGAAGCGCAGGCGCGCAGGGACGAAGTGATGGCCAAGCTGGGTACCGATTCGGTGGTCATCCAGCACG  >  NC_002947/2797212‑2797503

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