Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,416,739	A→T	7.8%	D370V (GAC→GTC)	PP_1238 →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,416,739	0	A	T	7.8%	145.9 / 3.1	51	D370V (GAC→GTC)	PP_1238	lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (24/23);  new base T (2/2);  total (26/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 = 5.20e-01

AGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAGGTGGGTGAGAGCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCT  >  NC_002947/1416602‑1416887                                                                                                                                          |                                                                                                                                                     agagCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAGGTGGGTGAGAGCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCt                                                                                                                                           <  7:219107/149‑1 (MQ=255)                     aGCGTTACCAGGTTCGCCTGAGCAAGGTGGGTGAGAGCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTc                                                                                                                       <  1:111826/149‑1 (MQ=255)                       cGTTACCAGGTTCGCCTGAGCAAGGTGGGTGAGAGCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCgg                                                                                                                     <  5:160075/149‑1 (MQ=255)                            ccAGGTTCGCCTGAGCAAGGTGGGTGAGAGCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTcgcg                                                                                                                             <  3:97546/136‑1 (MQ=255)                            ccAGGTTCGCCTGAGCAAGGTGGGTGAGAGCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTcgcg                                                                                                                             >  4:97546/1‑136 (MQ=255)                             cAGGTTCGCCTGAGCAAGGTGGGTGAGAGCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCgg                                                                                                               >  3:76818/1‑149 (MQ=255)                              aGGTTCGCCTGAGCAAGGTGGGTGAGAGCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCgg                                                                                                               >  1:140668/1‑148 (MQ=255)                               ggTTCGCCTGAGCAAGGTGGGTGAGAGCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCa                                                                                                             <  3:203339/149‑1 (MQ=255)                                    gCCTGAGCAAGGTGGGTGAGAGCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGcc                                                                                                          >  5:49033/1‑147 (MQ=255)                                                   gtgAGAGCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCt                                                                                          >  4:99190/1‑148 (MQ=255)                                                      agagagTGCAGGTCACGGTCGAGAAAAACATCAACGCCTGGGCGCCGGCCGATGTCGGCCGCGGCGTGCTGAGATCTTTTCAGGACCACCTTGGCTTAGTGCGCTTCTCGGTACTCGCTATCGGCAGCCAAGGGAACGGCACGCTGAt                                                                                       <  6:260020/148‑1 (MQ=255)                                                      agagCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGTCCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATc                                                                                      >  4:124980/1‑149 (MQ=255)                                                      agagCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATc                                                                                      >  1:222952/1‑149 (MQ=255)                                                       gagCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAAc                                                                                                  >  3:54675/1‑136 (MQ=255)                                                       gagCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAAc                                                                                                  <  4:54675/136‑1 (MQ=255)                                                             gCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAgg                                                                                                      <  1:11139/126‑1 (MQ=255)                                                             gCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAgg                                                                                                      >  2:11139/1‑126 (MQ=255)                                                               aGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCCGCCAAGGGAACGGCACGCTGATCGCCAgtggt                                                                             >  6:209296/1‑149 (MQ=255)                                                               aGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGCCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGCGAACGGCACGCTGATCGCCAgtggt                                                                             >  2:326373/1‑149 (MQ=255)                                                                ggTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAgtggtg                                                                            <  2:222952/149‑1 (MQ=255)                                                                    aCGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGacac                                                                        <  6:49033/149‑1 (MQ=255)                                                                       gTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGtt                                                                     <  4:169702/149‑1 (MQ=255)                                                                          gagaAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCAt                                                                  >  6:94456/1‑149 (MQ=255)                                                                                aaCATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGTCCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGt                                                            <  3:124980/149‑1 (MQ=255)                                                                                      aaCACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTg                                                      <  7:88686/149‑1 (MQ=255)                                                                                       acacCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCg                                                                                     >  2:224739/1‑117 (MQ=255)                                                                                       acacCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCg                                                                                     <  1:224739/117‑1 (MQ=255)                                                                                            gTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCtt                                                <  2:140668/149‑1 (MQ=255)                                                                                                    ggCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCaa                                                                                                        <  6:60466/85‑1 (MQ=255)                                                                                                    ggCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGc                                        <  3:6580/149‑1 (MQ=255)                                                                                                    ggCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGc                                        <  5:94456/149‑1 (MQ=255)                                                                                                     gCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCg                                       >  8:142672/1‑149 (MQ=255)                                                                                                       cGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGt                                      >  8:19829/1‑148 (MQ=255)                                                                                                            tCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAAcc                                >  3:16831/1‑149 (MQ=255)                                                                                                              gcccgccGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAAcc                                >  4:167197/1‑147 (MQ=255)                                                                                                               cccgccGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAg                             >  1:201621/1‑149 (MQ=255)                                                                                                                     gcgTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGgcgc                       <  7:142672/149‑1 (MQ=255)                                                                                                                         gCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTgg                   >  6:86466/1‑149 (MQ=255)                                                                                                                             aGCGCTATTCAGGTCCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCgg                                                   >  1:42340/1‑113 (MQ=255)                                                                                                                             aGCGCTATTCAGGTCCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCgg                                                   <  2:42340/113‑1 (MQ=255)                                                                                                                               cgcTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGACAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCgctgc             >  1:81720/1‑149 (MQ=255)                                                                                                                                 cTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTc           <  3:167197/149‑1 (MQ=255)                                                                                                                                  tATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCg          >  1:35958/1‑149 (MQ=255)                                                                                                                                  tATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCg          >  4:225902/1‑149 (MQ=255)                                                                                                                                    ttCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCgg         <  3:225902/148‑1 (MQ=255)                                                                                                                                      cAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGt      >  5:287726/1‑149 (MQ=255)                                                                                                                                      cAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCAGCGCCAGGCGCCGCTCggtt       <  2:201621/148‑3 (MQ=255)                                                                                                                                       aGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGtt     >  4:192898/1‑149 (MQ=255)                                                                                                                                        ggACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGCACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGCCGCTGCTCGGGGttt    >  3:304883/1‑149 (MQ=255)                                                                                                                                        ggACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGttt    <  7:19829/149‑1 (MQ=255)                                                                                                                                        ggACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGttt    <  4:76818/149‑1 (MQ=255)                                                                                                                                          accaccTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCt  <  2:51680/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                     AGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAGGTGGGTGAGAGCGTGCAGGTCACGGTCGAGAAAAACATCAACACCGTGGCTCCGGCCGATGTCGCCCGCCGCGTGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCT  >  NC_002947/1416602‑1416887

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