Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	811,323	T→A	8.3%	L168M (TTG→ATG)	PP_5455 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	811,323	0	T	A	8.3%	133.0 / 2.1	46	L168M (TTG→ATG)	PP_5455	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/20);  new base A (2/2);  total (24/22)
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

AGTGCGCATACCCGCGCTTTCGCCGTCGATATACTGCTGCCGGCCTGGGCCGACCCCAATCCGCTGTGCCGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGAAGCGCTGGCGCTTGCCGCTGTCGGCC  >  NC_002947/811185‑811463                                                                                                                                           |                                                                                                                                             aGTGCGCATACCCGCGCTTTCGCCGTCGATATACTGCTGCCGGCCTGGGCCGACCCCAATCCGCTGTGCCGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATAccgcc                                                                                                                                   >  2:953203/1‑150 (MQ=255)       cATACCCGCGCTTTCGCCGTCGATATACTGCTGCCGGCCTGGGCCGACCCCAATCCGCTGTGCCGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCGGTATCCAAATGGATAcggcctgatgg                                                                                                                             <  1:1077267/150‑7 (MQ=255)                    tCGCCGTCGATATACTGCTGCCGGCCTGGGCCGACCCCAATCCGCTGTGCCGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGggtg                                                                                                                <  2:59950/150‑1 (MQ=255)                         gTCGATATACTGCTGCCGGCCTGGGCCGACCCCAATCCGCTGTGCCGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGAt                                                                                                           >  1:30552/1‑150 (MQ=255)                           cGATATACTGCTGCCGGCCTGGGCCCACCCCAATCCGCTGTGCCGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATgg                                                                                                         <  1:73252/150‑1 (MQ=255)                             atatAGTGCTGCCGGCCTGGGCCGACCCCAATCCGCTGTGCCGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCg                                                                                                       <  2:436873/150‑1 (MQ=255)                                   tgctgcCGGCCTGGGCCGACCCCAATCCGCTGTGCCGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCgg                                                                                                 >  2:153972/1‑150 (MQ=255)                                   tgctgcCGGCCTGGGCCGACCCCAATCCGCTGTGCCGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCgg                                                                                                 <  1:489440/150‑1 (MQ=255)                                   tgctgcCGGCCTGGGCCGACCCCAATCCGCTGTGCCGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCgg                                                                                                 <  2:590370/150‑1 (MQ=255)                                     ctgcCGGCCTGGGCCGACCCCAATCCGCTGTGCCGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCgggg                                                                                               <  2:1067691/150‑1 (MQ=255)                                                       cccAATCCGCTGTGCCGTTTCTGGCAGCAACACAGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAAc                                                                             >  1:1112313/1‑150 (MQ=255)                                                            tCCGCTGTGCCGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCaa                                                                        >  1:615113/1‑150 (MQ=255)                                                                cTGTGCCGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGa                                                                    <  2:424867/150‑1 (MQ=255)                                                                cTGTGCCGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCAGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGa                                                                    >  1:504313/1‑150 (MQ=255)                                                                   ggCCGTTGCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGcc                                                                 <  2:508519/149‑1 (MQ=255)                                                                    gCCGTTTCTGGCAGCAACACGGCCAAGCCTGGCACATGGGTCTTGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCa                                                                <  2:1112313/150‑1 (MQ=255)                                                                     caggtacTGGCAGAAACACGGCGAAGCCTGGAACACGTGGCTGGAGCGCTACCCATCAGGCCGTACCAATGTGGATACCGCCGGTTATCGCACCTGGGGTGGTGTTGGCGGTGCGGGGGCGTCGCGCTGTATCAACCTCAAATGAGGCAg                                                               >  2:518602/7‑150 (MQ=255)                                                                      cGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAgg                                                              <  2:1103695/150‑1 (MQ=255)                                                                               cagcaACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCCTCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGCGCGGGGGCGTCGTTCTGTATCACCCCCAAACGCGCCAGGCACTccctg                                                     <  2:14765/150‑3 (MQ=255)                                                                                 gcaACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCa                                                   >  2:318582/1‑150 (MQ=255)                                                                                    acacGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGGGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTc                                                >  2:37789/1‑150 (MQ=255)                                                                                       cgaccAAGCCTGGAACATGGGGTTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATGTCAGCTGGTGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTc                                             <  1:518602/147‑1 (MQ=255)                                                                                                cTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCa                                    >  1:31143/1‑150 (MQ=255)                                                                                                cTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGGTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTGCTGTATCAACCCCAAATGAGCCAGGCACTTCCGTCATTCGTCCATTATCCa                                    >  2:378440/1‑150 (MQ=255)                                                                                                     aCAGGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCgg                               <  1:966961/150‑1 (MQ=255)                                                                                                      cATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGc                              <  2:643867/150‑1 (MQ=255)                                                                                                         ggggCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGaa                           >  1:398613/1‑150 (MQ=255)                                                                                                           ggCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGCCATTCGTCCCTTAGCAACTCGGCGAAgc                         >  1:501278/1‑150 (MQ=255)                                                                                                              tGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGAAGCGCt                      >  2:690235/1‑150 (MQ=255)                                                                                                              tGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGAAGCGCt                      >  1:1052946/1‑150 (MQ=255)                                                                                                                 gCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCCAATGAGCCAGGCACTTCCTGCATTCGTCCATTATCCACTCGGCGAAGCGCTGgc                   <  1:423352/150‑1 (MQ=255)                                                                                                                        cccATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTggg                                                                                                                   <  1:590482/47‑1 (MQ=255)                                                                                                                        cccATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTggg                                                                                                                   >  2:590482/1‑47 (MQ=255)                                                                                                                        cccATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGAAGCGCTGGCGCTTGcc            >  1:497191/1‑150 (MQ=255)                                                                                                                         ccATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGAAGCGCTGGCGCTTGCCg           >  1:485306/1‑150 (MQ=255)                                                                                                                            tCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGAAGCGCTGGCGCTTGCCGCTg        >  1:969917/1‑150 (MQ=255)                                                                                                                             cAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGAAGCGCTGGCGCTTGCCGCTGt       <  1:153972/150‑1 (MQ=255)                                                                                                                              cggcggTATCCAAATGGATACGGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGAAGCGCTGGCGCTTGCCGCTGTc      >  2:885465/6‑150 (MQ=255)                                                                                                                              aGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGAAGCGCTGGCGCTTGCCGCTGTc      <  2:1052946/150‑1 (MQ=255)                                                                                                                               ggCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGAAGCGCTGGCGCTTGCCGCTGTCg     <  2:615113/150‑1 (MQ=255)                                                                                                                                gcggTATCCAAATGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGAAGCGCTGGCGCTTGCCGCTGTCgg    >  2:1050878/4‑150 (MQ=255)                                                                                                                                 ccGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGAAGCGCTGGCGCTTGCCGCTGTCGGc   >  2:827832/1‑150 (MQ=255)                                                                                                                                 ccGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAACGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGAAGCGCTGGCGCTTGCCGCTGTCGGc   >  2:270505/1‑150 (MQ=255)                                                                                                                                  ggTATCCAAATGGATACGGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGAAGCGCTGGCGCTTGCCGCTGTCGGcc  <  1:885465/149‑1 (MQ=255)                                                                                                                                         aTTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCa                                           <  1:607907/102‑1 (MQ=255)                                                                                                                                         aTTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCa                                           >  2:607907/1‑102 (MQ=255)                                                                                                                                           |                                                                                                                                             AGTGCGCATACCCGCGCTTTCGCCGTCGATATACTGCTGCCGGCCTGGGCCGACCCCAATCCGCTGTGCCGTTTCTGGCAGCAACACGGCCAAGCCTGGAACATGGGGCTGGGCCGCTACCCATCAGGCCGTATCCATTTGGATACCGCCGGTTATCGCACCTGGGGTGGTGATGGCGGTGCGGGGGCGTCGTTCTGTATCAACCCCAAATGAGCCAGGCACTTCCGGCATTCGTCCATTATCCACTCGGCGAAGCGCTGGCGCTTGCCGCTGTCGGCC  >  NC_002947/811185‑811463

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