Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	BS168_reference_genome	78,236	A→T	N418I (AAT→ATT)	ftsH →	ATP‑dependent cytoplasmic membrane protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	78,236	0	A	T	96.5%	165.6 / ‑3.8	58	N418I (AAT→ATT)	ftsH	ATP‑dependent cytoplasmic membrane protease
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (25/31);  total (26/32)
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

CTTGAATGAAGCTGCGCTTGTAGCGGCTCGTCAAAACAAGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCAATATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGCTGGCGGTTATGCTGTTATGCTGCCAAGAGAAGA  >  BS168_reference_genome/78102‑78368                                                                                                                                       |                                                                                                                                     ctTGAATGAAGCTGCGCTTGTAGCGGCTCGTCAAAACAAGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTGATATCAAAGAAAGAAAGCaa                                                                                                                                      >  7:141036/1‑135 (MQ=255)  ttGAATGAAGCTGCGCTTGTAGCGGCTCGTCAAAACAAGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCAtta                                                                                                                                    >  8:46373/1‑136 (MQ=255)    gaatgaaGCTGCGCTTGTAGCGGCTCGTCAAAACAAGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCAtta                                                                                                                                    <  6:68864/134‑1 (MQ=255)      atgaaGCTGCGCTTGTAGCGGCTCGTCAAAACAAGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGt                                                                                                                                <  7:138593/136‑1 (MQ=255)            cTGCGCTTGTAGCGGCTCGTCAAAACAAGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTa                                                                                                                          >  7:173784/1‑136 (MQ=255)             tGCGCTTGTAGCGGCTCGTCAAAACAAGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTAt                                                                                                                         >  1:56126/1‑136 (MQ=255)                gcTTGTCGCGGCTCGTCAAAACAAGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCAc                                                                                                                      <  4:237225/136‑1 (MQ=255)                gcTGGTAGCGGCTAGTCAAAACAAGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCAc                                                                                                                      <  8:166499/136‑1 (MQ=255)                  ttGTAGCGGCTCGTCAAAACAAGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACg                                                                                                                     >  3:155689/1‑135 (MQ=255)                  ttGTAGCGGCTCGTCAAAACAAGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGGTATCTCCAAGAAAGAACGCATTATCGTGGCGTATCACg                                                                                                                     >  3:92267/1‑135 (MQ=255)                       gCGGCTCGTCAAAACAAGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCg                                                                                                               >  6:110576/1‑136 (MQ=255)                         ggCTCGTCAAAACAAGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGa                                                                                                             <  5:313044/136‑1 (MQ=255)                                     cAAGAAAAAAATCGATGCGCGTGATATTTACCAAGCGACGGTCCGTGTAATTGCCTGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATc                                                                                                 <  6:272356/136‑1 (MQ=255)                                      aaGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCg                                                                                                <  8:173784/136‑1 (MQ=255)                                      aaGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCg                                                                                                <  2:322497/136‑1 (MQ=255)                                       aGAAAAAAATCGATGAGCGTGATATTGACGAAGCGACGGACAGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCgg                                                                                               <  8:300774/136‑1 (MQ=255)                                           aaaaaTCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTTTCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGtctc                                                                                           >  6:67227/1‑136 (MQ=255)                                                tCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGGTATCTCCAAGAAAGAACGCATTATCGGGGCTTATCACGAAGGCGGGCACCCCGGTTTCGGTGTCGtgtt                                                                                      >  8:329444/1‑136 (MQ=255)                                                  gATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTa                                                                                     <  3:238872/135‑1 (MQ=255)                                                  gATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTa                                                                                     <  5:110576/135‑1 (MQ=255)                                                    tGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCAATATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGAt                                                                                  <  7:256778/136‑1 (MQ=255)                                                      cgcgTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTTTCTCCAAGAAAGAACGCATTTTCGTGGGTTATCACGAAGGCGGACACCCCGGTATCGGTCGCGTGTTTGAt                                                                                  >  5:121169/1‑134 (MQ=255)                                                      cgcgTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGAt                                                                                  >  7:133372/1‑134 (MQ=255)                                                             tatTGACGCAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTATAGATGAGGCAGat                                                                         <  5:415072/136‑1 (MQ=255)                                                             tatTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGat                                                                         <  2:329826/136‑1 (MQ=255)                                                              atTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGat                                                                         <  3:118271/135‑1 (MQ=255)                                                                  aCGAAGCGACGGCCCGTGTAATTGCCGCACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGTTGAGGCAGATATgg                                                                     <  2:63642/135‑1 (MQ=255)                                                                          aCGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAg                                                                           >  8:100019/1‑121 (MQ=255)                                                                          aCGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAg                                                                           <  7:100019/121‑1 (MQ=255)                                                                          aCGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGGTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGCAGGCGGACACACCGGTATCGGTGTCGTTTTTGATGAGGCAGATATGGTTCATAAAg                                                            >  4:152728/1‑136 (MQ=255)                                                                                cGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGa                                                      <  4:155689/136‑1 (MQ=255)                                                                                    tAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACCCCGTTATCGGTCTCGTGTTAGTAGAGGCAGATTTGGTTCATAAAGTAAAGTTTGt                                                  >  5:412054/1‑136 (MQ=255)                                                                                     aaTTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGCACGCATTATCGTGGCTTATCACGAAGGCGGACCCACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGtt                                                 >  5:202092/1‑136 (MQ=255)                                                                                     aaTTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGtt                                                 >  7:241565/1‑136 (MQ=255)                                                                                     aaTTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGGTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGtt                                                 >  5:160477/1‑136 (MQ=255)                                                                                      aTTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGtt                                                 <  8:14284/135‑1 (MQ=255)                                                                                           cGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCTTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCgg                                           <  8:133372/136‑1 (MQ=255)                                                                                                 cGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAg                                     <  8:255507/136‑1 (MQ=255)                                                                                                 cGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGGGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAg                                     <  7:341656/136‑1 (MQ=255)                                                                                                   cTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGc                                   <  2:56126/136‑1 (MQ=255)                                                                                                   cTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGc                                   <  5:67227/136‑1 (MQ=255)                                                                                                   cTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGc                                   <  7:46373/136‑1 (MQ=255)                                                                                                                cgTTTTCTCCAATAAAGAACGCATTATAGTGGGTTTTGAAGAAGGCGGACACACCCTTATCGGTCTCGTTTTAGATGATGCAGCTATGGTTCATAAAGTAACGATTGTTGCTCGTGGCCAGGGGGGCGgttatgcg                      >  4:401102/1‑135 (MQ=255)                                                                                                                cgTTAACACCAAGAATCTACGCAATATCGTGGGATATCACGAAGGCGGCCACACCGTTATCGGTGTCGGTTTAGTTGAGGCAGAGATGGTTAATTTAGTAACCATTGTTCCTCGGGGCCAGGCAGGCGgtaatgct                      >  8:30895/1‑136 (MQ=255)                                                                                                                    aTCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGCTGGCGgttatgctgtt                   <  5:160701/135‑1 (MQ=255)                                                                                                                    aTCTCCAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGATTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGCTGGCGgttatgctgtt                   <  8:141036/135‑1 (MQ=255)                                                                                                                        ccAAGAAAGAACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGGCTCGGTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCACGGGGGCAGGCTGGCGgttatgctgttatgct              >  2:84304/1‑136 (MQ=255)                                                                                                                          aagaaagaACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGCTGGCGgttatgctgtt                   >  3:374104/1‑129 (MQ=255)                                                                                                                          aagaaagaACGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGCTGGCGgttatgctgtt                   <  4:374104/129‑1 (MQ=255)                                                                                                                            gaaagaaCGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGCTGGCGGTTATGCTGTTATGCTGCCa          <  8:241565/136‑1 (MQ=255)                                                                                                                              aagaaCGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGCTGGCGGTTATGCTGTTATGCTGCCAag        >  1:199804/1‑136 (MQ=255)                                                                                                                                 aaCGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGCTGGCGGTTATGCTGTTATGCTGCCAAgaga     <  4:261163/136‑1 (MQ=255)                                                                                                                                  aCGCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGCTGGCGGTTATGCTGTTATGCTGCCAagag      >  7:199892/1‑134 (MQ=255)                                                                                                                                    gCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCgg                                                                                               <  7:134745/43‑1 (MQ=255)                                                                                                                                    gCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCgg                                                                                               >  8:134745/1‑43 (MQ=255)                                                                                                                                    gCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCg                                                                                                <  6:367349/42‑1 (MQ=255)                                                                                                                                    gCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCg                                                                                                >  5:367349/1‑42 (MQ=255)                                                                                                                                    gCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGCTGGCgg                             <  2:224471/109‑1 (MQ=255)                                                                                                                                    gCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGCTGGCgg                             >  1:224471/1‑109 (MQ=255)                                                                                                                                    gCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGCTGGCGGTTATGCTGTTATGCTGCCAAGagaaga  >  2:388522/1‑136 (MQ=255)                                                                                                                                    gCATTATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGCTGGCGGTTATGCTGTTATGCTGCCAAGagaaga  >  4:47469/1‑136 (MQ=255)                                                                                                                                       |                                                                                                                                     CTTGAATGAAGCTGCGCTTGTAGCGGCTCGTCAAAACAAGAAAAAAATCGATGCGCGTGATATTGACGAAGCGACGGACCGTGTAATTGCCGGACCCGCTAAGAAGAGCCGCGTTATCTCCAAGAAAGAACGCAATATCGTGGCTTATCACGAAGGCGGACACACCGTTATCGGTCTCGTTTTAGATGAGGCAGATATGGTTCATAAAGTAACGATTGTTCCTCGGGGCCAGGCTGGCGGTTATGCTGTTATGCTGCCAAGAGAAGA  >  BS168_reference_genome/78102‑78368

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