Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,766,421	T→G	20.3%	S53S (TCA→TCC)	ywnA ←	putative transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,766,421	0	T	G	20.3%	77.1 / 4.2	46	S53S (TCA→TCC)	ywnA	putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (8/28);  new base G (8/2);  total (16/30)
Fisher's exact test for biased strand distribution p-value = 1.46e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTACCGGACATTTCGGATTCGGGTTTTCATGAACCGCAAAGAGTTCTTCTTGCTTTTGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGATTTCTGAAGATGTTTTTTCATCCATCGAGATGAGGGATAAAATATGAATGGCAAC  >  BS168_reference_genome/3766288‑3766558                                                                                                                                      |                                                                                                                                          cTACCGGACATTTCGGATTCGGGTTTTCATGAACCGCAAAGAGTTCTTCTTGCTTTTGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGAtgt                                                                                                                                      <  1:122349/139‑1 (MQ=255)                      ggTTTTCATGAACCGCAAAGAGTTCTTCTTGCTTTTGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTcag                                                                                                                 <  4:17549/139‑1 (MQ=255)                      ggTGTTCATGAACCGCAAAGAGTTCTCCTTGCTTTTGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTcag                                                                                                                 <  6:261987/139‑1 (MQ=255)                                        agagTTCTTCTTGCTTTTGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCg                                                                                               <  3:5937/139‑1 (MQ=255)                                          agTTCTTCTTGCTTTTGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGca                                                                                             <  6:129168/139‑1 (MQ=255)                                            ttcttcTTGCTTTTGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGGGGTGTGTGAATTTAAGCGTTTTTTAGCAGGCGGCTCATTGTTGGcaca                                                                                           >  4:237076/1‑139 (MQ=255)                                               ttcttGCTTTTGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATTTGACAAAATAATCGTTTTTAACCAGGCTGAAAATTATTCGCACAAc                                                                                         >  6:392018/1‑138 (MQ=255)                                                    gCTTTTGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGGGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTATTCGCACAACTACCgg                                                                                   >  3:397775/1‑139 (MQ=255)                                                       tttGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCg                                                                                               <  7:32750/124‑1 (MQ=255)                                                       tttGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCg                                                                                               >  8:32750/1‑124 (MQ=255)                                                        ttGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGAtt                                                                                <  7:244266/138‑1 (MQ=255)                                                        ttGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGAtt                                                                                <  5:48145/138‑1 (MQ=255)                                                        ttGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGAtt                                                                                <  8:197207/138‑1 (MQ=255)                                                        ttGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGAtt                                                                                <  4:153881/138‑1 (MQ=255)                                                                   ggTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACt                                                                    <  8:283350/139‑1 (MQ=255)                                                                     tACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGt                                                                  <  2:66771/139‑1 (MQ=255)                                                                          ttCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGGAGTGAGAAAATAAGCATTTTTCAAGAGGGTGTTCATTGTTGGCCCAACTACCGGCTTCGTGTTTTCCCTGTCTGc                                                              >  3:161683/1‑138 (MQ=255)                                                                          ttCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGGGATGTGAGAATATCAGCTTTTTTcagcag                                                                                                              <  4:4794/90‑1 (MQ=255)                                                                          ttCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGGGATGTGAGAATATAAGCTTTTTTcagcag                                                                                                              >  3:4794/1‑90 (MQ=255)                                                                          ttCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGGGATGTGAGAATATAAGCTTTTTTCAGCAGGCTGATCATTCTTCGCCCCACTACCGGGTTCGTGTTGACACTGGCTGc                                                              >  3:385804/1‑138 (MQ=255)                                                                           tCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGGTGTGGGAATATAAGCTTTTTTGCGCCGGGGGGTCATTGTTCGCCCCACAACCGGCTTCGTGTTGGCCCTggagtga                                                             >  2:352496/1‑132 (MQ=255)                                                                            cTAAACGAGAAATATCAGGAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCa                                                             <  4:285888/137‑1 (MQ=255)                                                                            cTAAACGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCa                                                             <  4:291591/137‑1 (MQ=255)                                                                            cTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCa                                                             <  4:131064/137‑1 (MQ=255)                                                                            cTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGGGATGTGAGAATATCAGCTTTTTTAAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACCCTGTCTGCAAt                                                           >  3:131064/1‑139 (MQ=255)                                                                               aaaGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGGGATGTGAGAATATCAGCTTTTTTAAGCAGGCGGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGAt                                                        >  7:60711/1‑139 (MQ=255)                                                                                   agaAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACt                                                                    <  7:395666/123‑1 (MQ=255)                                                                                   agaAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACt                                                                    >  8:395666/1‑123 (MQ=255)                                                                                    gaAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCGGATCATTCTTCGCACAACTACCCGATTCGTGTTGACACGGTCTGCAATTATTTCt                                                    >  7:60327/1‑138 (MQ=255)                                                                                    gaAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAAAATAAGCTTTTTTAATCAGGCTGATCATTATTAGCACAACAACCCGATTAGTGTTGACACAGTCTGCAATGAtttat                                                    >  6:92566/1‑136 (MQ=255)                                                                                              cagGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGGGTTGTGATAATATAAGCTTTTTTAAGCAGGCGGATCATTGTTCGCCCCACAAACCGCTTCGTGTTGACACTGTGTGCAAATATTTATGAATATGtttt                                         >  6:364049/1‑139 (MQ=255)                                                                                                 gATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGATTTCTGAAGATGTTTTTTc                                      <  5:392018/139‑1 (MQ=255)                                                                                                     ttctttGAGGCCTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGATTTCTGAAGATGTTTTTTcat                                    <  2:327687/134‑1 (MQ=255)                                                                                                             ggCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGATTTCTGAAGATGTTTTTTCATCCATCGAGAt                          <  6:307763/139‑1 (MQ=255)                                                                                                              gCTTGCACCCGGCACCCCGGCGCGTTATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGTTTCGTGTTGACACTGTCTGCAATGATTTCTGAAGATGTTTTTTCATCCATCGAGAt                          <  4:161683/138‑1 (MQ=255)                                                                                                              gCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGATTTCTGAAGATGTTTTTTCATCCATCGAGATg                         <  8:8930/139‑1 (MQ=255)                                                                                                              gCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGATTTCTGAAGATGTTTTTTCATCCATCGAGATg                         <  5:163065/139‑1 (MQ=255)                                                                                                                 tGCACCCGGCACCCCGGCGCGGGATGTGGGAATATAAGGGTTTTTAAGCAGGCTGATCATTGTTAGCCCCACTACCGGCTTCGTGTTGACCCTGGCGGGAATGATTTCTGAAGATGTTTTTTAAACCAGAGAGaggagg                      >  2:36878/1‑139 (MQ=255)                                                                                                                  gCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGATTTCTGAAGATGTTTTTTCATCCATCGAGATGAg                       <  5:364049/137‑1 (MQ=255)                                                                                                                     cccGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGATTTCTGAAGATGTTTTTTCATCCATCGAGATGAGGGATa                  <  6:280954/139‑1 (MQ=255)                                                                                                                            ccccGGCGCGGGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGATTTCTGAAGATGTTTTTTCATCCATCGAGATGAGGGATAAAATATg           <  2:216845/139‑1 (MQ=255)                                                                                                                            ccccGGCGAGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGATTTCTGAAGATGTTTTTTCATCCATCGAGATGAGGGATAAAATATg           <  1:59248/139‑1 (MQ=255)                                                                                                                             cccGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGATTTCTGAAGATGTTTTTTCATCCATCGAGATGAGGGATAAAATATGa          <  1:117510/139‑1 (MQ=255)                                                                                                                                ggCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTATTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGATTTATGAATATTTTTTTTCATCCATCGAGATGAGGGATAAAATATGAAtt       >  8:109358/1‑138 (MQ=255)                                                                                                                                  cgcgTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGATTTCTGAAGATGTTTTTTCATCCATCGAGATGAGGGATAAAATATGAATGGc     <  5:92566/139‑1 (MQ=255)                                                                                                                                     gTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGATTTCTGAAGATGTTTTTTCATCCATCGAGATGAGGGATAAAATATGAATGGCAAc  <  8:60711/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                          CTACCGGACATTTCGGATTCGGGTTTTCATGAACCGCAAAGAGTTCTTCTTGCTTTTGGACAGCCCGGTACACTTCTAAAAGAGAAATATCAGCAGGATCTTTTTTGAGGCTTGCACCCGGCACCCCGGCGCGTGATGTGAGAATATCAGCTTTTTTCAGCAGGCTGATCATTCTTCGCACAACTACCGGATTCGTGTTGACACTGTCTGCAATGATTTCTGAAGATGTTTTTTCATCCATCGAGATGAGGGATAAAATATGAATGGCAAC  >  BS168_reference_genome/3766288‑3766558

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