Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	1,511,878	C→A	24.8%	Q191K (CAA→AAA)	sipT →	type I signal peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	1,511,878	0	C	A	24.8%	93.7 / 8.9	53	Q191K (CAA→AAA)	sipT	type I signal peptidase
Reads supporting (aligned to +/- strand):  ref base C (28/11);  new base A (2/12);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 3.20e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAAGGCAAATACTTTGTCATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAAAGAAGCGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCACACTGAAGAAAAGGTGAACCTG  >  BS168_reference_genome/1511744‑1512015                                                                                                                                       |                                                                                                                                          cGAAAGGCAAATACTTTGTCATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTGCCGTTTAACGAAATGCGTcaaa                                                                                                                                        >  1:277674/1‑138 (MQ=255) cGAAAGGCAAATACTTTGTCATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTcaaa                                                                                                                                        >  1:134810/1‑138 (MQ=255) cGAAAGGCAAATACTTTGTCATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTcaaa                                                                                                                                        >  4:183964/1‑138 (MQ=255) cGAAAGGCAAATACTTTGTCATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTcaaa                                                                                                                                        >  5:10869/1‑138 (MQ=255) cGAAAGGCAAATACTTTGTCATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTcaaa                                                                                                                                        <  7:189359/138‑1 (MQ=255) cGAAAGGCAAATACTTTGTCATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTcaaa                                                                                                                                        >  7:94259/1‑138 (MQ=255) cGAAAGGCAAATACTTTGACATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGTAATGCGTcaaa                                                                                                                                        >  8:116775/1‑138 (MQ=255)           aTACTTTGTCATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATaaaaa                                                                                                                             >  8:211731/1‑139 (MQ=255)             gacaggTCATGGTAGATAACCGGCTTAATTCTATGGACAGCCGAAACGGGCTGGTACTTATCGCGTAATATCGAATTGTCGCCACATCGAAGTTTGTCTTTTTCCCTTTTAACTAAATGCGTCAAAAAAAATAAAAACg                                                                                                                           <  6:29516/134‑1 (MQ=255)                 tGTCATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTGCCGTTTAACGAAATGCGTAAAACAAAATAAAAACGCCtt                                                                                                                       >  7:88504/1‑139 (MQ=255)                 tGTCATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCtt                                                                                                                       >  7:237137/1‑139 (MQ=255)                 tGTCATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCtt                                                                                                                       >  7:156535/1‑139 (MQ=255)                 tGTCATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCtt                                                                                                                       >  3:118550/1‑139 (MQ=255)                 tGTCATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAAAAAAAACGccat                                                                                                                       >  5:254695/1‑137 (MQ=255)                   tCATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGc                                                                                                                     <  7:199116/139‑1 (MQ=255)                           gaTAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACCAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAg                                                                                                             <  7:201996/139‑1 (MQ=255)                             tAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGccctgtc                                                                                                           >  5:29516/1‑134 (MQ=255)                             tAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAAAGAAATGCGTCAAACAAAATAAAAACGACTTGCTGGCCTAGGa                                                                                                           >  3:237846/1‑139 (MQ=255)                                   gCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGCAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGAcagcag                                                                                                     >  4:32555/1‑139 (MQ=255)                                   gCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGCcagcag                                                                                                     >  8:222890/1‑139 (MQ=255)                                   gCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGAcagcag                                                                                                     >  8:69260/1‑139 (MQ=255)                                           cTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGttttt                                                                                             >  2:171560/1‑139 (MQ=255)                                             aTGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTAt                                                                                           >  4:122506/1‑139 (MQ=255)                                              tGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTGtttttt                                                                                          >  3:103611/1‑139 (MQ=255)                                              tGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTGtttttt                                                                                          >  3:156491/1‑139 (MQ=255)                                                   aGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTAAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTATTTTTgt                                                                                      >  7:110287/1‑137 (MQ=255)                                                     ccGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGCCAGCAGGCGTGttttttttttaa                                                                                    >  5:16631/1‑134 (MQ=255)                                                        aaaCGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTTTTTTTtaaaat                                                                                 >  6:5817/1‑137 (MQ=255)                                                         aaCGGGCGGGTACTGATCGCGCAAGATGTAATTGGCGGCACATCGATGTTTGTGTTTTACCCGTTTAACTAAAGGCGAAAAACAAAATAAAACCGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAaagaa                                                                               <  7:290914/139‑1 (MQ=255)                                                             ggCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATCCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAAAGAAGCGa                                                                           >  8:203762/1‑139 (MQ=255)                                                             ggCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCATTTTACCGTTTAACGAAATGCGTCAAACAAAAAAAAAACGCCATGCTGGCCTAGGACAGCAGGCGTTTTTTTTTTTATaaaaaaaacg                                                                           >  6:42584/1‑135 (MQ=255)                                                              gCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAAAAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAAAGAAGCGAc                                                                          <  1:24902/139‑1 (MQ=255)                                                                     cttacggcgTAATATCTAATTTCCGGCACATCTATGTTTTTCTTTTTCCCGTTAAACGAAATGCGACAAAAAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAAAGAAGCGACcaatgca                                                                   <  3:212022/133‑1 (MQ=255)                                                                     cttaTCGGGCAAGACCTAATTGCCGGCAATACCAAGTTTGTCTTTTACCCGTTAAACTAAATGCGTCAAAAAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAAAGAAGCGACcaatgca                                                                   <  5:5817/136‑1 (MQ=255)                                                                     cTGATCGCGGAATAGAGAATTGGCGGCACATCGATGTTTGTCTTTTTCCCGTTTAACGAAATGCGTAAAACAAAAAAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAAAGAAGCGACcaatgca                                                                   <  3:189879/139‑1 (MQ=255)                                                                     cTGATCCCGTAATATGTAATTGTCGCCAAATCGAAGTTTTAATTTTACCCTTTAAACCAAATCCGTAAAAAAAAAAAAAAACGCCTTGCTGGCCTACGACAGCAGGCTTTTTTATTTTGAAAAAGAAGCGACcaatgca                                                                   <  8:88504/139‑1 (MQ=255)                                                                        aTCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTTTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCGTAGGACAGCAGGCGTTTTTATTTTGTAAAATAAGCGCCcaattcaat                                                                 >  5:140839/1‑138 (MQ=255)                                                                        aTCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGc                                                                                                                     <  1:69748/86‑1 (MQ=255)                                                                        aTCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGc                                                                                                                     >  2:69748/1‑86 (MQ=255)                                                                        aTCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTTTAAAATAAACGACcaattcaaag                                                                >  4:212022/1‑137 (MQ=255)                                                                        aTCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTGTTTATTTTTAAAAATAATCGCCcaaaggaaa                                                                 >  1:208332/1‑132 (MQ=255)                                                                        aTCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCATGCTGGCCTAGGACAGCAGGCGTGTTTTTTTTTAAAAAGAAAGGCCcaaagcaaaa                                                                >  4:204763/1‑137 (MQ=255)                                                                        aTCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAAAGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTGTTTATTTTGGAAAATAAGCGACcaatgcaagg                                                                >  5:49299/1‑137 (MQ=255)                                                                        aTCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGGTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGGTTTTATTTTGTAAAAGCAACGACcaaagcaatg                                                                >  2:163652/1‑139 (MQ=255)                                                                           ggggAAGTTCTATTTGTCGGAACATCGAAGTTTGTCTTTTTCCCTTTTAACTAAATGCTTCAAAAAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAAAGAAGCGACCAatgcaatgca                                                              <  2:208332/136‑1 (MQ=255)                                                                                 gATCGAATTGGCGGCACATCGATGTTTGTCTTTTTCCCTTTTAAGTAAATGCTTCAAACAAAAAAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAAAGAAGCGACCAATGCAATGCATAAAAcc                                                       <  4:118550/139‑1 (MQ=255)                                                                                  aTCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTGTTTTTTTTGAAAAAGAAGCGCCCAAAGCAATGCATAAAAccc                                                      >  3:287269/1‑138 (MQ=255)                                                                                           tCGGCAAATCGAAGTTTGTTTTTTTCCCTTTAAACTAAATGCGAAAAAAAAAAAAAAAACGCCGGGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAAAGAAGCGACCAATGCAATGCATAAAACCGTGATGAaca                                             <  3:122506/139‑1 (MQ=255)                                                                                           gcGGCACATCTAAGTTTGTATTTTACCATTTAAAACAAAGGCAAAAAAAAAAAAAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAAAGAAGCGACCTATGCAATGCATAAAACCGTGATGAaca                                             <  3:164286/138‑1 (MQ=255)                                                                                                 caTCGATTTTTGTATTTTTCCCGTTAAACTAAATGCCTAAAACAAAAAAAAAACGCCTTGCTGGCGTAGGACAGCAGGCGTTTTTATTTTTAAAAAGAAGCGACCAATGCAATGCATAAAACCGTGATGAACACTGAAc                                       <  7:211731/139‑1 (MQ=255)                                                                                                 caTCGAAGTTTGTCTTTTTCCCGTTTACCGAAATGCGTCAAACAAAATAACACCGCCTTGGTGGCCTAGGCCAGCAGGCGTTTTTATTTTGAAAAAGAAGCGACCAATGCAATGCATAAAACCGTGATGAACACTGAAc                                       <  8:191372/139‑1 (MQ=255)                                                                                                      acgTTTGTCTTTTCCCCTTTTAACTAAATGCGTAAAACAAAAAAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAAAGAAGCGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGata                                  <  7:222890/137‑1 (MQ=255)                                                                                                      aaGTTTGTCTTTTTCCCGTTAAACGAAATGCGTCAAAAAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAAAGAAGCGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGata                                  <  8:94259/139‑1 (MQ=255)                                                                                                         tttgagtttttCCCGCTTAACTAAAGTCCAAAAACAAAATAAAAACGCCTTGCGTGCATATGACAGAAGGCGTTTTTATTTTGAAAAAGAAGCGACCAATTCAATGCATAAAACCGTGATGAACACTGAACCGataat                                <  8:153640/132‑1 (MQ=255)                                                                                                         tttGTCTTTTCCCTTTTAAACTAAATGCCAAAAAAAAAAAAAAAACGCCTTGGTGGCCTAGTACAGCTGGGTTTTTTATTTTGAAAAAGAAGCGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGataat                                <  6:12161/138‑1 (MQ=255)                                                                                                         tttGTCTTTTACCCTTAAAACAAATGGCCAAAAAAAAAAAAAACACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAAAGAAGCGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGataat                                <  1:171560/138‑1 (MQ=255)                                                                                                                  ttCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTGTTTTTTTTTTAAAAGAACCGACCAAAGCAAAGAAAAAAAACGGGTTGTAAACTGACACCCTAATAGACAGcacc                      >  5:12161/1‑138 (MQ=255)                                                                                                                          tAACTAAATGCGTAAAAAAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTTAAAAAGAAGCGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCACACTgaaga               <  7:10457/138‑1 (MQ=255)                                                                                                                                    cGTAAAAAAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGGGTTTTTATTTTGAAAAAGAAGCGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCACACTGAAGAAAAGGTGAAcc    <  3:101506/139‑1 (MQ=255)                                                                                                                                      tCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAAAGAAGCGCCCAATGCAATGCATAAAACCCTGATGAACACAGAACCGATAATAGCCAGCACCCAGAAGAAAAGGTGAACCtg  >  4:101506/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                          CGAAAGGCAAATACTTTGTCATGGGAGATAACCGGCTGAATTCTATGGACAGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTTTTTATTTTGAAAAAGAAGCGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCACACTGAAGAAAAGGTGAACCTG  >  BS168_reference_genome/1511744‑1512015

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