Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,153,770	G→A	8.1%	G271G (GGG→GGA)	PP_1908 →	putative Peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,153,770	0	G	A	8.1%	139.1 / 2.6	49	G271G (GGG→GGA)	PP_1908	putative Peptidase
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base A (2/2);  total (24/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 = 1.00e+00

CCGTTTCTGGCAGGGTGTGCTGGATACCACGCACCGGCAGTTCATCGCCATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGTTCAGGAGTCGCCTTTCGATCGCTTCTCCAAGCGCGTGGGTGCCAGCGTGGCTGAGCACCTGGCGATGT  >  NC_002947/2153623‑2153915                                                                                                                                                     |                                                                                                                                                  ccGTTTCTGGCAGGGTGTGCTGGATACCACGCACCGGCAGTTCATCGCCATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGc                                                                                                                                                  >  1:207406/1‑149 (MQ=255) ccGTTTCTGGCAGGGTGTGCTGGATACCACGCACCGGCAGTTCATCGCCATGGACAAGAAGGGGCGGGGTGCGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGCCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGc                                                                                                                                                  >  1:100640/1‑149 (MQ=255)      tCTGGCAGGGTGTGCTGGATACCACGCACCGGCAGTTCATCGCCATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCtggt                                                                                                                                              <  5:45153/148‑1 (MQ=255)      tCTGGCAGGGTGTGAGGGATACCACGCACCGGCAGTGCATCGCCATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCtggt                                                                                                                                              <  2:323672/148‑1 (MQ=255)          gCAGGGTGTGCTGGATACCACGGACCGGCAGTTCATCGCCATGGTCAAGCAGGGGCGGGGTTAGCGGCTGATGGATAAGGAGCACGCGGAGTTGTTCAGCTGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCtggtgg                                                                                                                                            <  7:105184/146‑1 (MQ=255)             gggTGTGCTGGATACCACGCACCGGCAGTTCATCGCCATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCt                                                                                                                                      <  1:25154/149‑1 (MQ=255)                tgtgCTGGATACCACGCACCGGCAGTTCATCGCCATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAGAGCGCTGGGGCTGGTGGATGGCTTg                                                                                                                                    <  7:240877/149‑1 (MQ=255)                  tgCTGGATACCACGCACCGGCAGTTCATCGCCATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTggg                                                                                                                                  >  2:175062/1‑148 (MQ=255)                        aTACCACGCACCGGCAGTTCATCGACATGGTCAAGCAGGTGCTGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGcc                                                                                                                           <  2:89977/149‑1 (MQ=255)                        aTACCACGCACCGGCAGTTCATCGACATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAATGATACGGAGCACCCGGAGTTGTTCAGCGGCATGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGcc                                                                                                                           <  2:84483/149‑1 (MQ=255)                            cACGCACCGGCAGTTCATCGCCATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGc                                                                                                                        >  2:141707/1‑148 (MQ=255)                             aCGCACCGGCAGTTCATCGCCATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTa                                                                                                                      >  2:8940/1‑149 (MQ=255)                              cGCACCGGCAGTTCATCGCCATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTAt                                                                                                                     >  2:97308/1‑149 (MQ=255)                               gCACCGGCAGTTCATCGCCATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGTCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTAtg                                                                                                                    >  7:87277/1‑149 (MQ=255)                                 aCCGGCAGTTCATCGCCATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGACTGGTGGATGGCTTGGGCAGTGCCAGCTAtgtg                                                                                                                  >  5:375051/1‑149 (MQ=255)                                 aCCGGCAGTTCATCGCCATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGACTGGTGGATGGCTTGGGCAGTGCCAGCTAtgtg                                                                                                                  >  5:380298/1‑149 (MQ=255)                                       aGTTCATCGCCATGTTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGGAGGCGAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATTTGGCGCGt                                                                                                            <  2:100640/149‑1 (MQ=255)                                       aGTTCATCGCCATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGGGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGt                                                                                                            <  5:197443/149‑1 (MQ=255)                                                 cATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGACTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTgg                                                                                                  <  6:375051/149‑1 (MQ=255)                                                 cATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGACTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGAGTGAGATCGTgg                                                                                                  <  6:380298/149‑1 (MQ=255)                                                         aGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCACCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAg                                                                                          <  2:151827/149‑1 (MQ=255)                                                          gCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAg                                                                                          >  5:155326/1‑148 (MQ=255)                                                             ggggCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAgt                                                                                      <  8:350079/149‑1 (MQ=255)                                                              gggCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAgtt                                                                                     >  7:327165/1‑149 (MQ=255)                                                                gCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGTCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAgttgg                                                                                   <  8:87277/149‑1 (MQ=255)                                                                gCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAgttgg                                                                                   <  1:141707/149‑1 (MQ=255)                                                                 cGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAgttggt                                                                                  <  5:246330/149‑1 (MQ=255)                                                                    ggTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGa                                                                               <  1:175062/149‑1 (MQ=255)                                                                         gcgcCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCa                                                                          >  3:132782/1‑149 (MQ=255)                                                                             cTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGCCTTCACCGt                                                                      >  1:233570/1‑149 (MQ=255)                                                                                   gATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGt                                                                      >  8:62140/1‑143 (MQ=255)                                                                                   gATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGt                                                                      <  7:62140/143‑1 (MQ=255)                                                                                    aTAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGTTCAGGAg                                                               >  6:349369/1‑149 (MQ=255)                                                                                       aGGAGCACCCTGAGTTGTTCAGCGTCCTGTTCTTTTCGGGTGATCAGGCCAAG‑GCTGGGGGGCTGGTTGATGGCTTTGGCAGTGCCAGCTATGTGGCGCGTGTGATCGTGTTCGAGAAGGAGTTTGTTGACTTTACTGTTCAGGAGTc                                                             <  8:242385/148‑1 (MQ=255)                                                                                        ggAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGTTCAGGAGTCGc                                                           >  5:284589/1‑149 (MQ=255)                                                                                        ggAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGGCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGTTCAGGAGTCGc                                                           <  1:242446/149‑1 (MQ=255)                                                                                            cACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGTTCAGGAGTCGCCttt                                                       <  2:207406/149‑1 (MQ=255)                                                                                              cccGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGTTCAGGAGTCGCCTTTc                                                      >  3:177926/1‑148 (MQ=255)                                                                                               ccGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAg                                                                                          <  6:169312/111‑1 (MQ=255)                                                                                               ccGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAg                                                                                          >  5:169312/1‑111 (MQ=255)                                                                                                           aGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGTTCAGGAGTCGCCTTTCGATCGCTTCTcc                                          >  1:185063/1‑147 (MQ=255)                                                                                                                ccTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTg                                                                                                           <  3:267379/77‑1 (MQ=255)                                                                                                                ccTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTg                                                                                                           >  4:267379/1‑77 (MQ=255)                                                                                                                     tCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGTTCAGGAGTCGCCTTTCGATCGCTTCTCCAAGCGCGTggg                               <  6:155326/148‑1 (MQ=255)                                                                                                                         gTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGTTCAGGAGTCGCCTTTCGATCGCTTCTCCAAGCGCGTGGGt                              <  4:177926/145‑1 (MQ=255)                                                                                                                                 aGCAGGCCAAG‑GCGCTGGGGCTGGTGGAGGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGTTCAGGAGTCGCCGTTCGATCGCTTCTCCAAGCGCGTGTGTGCCAGCGTGGCt                  >  6:275584/1‑149 (MQ=255)                                                                                                                                         aaG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGTTCAGGAGTCGCCTTTCGATCGCTTCTCCAAGCGCGTGGGTGCCAGCGTGGCTGAGCACCt          >  7:4722/1‑149 (MQ=255)                                                                                                                                               gcTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGTTCAGGAGTCGCCTTTCGATCGCTTCTCCAAGCGCGTGGGTGCCAGCGTGGCTGAGCACCTGGCGa     >  5:116271/1‑149 (MQ=255)                                                                                                                                               gcTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGTTAAGGAGTCGCCTTTCGATCGCTTCTCCAAGCGCGTGGGTGCCCGCGTGGCTGAGGACCTGGCGa     >  7:96982/1‑149 (MQ=255)                                                                                                                                                   gggCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGTTCAGGAGTCGCCTTTCGATCGCTTCTCCAAGCGCGTGGGTGCCAGCGTGGCTGAGCACCTGGCGAtgt  <  1:8940/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CCGTTTCTGGCAGGGTGTGCTGGATACCACGCACCGGCAGTTCATCGCCATGGTCAAGCAGGGGCGGGGTGAGCGCCTGAAGGATAAGGAGCACCCGGAGTTGTTCAGCGGCCTGATCTGGTCGGGTGAGCAGGCCAAG‑GCGCTGGGGCTGGTGGATGGCTTGGGCAGTGCCAGCTATGTGGCGCGTGAGATCGTGGGCGAGAAGGAGTTGGTTGACTTCACCGTTCAGGAGTCGCCTTTCGATCGCTTCTCCAAGCGCGTGGGTGCCAGCGTGGCTGAGCACCTGGCGATGT  >  NC_002947/2153623‑2153915

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