Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	401,700	T→C	S93P (TCG→CCG)	ybaD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	401,700	0	T	C	98.1%	200.2 / ‑4.6	52	S93P (TCG→CCG)	ybaD	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (23/28);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

GGAATGCTGCGGGCGCTGGAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAATCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAGATGATTGGCAATCTGGTGATGGAGCAATTGA  >  REL606/401628‑401772                                                                         |                                                                         ggAATGCTGCGGGCGCTGGAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGc                                                                       <  2:1027348/76‑1 (MQ=255) ggAATGCTGCGGGCGCTGGAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGc                                                                       <  2:113181/76‑1 (MQ=255)   aaTGCTGCGGGCGCTGGAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAg                                                                     <  2:563020/76‑1 (MQ=255)     tgctgcGGGCGCTGGAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCt                                                                   >  1:470285/1‑76 (MQ=255)       ctgcGGGCGCTGGAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTgc                                                                 >  1:1179171/1‑76 (MQ=255)       ctgcGGGCGCTGGAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTgc                                                                 >  1:1426889/1‑76 (MQ=255)         gcGGGCGCTGGAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTgcgc                                                               >  1:616312/1‑76 (MQ=255)         gcGGGCGCTGGAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTgcgc                                                               >  2:1820431/1‑76 (MQ=255)          cGGGCGCTGGAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTgcgcg                                                              <  1:695737/76‑1 (MQ=255)          cGGGCGCTGGAAAAACGTCCGGTGAGTTCCGATGACGTAGAAATGGCAATCAATCATATTAAACCGCAGCTgcgcg                                                              <  1:348371/76‑1 (MQ=255)             gcgcTGGAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCa                                                           <  2:205990/76‑1 (MQ=255)              cgcTGGAAAAACGTCCGGTGAGTTCCGATGACGTCGATATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCAc                                                          >  2:1064107/1‑76 (MQ=255)              cgcTGGAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCAc                                                          >  2:1268793/1‑76 (MQ=255)                  ggAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGt                                                      >  1:1111359/1‑76 (MQ=255)                  ggAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGt                                                      <  1:478410/76‑1 (MQ=255)                  ggAAAAACGTCCGGTGAGTTCCGACGACNCCGAAATGGCAATCAATCACANTACACCGCAGCTGCGCGCCACCggc                                                      >  1:288978/1‑75 (MQ=255)                   gAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTg                                                     <  2:1320076/76‑1 (MQ=255)                     aaaaCGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAg                                                   >  2:1304832/1‑76 (MQ=255)                          gTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGa                                              <  2:413477/76‑1 (MQ=255)                          gTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGa                                              <  2:471634/76‑1 (MQ=255)                           tCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGaa                                             >  2:1915064/1‑76 (MQ=255)                             cGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGt                                           >  1:1023789/1‑76 (MQ=255)                               gTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGc                                         <  2:894156/76‑1 (MQ=255)                                 gAGTTGCGATGGCGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCg                                       <  2:1143585/76‑1 (MQ=255)                                 gAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCg                                       <  2:213467/76‑1 (MQ=255)                                    ttCCGATGACGTCGAAATGGCAATCAATCATATTAAACTGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGc                                    <  1:1698007/76‑1 (MQ=255)                                    ttCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGc                                    >  1:1931347/1‑76 (MQ=255)                                     tCCGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCa                                   <  2:1544902/76‑1 (MQ=255)                                       cGNTGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGTCGAGGNAg                                 <  2:1665413/76‑1 (MQ=255)                                       cGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCaac                                 >  2:1745305/1‑75 (MQ=255)                                       cGATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAg                                 >  1:405949/1‑76 (MQ=255)                                        gATGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAga                                >  1:673966/1‑76 (MQ=255)                                         aTGACGTCGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAgat                               <  2:1023788/76‑1 (MQ=255)                                          tggcGTCGAAATGGCAATAAATAATATNGANCTNCATCTGCNCGGCACTGGTGAGCGCGAAGTGCCGAGCAAgatg                              <  1:1762162/73‑1 (MQ=255)                                           gACGTCGAAATGGCAATCCATNATNTTAAACCTCAGCTGCGCGCCANCGGTGAGNGNGAAGTGCCGAGCAAgatga                             >  2:298795/1‑76 (MQ=255)                                               tCGAAATGGCANTCAATCATATTAAACCGCAGCTTCACGTCACGNGTGAGCGCGAAGTGCCGCGGACGATGATTgg                         <  1:245338/76‑1 (MQ=255)                                                cGAAATGGCAATCAATCATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAGATGATTGGc                        <  1:1263892/76‑1 (MQ=255)                                                 gAAATGGCAATCAATCATATTAAATCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAGATGATTGGCa                       <  2:1424379/76‑1 (MQ=255)                                                 gAAATGGCAAACAGTCATATNAAACCGCAGTTGCGCNCCACCGGTGAGCGCGCAGTGCCGAGCAAGATGATTGGCa                       <  2:2015257/76‑1 (MQ=255)                                                       gCAATCAATCATATTTAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAGATGATTGGCAATCTgg                 <  1:1691917/76‑1 (MQ=255)                                                         aatcaatcATATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAGATGATTGGCAATCTGGTg               <  1:1521964/76‑1 (MQ=255)                                                          atcaatcaTATTAAACCGCAGCTGCGCGCCACCGGTGAGCTCGAAGTGCCGAGCAAGATGATTGGCACTCTGGTGa              <  2:42338/76‑1 (MQ=255)                                                            caatcaTATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAGATGATTGGCAATCTGGTGATg            <  2:1288056/76‑1 (MQ=255)                                                             aatcaTATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAGATGATTGGCAATCTGGTGATgg           <  2:1728866/76‑1 (MQ=255)                                                             aatcaTATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAGATGATTGGCAATCTGGTGATgg           <  2:1370608/76‑1 (MQ=255)                                                                caTATTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAGATGATTGGCAATCTGGTGATGGAGc        >  2:1547954/1‑76 (MQ=255)                                                                 atatTAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAGATGATTGGCAATCTGGTGATGGAGCa       >  2:1810352/1‑76 (MQ=255)                                                                 atatTAAACCGCAGCTGCGCGCCACCGGGGAGCGCGAAGTGCCGAGCAAGATGATTGGCAATCTGGGGATGGGGCa       >  1:454568/1‑76 (MQ=255)                                                                     tAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAGATGATTGGCAATCTGGTGATGGAGCAATTg   >  2:1712288/1‑76 (MQ=255)                                                                     tAAACCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAGATGATTGGCAATCTGGTGATGGAGCAATTg   <  2:1001669/76‑1 (MQ=255)                                                                      aaaCCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGc                                    >  1:434413/1‑42 (MQ=255)                                                                      aaaCCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGc                                    <  2:434412/42‑1 (MQ=255)                                                                      aaaCCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAGATGATTGGCAATCTGGTGATGGAGCAATTGa  >  2:356324/1‑76 (MQ=255)                                                                      aaaCCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAGATGATTGGCAATCTGGTGATGGAGCAATTGa  <  1:1272264/76‑1 (MQ=255)                                                                         |                                                                         GGAATGCTGCGGGCGCTGGAAAAACGTCCGGTGAGTTCCGATGACGTCGAAATGGCAATCAATCATATTAAATCGCAGCTGCGCGCCACCGGTGAGCGCGAAGTGCCGAGCAAGATGATTGGCAATCTGGTGATGGAGCAATTGA  >  REL606/401628‑401772

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑