Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	990,700	T→A	7.1%	L193Q (CTG→CAG)	PP_0856 →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	990,700	0	T	A	7.1%	153.3 / 2.5	55	L193Q (CTG→CAG)	PP_0856	lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (23/28);  new base A (2/2);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

AACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCTTGTTGCTGTCTG  >  NC_002947/990552‑990846                                                                                                                                                     |                                                                                                                                                    aacaacGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTg                                                                                                                                                    >  1:288736/1‑150 (MQ=255) aacaacAGTGACGGGGGGGTCGTGCAGCTCCCGCACGATCGAATGATCGGCATGGAGGCCGCCACTGGCGACCGCCGTTGGAGCTACGAAAACTCCCCGGCCGTGCTGATGGTGCGTGGCACCGGTGCGCCGATCGCCCACAAACGCCTg                                                                                                                                                    >  2:472760/1‑150 (MQ=255)  acaacGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAATCTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTgg                                                                                                                                                   <  1:358388/150‑1 (MQ=255)  acaacGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTgg                                                                                                                                                   <  2:714701/150‑1 (MQ=255)    aacGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCACCACCAACCGCCTGGCg                                                                                                                                                 >  1:188951/1‑150 (MQ=255)         tGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGc                                                                                                                                            <  1:173396/150‑1 (MQ=255)          gACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCAGGCGGTTggt                                                                                                                                           <  2:201623/150‑3 (MQ=255)              tggtggtggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGcc                                                                                                                                       >  2:751447/1‑150 (MQ=255)               ggtggtggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCt                                                                                                                                      <  1:469579/150‑1 (MQ=255)                     ggtgCCGACGCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCCCCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGGCTAc                                                                                                                                <  2:212693/150‑1 (MQ=255)                       tgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCg                                                                                                                              <  2:352823/150‑1 (MQ=255)                        gCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCg‑g                                                                                                                            >  2:624876/1‑150 (MQ=255)                               cAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCt                                                                                                                                      <  2:755083/134‑1 (MQ=255)                               cAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCt                                                                                                                                      >  1:755083/1‑134 (MQ=255)                                   aCGACCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTgg                                                                                                                 <  2:338350/150‑1 (MQ=255)                                      aTCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCgg                                                                                                              >  2:355132/1‑150 (MQ=255)                                          ccTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCCAGGTAGTGGCGGTGGa                                                                                                          >  2:362317/1‑150 (MQ=255)                                          ccTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGa                                                                                                          >  2:969876/1‑150 (MQ=255)                                                    cTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCaacaac                                                                                                >  2:909622/1‑150 (MQ=255)                                                            cgccACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGcc                                                                                        >  2:924055/1‑150 (MQ=255)                                                            cgccACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGcc                                                                                        <  2:842461/150‑1 (MQ=255)                                                             gccACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCAGGCGGTTggtgg                                                                                                                                         <  1:361701/101‑5 (MQ=255)                                                             gccACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCAGGCGGTTggtgg                                                                                                                                         >  2:361701/1‑97 (MQ=255)                                                               cACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGgt                                                                                     <  2:623146/150‑1 (MQ=255)                                                                   ggCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGCGGACATCAACAACGGTGTGCCGGTGTgg                                                                                 <  2:581904/150‑1 (MQ=255)                                                                            cgTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCcgc                                                                        <  1:355132/150‑1 (MQ=255)                                                                              ttGGATCTACGAAAACTCGCCGGCCGTGCTTACGCTGCGTGGCACCGGTGCGCCGATCGCCCCCCACCGCCTGGCGGTTGCCGGCATGTCTGACG‑GCAAGGTAGTGGCGGTTGACATCTACAACTTTTTGCCAGTGTTTGAGAGCCGCGt                                                                      <  2:852145/150‑1 (MQ=255)                                                                              ttGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGt                                                                      <  1:909622/150‑1 (MQ=255)                                                                               ctgATCTACGAAAACTCGCCGGCCGTGCTGAACAAGAGTGGCGCCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTg                                                                     <  2:819755/148‑1 (MQ=255)                                                                                 gATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGc                                                                   >  1:440958/1‑150 (MQ=255)                                                                                  aTCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGcc                                                                  >  1:209708/1‑150 (MQ=255)                                                                                  aTCTACGAAAACTCGCCGGCCGTGCTGACGCTGAGTGTCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTggag                                                                  >  1:819755/1‑148 (MQ=255)                                                                                           aaCTCGCCGGCCGTGATGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAg                                                         <  1:224416/150‑1 (MQ=255)                                                                                                gccggccgTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTCCCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGACGGGGTGGGAGAGCCGCGTGGCCATCCAGAGGGGCCg                                                    >  1:331054/1‑150 (MQ=255)                                                                                                   ggccgTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGCAGGGCCGTTc                                                 >  2:390376/1‑150 (MQ=255)                                                                                                             ccGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGAc                                       <  1:296200/149‑1 (MQ=255)                                                                                                                    gTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTgg                                <  1:362317/150‑1 (MQ=255)                                                                                                                      ggCGCCGGGGCGCCGATCGCCACCAACCGCCTGGCGGCTCCAGACGCGCCTGGCC‑GCAATGTAGGCTCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGtt                              <  1:806056/150‑1 (MQ=255)                                                                                                                       gCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTTGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTg                             <  2:288736/150‑1 (MQ=255)                                                                                                                         aCCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGat                           >  2:622377/1‑150 (MQ=255)                                                                                                                         aCCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTCT‑TACTGCGGAGGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGat                           <  2:30048/150‑1 (MQ=255)                                                                                                                         aCCGGTGCGCCGATCGCCACCAACCCTGAATCTGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGat                           >  1:30048/1‑150 (MQ=255)                                                                                                                            ggTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATc                        <  1:751447/150‑1 (MQ=255)                                                                                                                                   cGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTTTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGAcggcg                 <  1:624876/150‑1 (MQ=255)                                                                                                                                   cGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGGCCGTGTGGTTGATATCGAcggcg                 >  1:254981/1‑150 (MQ=255)                                                                                                                                   cGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGAcggcg                 <  1:969876/150‑1 (MQ=255)                                                                                                                                         ccaccaACCGCCTGGCGGTTGCGGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCttgt           <  1:472760/150‑1 (MQ=255)                                                                                                                                         ccaccaACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCTGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCttgt           >  2:312324/1‑150 (MQ=255)                                                                                                                                          caccaACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCttgtt          >  1:369960/1‑150 (MQ=255)                                                                                                                                           gcgatccGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCttgttg         <  2:171476/145‑1 (MQ=255)                                                                                                                                           accaACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCttgttg         <  2:209708/150‑1 (MQ=255)                                                                                                                                               aCCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCTTGTTGCTGt     >  1:798542/1‑150 (MQ=255)                                                                                                                                                ccGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTg                                                                                                            >  2:128814/1‑46 (MQ=255)                                                                                                                                                ccGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTg                                                                                                            <  1:128814/46‑1 (MQ=255)                                                                                                                                                  gCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCTTGTTGCTGTCtg  >  1:454973/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCG‑GCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCTTGTTGCTGTCTG  >  NC_002947/990552‑990846

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