Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	2,436,674	A→T	6.7%	Q210H (CAA→CAT)	ypuA →	putative exported protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	2,436,674	0	A	T	6.7%	171.8 / 3.7	60	Q210H (CAA→CAT)	ypuA	putative exported protein
Reads supporting (aligned to +/- strand):  ref base A (17/39);  new base T (2/2);  total (19/41)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

CCAAGAGCTTGTCACGACATCAGAGCTTGGAGACAAAATCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACATTGACTGGGGACAGGTCAGTGACCAGC  >  BS168_reference_genome/2436542‑2436795                                                                                                                                     |                                                                                                                          ccaaGAGCTTGTCACGACATCAGAGCTTGGAGACAAAATCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGAt                                                                                                                     <  1:51054/139‑1 (MQ=255)   aaGAGCTTGTCACGACATCTGTCCTTGGTGACAAAATCGGTAACGAAATTGCGGCAGCCCTGATTGAAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGa                                                                                                                   <  2:160977/139‑1 (MQ=255)      agCTTGTCACGACATCAGAGCTTGGAGACAAAATCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGc                                                                                                                <  5:74041/139‑1 (MQ=255)           gTCACGACATCAGATCTTGGAGACAAAATCGGTAACGAAAATGCGGCAGCCATGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAAAGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCtt                                                                                                           >  7:29522/1‑139 (MQ=255)           gTCACGACATCAGAGCTTGGAGACAAAATCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCtt                                                                                                           >  5:405664/1‑139 (MQ=255)           gTCACGACATCAGAGCTTGGAGACAAAATCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCtt                                                                                                           >  5:405669/1‑139 (MQ=255)           gTCACGACATCAGAGCTTGGAGACAAAATCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCtt                                                                                                           >  7:278894/1‑139 (MQ=255)            tCACGACATCAGAGCTTGGAGACAAAATCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACATGTGGATGACGCAGCTTc                                                                                                          >  7:60839/1‑139 (MQ=255)            tCACGACATCAGAGCTTGGAGACAAAATCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACATGTGGATGACGCAGCTTc                                                                                                          >  7:67471/1‑139 (MQ=255)            tCACGACATCAGAGCTTGGAGAAAAAATCGTTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGTGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTc                                                                                                          <  6:347703/139‑1 (MQ=255)                     cAGAGCTTGGAGACAAAATCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTaaa                                                                                                 >  1:136891/1‑139 (MQ=255)                      agagCTTGGAGACAAAATCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAAc                                                                                                >  3:240748/1‑139 (MQ=255)                        agCTTGGAGACAAAATCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGt                                                                                              <  4:374535/139‑1 (MQ=255)                        agCTTGGAGACAAAATCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGGGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGt                                                                                              <  2:152401/139‑1 (MQ=255)                          cTTGGAGACAAAATCCGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTca                                                                                            <  8:401684/139‑1 (MQ=255)                              gagaCAAAATCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCacac                                                                                         <  7:96193/138‑1 (MQ=255)                                    aaaTCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGATTTTGCTAAAAACGGGGGGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTacaga                                                                                  <  6:183116/139‑1 (MQ=255)                                       tCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTacagacag                                                                               <  6:355235/139‑1 (MQ=255)                                           tcaCGAAAATGCGGCAGCCCTGATTGCAAAAAACCAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCaa                                                                           <  6:318444/137‑1 (MQ=255)                                           tAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCa                                                                            >  1:421447/1‑138 (MQ=255)                                           gaaCGAAAAAGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGGGCCTGAATATAAAGCTGACATCGAAAAACAACTTGAAGACGCAGCTTCAGAATTAAACCTCACACTTACAGAGAGGCaa                                                                           <  6:324293/138‑1 (MQ=255)                                                aaaaTGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGaa                                                                      <  1:273116/139‑1 (MQ=255)                                                   ctGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCa                                                                   <  2:31267/138‑1 (MQ=255)                                                         cAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGt                                                             >  4:11558/1‑139 (MQ=255)                                                         cAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGt                                                             >  4:174978/1‑139 (MQ=255)                                                           gCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACTGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGAAGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGgct                                                           <  3:400110/139‑1 (MQ=255)                                                             ccTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGCCATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTcc                                                         <  6:405669/139‑1 (MQ=255)                                                             ccTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTcc                                                         <  6:405664/139‑1 (MQ=255)                                                             ccTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTcc                                                         <  2:136891/139‑1 (MQ=255)                                                                      aaaaaTCAAAGTAGAGTTTGCTAAAAACGGGCTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAAt                                                <  3:240498/139‑1 (MQ=255)                                                                      aaaaaTCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCaa                                                 >  8:414487/1‑138 (MQ=255)                                                                      aaaaaTCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCaa                                                 >  8:420106/1‑138 (MQ=255)                                                                      aaaaaTCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAAt                                                <  5:131487/139‑1 (MQ=255)                                                                              aagaagAGTTTGCTAACAACGGGGTTCCTGATAATAAAGCTGACATCGAAAAACAATTGGATGACGCAGCTTCTGATTTAAACGTAACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGaa                                        <  1:119167/139‑1 (MQ=255)                                                                              aagaagAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGaa                                        <  1:119168/139‑1 (MQ=255)                                                                                   gagTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACATGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACg                                   <  8:67471/139‑1 (MQ=255)                                                                                   gagTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACATGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACg                                   <  8:60839/139‑1 (MQ=255)                                                                                   gagTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTTTTCAATAAGATGAAAAACg                                   <  5:170722/139‑1 (MQ=255)                                                                                   gagTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACg                                   <  5:350832/139‑1 (MQ=255)                                                                                          cTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACAt                            <  3:11558/139‑1 (MQ=255)                                                                                          cTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACAt                            <  4:240748/139‑1 (MQ=255)                                                                                               aaCGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACATTGACt                       <  2:421447/139‑1 (MQ=255)                                                                                               aaCGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACATTGACt                       <  3:24451/139‑1 (MQ=255)                                                                                                    ggTGCCTGATACTTAAGATGAGATCGTAAATATACTGGATGTCGCAGCATGTTTTTTAATCGTCTCTCATACAGTCTTGCAAAAGTATCTGGTTTTCTCCCTTATCAATATGATGCAAAACGCTGACATTGACTGGGGa                  <  4:425530/139‑1 (MQ=255)                                                                                                    ggTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACATTGACTGGGGa                  <  4:370200/139‑1 (MQ=255)                                                                                                    ggTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACATTGACTGGGGa                  <  8:68588/139‑1 (MQ=255)                                                                                                    ggTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACATTGACTGGGGa                  >  1:284386/1‑139 (MQ=255)                                                                                                    ggTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACATTGACTGGGGa                  <  7:99101/139‑1 (MQ=255)                                                                                                               ataaAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACATTGACTGGGGACAGGTCAGTg        <  3:174978/138‑1 (MQ=255)                                                                                                                taaAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACATTGACTGGGGACAGGTCAGTGAc      <  6:259081/139‑1 (MQ=255)                                                                                                                taaAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACATTGACTGGGGACAGGTCAGTGAc      <  8:354939/139‑1 (MQ=255)                                                                                                                   aGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTg                                                              >  6:40805/1‑80 (MQ=255)                                                                                                                   aGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTg                                                              <  5:40805/80‑1 (MQ=255)                                                                                                                        aCATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGtct                                                           <  6:376643/78‑1 (MQ=255)                                                                                                                        aCATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGtct                                                           >  5:376643/1‑78 (MQ=255)                                                                                                                        aCATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTa                                                      <  5:2213/83‑1 (MQ=255)                                                                                                                        aCATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTa                                                      >  6:2213/1‑83 (MQ=255)                                                                                                                        aCATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACATTGACTGGGGACAGGt             >  7:276551/1‑124 (MQ=255)                                                                                                                        aCATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACATTGACTGGGGACAGGt             <  8:276551/124‑1 (MQ=255)                                                                                                                        aCATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACATTGACTGGGGACAGGTCAGTGACCAGc  >  7:218429/1‑135 (MQ=255)                                                                                                                                     |                                                                                                                          CCAAGAGCTTGTCACGACATCAGAGCTTGGAGACAAAATCGGTAACGAAAATGCGGCAGCCCTGATTGCAAAAATCAAAGAAGAGTTTGCTAAAAACGGGGTGCCTGATAATAAAGCTGACATCGAAAAACAAGTGGATGACGCAGCTTCTGATTTAAACGTCACACTTACAGACAGCCAAAAGAATCAGCTTGTCTCCCTATTCAATAAGATGAAAAACGCTGACATTGACTGGGGACAGGTCAGTGACCAGC  >  BS168_reference_genome/2436542‑2436795

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