Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,867,779	T→C	8.0%	K262R (AAG→AGG)	lgt ←	prolipoprotein diacylglyceryl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,867,779	0	T	C	8.0%	137.4 / 2.4	50	K262R (AAG→AGG)	lgt	prolipoprotein diacylglyceryl transferase
Reads supporting (aligned to +/- strand):  ref base T (27/19);  new base C (2/2);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

GCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTTGTGCATCCGGCACCCGGACAAACTCCACGGCGAAGCGGAAGATGC  >  NC_002947/5867631‑5867926                                                                                                                                                      |                                                                                                                                                    gCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCt                                                                                                                                                     >  7:159882/1‑149 (MQ=255)       tccGGTTCGGTTTCCCTTCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGTCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGTCGGGTTTCGCCGTAGGCTTGCGa                                                                                                                                                <  7:50356/147‑1 (MQ=255)          ggTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATtat                                                                                                                                            <  6:19806/149‑1 (MQ=255)                   tCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCccac                                                                                                                                    >  7:57532/1‑148 (MQ=255)                        gCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCCTGCGATTATAAGCCCACCAGAt                                                                                                                               >  2:149019/1‑148 (MQ=255)                                 cGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAgg                                                                                                                     >  6:309163/1‑149 (MQ=255)                                   ttATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGc                                                                                                                    >  2:258987/1‑148 (MQ=255)                                        aaTGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCCTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGa                                                                                                              >  8:213310/1‑149 (MQ=255)                                              gTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCTGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACt                                                                                                            <  2:247386/145‑1 (MQ=255)                                              gTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACt                                                                                                            >  1:247386/1‑145 (MQ=255)                                              gTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatca                                                                                                        >  1:298369/1‑149 (MQ=255)                                              gTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATAAATCAGGCGGGTTTCGCCGTAGGCTTGAGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatca                                                                                                        >  2:13319/1‑149 (MQ=255)                                                cgcgGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTatcatc                                                                                                      >  6:218617/1‑149 (MQ=255)                                                 gcgGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCg                                                                                                     <  8:57532/149‑1 (MQ=255)                                                    gACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCCTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGaa                                                                                                  <  7:213310/149‑1 (MQ=255)                                                       cAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCg                                                                                                     >  4:267408/1‑143 (MQ=255)                                                       cAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCg                                                                                                     <  3:267408/143‑1 (MQ=255)                                                       cAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGc                                                                                               <  4:86540/149‑1 (MQ=255)                                                           tCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGa                                                                                                                       <  6:58070/121‑1 (MQ=255)                                                           tCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGa                                                                                                                       >  5:58070/1‑121 (MQ=255)                                                               gATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAAc                                                                                       <  3:78419/149‑1 (MQ=255)                                                               gATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAAc                                                                                       >  5:16585/1‑149 (MQ=255)                                                                  aCTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATATCCCAACCAGAACAAAACGAGGCAACCGACTATCATCGGAATGAAAAGCAAATgg                                                                                    >  5:220217/1‑149 (MQ=255)                                                                     gTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCt                                                                                 >  2:243752/1‑149 (MQ=255)                                                                      tttCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATTAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCtt                                                                                >  4:157473/1‑149 (MQ=255)                                                                      tttCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCtt                                                                                >  6:130345/1‑149 (MQ=255)                                                                         aaatggttACCTGTATCAAAAAAGGCATGGG‑GCTCGGCCCTTGCGTGGGTTTTATCATGCGGGTTTCGGCGTTTGCTTGCGTTTAAAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtga                                                                              <  6:220217/142‑1 (MQ=255)                                                                             gCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTtgagtgag                                                                         <  8:159882/149‑1 (MQ=255)                                                                                     gTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCg                                                                  >  5:47467/1‑148 (MQ=255)                                                                                      tAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGaa                                                                >  4:125943/1‑149 (MQ=255)                                                                                        gCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGTCAACCGACTATCATCCGAATGCACAGCAACTGTCCTTGAGTGAGCCAGCCGAAcg                                                              >  6:36591/1‑149 (MQ=255)                                                                                          aaaaaCGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAAcgc                                                             >  1:168944/1‑148 (MQ=255)                                                                                           aaaaCGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGa                                                           <  1:243752/149‑1 (MQ=255)                                                                                            aaaCGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGAt                                                          <  1:122673/149‑1 (MQ=255)                                                                                                ggCAGGGGCGATC‑CCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGa                                                           >  1:15386/1‑144 (MQ=255)                                                                                                ggCAGGGGCGATC‑CCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGa                                                           <  2:15386/144‑1 (MQ=255)                                                                                                ggCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTaa                                                      <  2:168944/149‑1 (MQ=255)                                                                                                  cAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAAcc                                                    <  5:309163/149‑1 (MQ=255)                                                                                                  cAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAAcc                                                    <  6:16585/149‑1 (MQ=255)                                                                                                    gggg‑gATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAAcc                                                    >  4:159612/1‑147 (MQ=255)                                                                                                          aTCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTtgt                                             <  3:125943/149‑1 (MQ=255)                                                                                                              ccccTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCg                                                                                                     >  5:90168/1‑89 (MQ=255)                                                                                                              ccccTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCg                                                                                                     <  6:90168/89‑1 (MQ=255)                                                                                                               cccTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTTGTGCATc                                        >  6:79930/1‑149 (MQ=255)                                                                                                               cccTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCCTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTTGTGCATc                                        <  1:149019/149‑1 (MQ=255)                                                                                                                ccTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTTGTGCATcc                                       >  1:282917/1‑149 (MQ=255)                                                                                                                      gTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTTGTGCATCCGGCAcc                                 >  6:140941/1‑149 (MQ=255)                                                                                                                                        gTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTTGTGCATCCGGCACCCGGACAAACTCCACGGc                >  7:159717/1‑148 (MQ=255)                                                                                                                                            cGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTTGTGCATCCGGCACCCGGACAAACTCCACGGCGAAGc           >  2:177211/1‑149 (MQ=255)                                                                                                                                                gTAGGCTTGCGATTATAAGCACACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTTGTGCATCCGGCACCCGGACAAACTCCACGGCGAAGCgg         >  1:138573/1‑147 (MQ=255)                                                                                                                                                    gCTTGCGATTATAAGCCCACCAGATTAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTTGTGCATCCGGCACCCGGACAAACTCCACGGCGAAGCGGAAGATg   <  3:157473/149‑1 (MQ=255)                                                                                                                                                      ttGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTTGTGCATCCGGCACCCGGACAAACTCCACGGCGAAGCGGAAGATGc  >  7:206341/1‑148 (MQ=255)                                                                                                                                                      ttGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTTGTGCATCCGGCACCCGGACAAACTCCACGGCGAAGCGGAAGATGc  >  3:180665/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                    GCGGATGCCGGTGCGGTTTCCCTGCAACGTGCCGTTATCAATGACGTCGCGGACCAGGTCCAGATACTGTTTCATGGCTTACCTGTAGCAAAAACGGCAGGGG‑GATCGCCCCTGCCGTGGGATTAATCAGGCGGGTTTCGCCGTAGGCTTGCGATTATAAGCCCACCAGATCAAGACGAGGCCACCGACTATCATCGGAATGCACAGCAACTGGCCTTGAGTGAGCCAGCCGAACGCGATGTAACCGAGTTGTGCATCCGGCACCCGGACAAACTCCACGGCGAAGCGGAAGATGC  >  NC_002947/5867631‑5867926

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