Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	2,831,221	A→G	12.9%	L187S (TTG→TCG)	yrbG ←	conserved hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	2,831,221	0	A	G	12.9%	72.1 / 5.4	31	L187S (TTG→TCG)	yrbG	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/16);  new base G (2/2);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

TTCCCATTATCGGACGGTTGATCATCCGATAATGGTCACGTGCGGTGCCCATCTTTTTCATCAATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGTCTATAATGAGCGGCATTTCGAGCTGCCGATGCTCTCCGCTGTTTTCTTTCTTAACAATA  >  BS168_reference_genome/2831089‑2831346                                                                                                                                     |                                                                                                                              ttCCCATTATCGGACGGTTGATCATCCGATAATGGTCACGTGCGGTGCCCATCTTTTTCATCAATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACttt                                                                                                                         >  5:137235/1‑139 (MQ=255)    ccATTATCGGACGGTTGATCATCCGATAATGGTCACGTGCGGTGCCCATCTTTTTCATCAATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCa                                                                                                                       <  6:206088/138‑1 (MQ=255)              aCGGTTGATCATCCGATAATGTTCACGTGCGGTGCACATCTTTTTCATCAATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCa                                                                                                                               >  4:254787/1‑120 (MQ=255)              aCGGTTGATCATCCGATAATGGTCACGTGCGGTGCCCATCTTTTTCATCAATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCa                                                                                                                               <  3:254787/120‑1 (MQ=255)                ggTTGATCATCCGATAATGGTCACGTGCGGTGCCCATCTTTTTCATCAATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTc                                                                                                          <  6:80813/139‑1 (MQ=255)                          ccGATAATGGTCACGTGCGGTGCCCATCTTTTTCATCAATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAagcag                                                                                                                  >  8:55395/1‑121 (MQ=255)                          ccGATAATGGTCACGTGCGGTGCCCATCTTTTTCATCAATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAagcag                                                                                                                  <  7:55395/121‑1 (MQ=255)                            gATAATGGTCACGTGCGGTGCCCATCTTTTTCATCAATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGa                                                                                              >  3:104680/1‑139 (MQ=255)                                   gTCTCGTGCGGTGCCCATCTTTTTCATCTATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGa                                                                                       <  2:148237/139‑1 (MQ=255)                                      aCGTGCGGTGCCCATCTTTTTCATCAATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAg                                                                                    >  2:126893/1‑139 (MQ=255)                                           cGGTGCCCATCTTTTTCATCAATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGtttt                                                                               <  4:52159/139‑1 (MQ=255)                                            ggTGCCCATCTTTTTCATCAATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTc                                                                              >  2:221756/1‑139 (MQ=255)                                                         ttCATCAATATATATTTATCCTTCTGTACACGAGCAAAATGAAATATCTGACGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTtgaatgaa                                                                  <  8:29971/138‑1 (MQ=255)                                                         ttCATCAATATATATTTATACGTCTGTACACGAGCAAAATGAAATATCTGAAGGGGTTGTGTATCCGTGTTTTTGCAAACTTTAAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTtgaatgaa                                                                  <  8:25627/138‑1 (MQ=255)                                                          tcatcaATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTTTCTTTACTGATTCTTGACAGGTTTTCTGTTTGAATGAAg                                                                 <  5:149803/138‑1 (MQ=255)                                                             tcaATATATATTTCTCCGTCTGTAAACGAGCAAAATGAACTATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGt                                                             <  8:27748/139‑1 (MQ=255)                                                             tcaATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGt                                                             <  1:221756/139‑1 (MQ=255)                                                               aaTATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCGAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGTc                                                            >  1:264030/1‑138 (MQ=255)                                                                atatatatTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCGAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGc                                                                >  6:228080/1‑133 (MQ=255)                                                                atatatatTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCGAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGc                                                                <  5:228080/133‑1 (MQ=255)                                                                 tatatatTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGTCta                                                          >  4:230057/1‑138 (MQ=255)                                                                     tatTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGTCTATAATg                                                     <  5:223735/139‑1 (MQ=255)                                                                          ctcCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGTCTATAATGAGCg                                                 <  3:128223/138‑1 (MQ=255)                                                                           tcCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGTCTATAATGAGCgg                                                >  6:133817/1‑138 (MQ=255)                                                                           tcCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGTCTATAATGAGCgg                                                >  1:21724/1‑138 (MQ=255)                                                                             cGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGtttt                                                                               >  7:187479/1‑105 (MQ=255)                                                                             cGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGTCTATAATGAGCGGCAt                                             >  5:48187/1‑139 (MQ=255)                                                                             cGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGACATTTTCTATCTTTACTGATTCTGGACATGtttt                                                                               <  8:187479/105‑1 (MQ=255)                                                                                     aaCGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGTCTATAATGAGCGGCATTTCGAGc                                      <  6:137235/138‑1 (MQ=255)                                                                                                                tttGTGTATCCGTGTTTTTGCGAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGTCTATAATGAGCGGCATTTCGAGCTGCCGATGCTCTCCGCTGTTttctttct          <  2:264030/139‑1 (MQ=255)                                                                                                                tttGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGTCTATAATGAGCGGCATTTCGAGCTGCCGATGCTCTCCGCTGTTttctttct          <  3:284727/139‑1 (MQ=255)                                                                                                                        tCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGTCTATAATGAGCGGCATTTCGAGCTGCCGATGCTCTCCGCTGTTTTCTTTCTTAACaa    >  5:203276/1‑137 (MQ=255)                                                                                                                        tCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGTCTATAATGAGCGGCATTTCGAGCTGCCGATGCTCTCCGCTGTTTTCTTTCTTAACAATa  <  4:104680/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                              TTCCCATTATCGGACGGTTGATCATCCGATAATGGTCACGTGCGGTGCCCATCTTTTTCATCAATATATATTTCTCCGTCTGTAAACGAGCAAAATGAAATATCTGAAGGGTTTGTGTATCCGTGTTTTTGCAAACTTTCAAGCAGCCATTTTCTATCTTTACTGATTCTGGACAGGTTTTCTGTTTGAATGAAGCCGTCTATAATGAGCGGCATTTCGAGCTGCCGATGCTCTCCGCTGTTTTCTTTCTTAACAATA  >  BS168_reference_genome/2831089‑2831346

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