Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,768,303	T→G	7.7%	L92R (CTG→CGG)	PP_5600 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,768,303	0	T	G	7.7%	132.1 / 3.8	52	L92R (CTG→CGG)	PP_5600	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/23);  new base G (2/2);  total (27/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.30e-01

CCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGC  >  NC_002947/3768159‑3768437                                                                                                                                                 |                                                                                                                                        ccgccTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTTGCCGTCTTTTTCGGCCTGCTGACCCTg                                                                                                                                        >  2:181653/1‑146 (MQ=255) ccgccTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTg                                                                                                                                        <  1:181653/146‑1 (MQ=255) ccgccTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCg                                                                                                                                    >  2:556165/1‑150 (MQ=255)    ccTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCg                                                                                                                                 <  2:736003/150‑1 (MQ=255)        gTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGa                                                                                                                             <  2:1088996/150‑1 (MQ=255)              ggCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCggg                                                                                                                       >  1:614391/1‑150 (MQ=255)              ggCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCggg                                                                                                                       >  1:712888/1‑150 (MQ=255)               gCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCGGCTCGCCGGCGAGCAgggt                                                                                                                      <  2:934204/150‑2 (MQ=255)                ccGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑Ca                                                                                                                    <  2:776645/150‑1 (MQ=255)                              aTGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGAACCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACccaac                                                                                                      <  1:365800/150‑1 (MQ=255)                                   ccAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCaaccaacc                                                                                                 >  1:362254/1‑150 (MQ=255)                                      aGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCcaaccaac                                                                                                  >  2:579677/1‑146 (MQ=255)                                      aGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCcaaccaac                                                                                                  <  1:579677/146‑1 (MQ=255)                                        gTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGc                                                                                            <  2:667637/150‑1 (MQ=255)                                             cAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGc                                                                                       >  1:703417/1‑150 (MQ=255)                                             cAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGc                                                                                       <  1:556165/150‑1 (MQ=255)                                                   ttCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGccc                                                                                 >  1:874710/1‑150 (MQ=255)                                                      aGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAag                                                                              >  1:632995/1‑150 (MQ=255)                                                            ccGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAgcgc                                                                        <  1:417063/150‑1 (MQ=255)                                                                     tCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTg                                                               <  2:632995/150‑1 (MQ=255)                                                                         ccGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCGGCTCGCCg                                                                                                                                 <  2:772544/81‑1 (MQ=255)                                                                         ccGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCGGCTCGCCg                                                                                                                                 >  1:772544/1‑81 (MQ=255)                                                                                 gACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCcgc                                                   <  1:169359/150‑1 (MQ=255)                                                                                   cTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGGGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCgcgc                                                 >  1:131856/1‑150 (MQ=255)                                                                                   cTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGATGAACCTGCTCGCCGGAGTGCTGGG‑CGTCCGCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCgcgc                                                 >  1:134868/1‑150 (MQ=255)                                                                                            cGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgc                                        <  1:254373/150‑1 (MQ=255)                                                                                             gCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgct                                       >  2:945802/1‑150 (MQ=255)                                                                                               cAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGTCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttc                                     >  1:452845/1‑150 (MQ=255)                                                                                                aGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttcg                                    <  1:830753/150‑1 (MQ=255)                                                                                                 gCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttcgc                                   >  2:291910/1‑150 (MQ=255)                                                                                                   cTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttcgctt                                 >  1:611244/1‑150 (MQ=255)                                                                                                      ttCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTtctc                              <  2:896711/150‑1 (MQ=255)                                                                                                         aGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTcgcg                           >  1:632815/1‑150 (MQ=255)                                                                                                            gcCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTc                        >  2:482285/1‑150 (MQ=255)                                                                                                             ccACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCtctc                                            <  2:721659/129‑1 (MQ=255)                                                                                                             ccACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCtctc                                            >  1:721659/1‑129 (MQ=255)                                                                                                             ccACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCt                       >  2:353910/1‑150 (MQ=255)                                                                                                             ccACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCt                       >  1:1063370/1‑150 (MQ=255)                                                                                                               acacTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGa                     <  1:291910/150‑1 (MQ=255)                                                                                                                  ctgctgGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACAt                  <  2:1063370/150‑1 (MQ=255)                                                                                                                   tgctgGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCACCGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATc                 >  1:457175/1‑150 (MQ=255)                                                                                                                     ctgGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGc               <  2:632815/150‑1 (MQ=255)                                                                                                                       ggCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGCG‑CAGCCCCTTACCCACACCACTTCAGACGCTCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCa             <  1:748531/150‑1 (MQ=255)                                                                                                                          cGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCGGGCCACGGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGt          >  1:12479/1‑150 (MQ=255)                                                                                                                                  tCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTgc  <  2:362254/150‑1 (MQ=255)                                                                                                                                            gACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGt    >  2:527047/1‑138 (MQ=255)                                                                                                                                            gACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGt    >  2:526979/1‑138 (MQ=255)                                                                                                                                            gACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGt    <  1:527047/138‑1 (MQ=255)                                                                                                                                            gACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGt    <  1:526979/138‑1 (MQ=255)                                                                                                                                             aCCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTccc                                          <  2:396731/99‑1 (MQ=255)                                                                                                                                             aCCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTccc                                          >  1:396731/1‑99 (MQ=255)                                                                                                                                              cccGGCTCGCCGGCGAGCAGGGTCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTtctc                              >  1:552600/1‑111 (MQ=255)                                                                                                                                                 |                                                                                                                                        CCGCCTGGTGCACGGCCGACAGCAGCGGCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGG‑CAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGC  >  NC_002947/3768159‑3768437

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