Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,902,276	(T)7→8	coding (7/342 nt)	yebY ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,902,269	1	.	T	96.2%	187.8 / 1.3	52	coding (14/342 nt)	yebY	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (23/27);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CGCCCATTTGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAGT  >  REL606/1902169‑1902364                                                                                                      |                                                                                                cGCCCATTTGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAc                                                                                                  >  1:1195528/1‑101 (MQ=255)   cccATTTGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑tt                                                                                               <  1:517383/101‑1 (MQ=255)   cccATTTGTCTTTACCCACTTCAAAACGGATGACGGGAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGCAACGCCAGAATAC‑tt                                                                                               >  1:923970/1‑101 (MQ=255)    ccATTTGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCGAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑ttt                                                                                              >  1:1103816/1‑101 (MQ=255)     cATTTGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑tttt                                                                                             <  1:379235/101‑1 (MQ=255)      aTTTGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑ttttt                                                                                            <  1:483605/101‑1 (MQ=255)        ttGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑ttttttt                                                                                          <  1:906466/101‑1 (MQ=255)         tGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑tttttttt                                                                                         <  2:1147605/101‑2 (MQ=255)         tGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑tttttttt                                                                                         >  1:1407051/1‑100 (MQ=255)         tGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑tttttttt                                                                                         <  1:716945/101‑2 (MQ=255)            cTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑ttttttttcat                                                                                      <  1:1056600/101‑5 (MQ=255)              ttACCCACTTCAAAACGGCTGACGGTAATCCCCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTtcatca                                                                                     <  2:1542924/101‑1 (MQ=255)              ttACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTtcatca                                                                                     <  1:27412/101‑1 (MQ=255)                 cccACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCAttt                                                                                  >  2:1166000/1‑101 (MQ=255)                  ccACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCAtttt                                                                                 <  2:286068/101‑1 (MQ=255)                    aCTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCg                                                                               >  1:227700/1‑101 (MQ=255)                    aCTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCg                                                                               <  1:1337707/101‑1 (MQ=255)                     cTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGt                                                                              <  1:414889/101‑1 (MQ=255)                       tCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTaa                                                                            <  1:994988/101‑1 (MQ=255)                        cAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAAt                                                                           <  2:1344734/101‑1 (MQ=255)                            aCGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCAt                                                                       <  1:1585221/101‑1 (MQ=255)                             cGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATc                                                                      >  2:1332840/1‑101 (MQ=255)                                 tGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCaaaa                                                                  >  2:1598922/1‑101 (MQ=255)                                   aCGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAg                                                                >  1:980014/1‑101 (MQ=255)                                   aCGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAg                                                                <  1:829714/101‑1 (MQ=255)                                      gTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGtt                                                             <  2:355725/101‑1 (MQ=255)                                      gTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGtt                                                             >  1:426650/1‑101 (MQ=255)                                       tAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCg                                                           >  2:305073/1‑101 (MQ=255)                                        aaTCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCg                                                           >  2:178559/1‑101 (MQ=255)                                         aTCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCga                                                          <  2:821083/101‑1 (MQ=255)                                              cTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTc                                                     >  1:171801/1‑101 (MQ=255)                                                   gTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAg                                                <  1:725259/101‑1 (MQ=255)                                                    tGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGc                                               <  2:980009/101‑1 (MQ=255)                                                       cgccAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAgg                                            <  1:1189290/101‑1 (MQ=255)                                                            gcgcCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCtgat                                       <  2:1103811/101‑1 (MQ=255)                                                            gcgcCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCtgat                                       <  2:153829/101‑1 (MQ=255)                                                             cgccgcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCtgatg                                      >  1:86747/1‑101 (MQ=255)                                                                 gcTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCa                                  <  1:2522/101‑1 (MQ=255)                                                                  cTGCAGAACTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCAc                                 >  1:1493691/1‑101 (MQ=255)                                                                   tGCAGAACTGGTGAGTAACAGAAACGCCAGAATTCTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCAcc                                <  1:1696077/101‑1 (MQ=255)                                                                        aaCTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGtac                           <  2:1098433/101‑1 (MQ=255)                                                                         aCTGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGtact                          <  1:987247/101‑1 (MQ=255)                                                                          cTGGTGAGTAACAGAAACGCCAGAATAC‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGtactac                        <  2:1281154/101‑1 (MQ=255)                                                                           tGGTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGtagtac                        >  1:578681/1‑101 (MQ=255)                                                                            ggTGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTAcc                       >  2:658296/1‑101 (MQ=255)                                                                              tGAGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGc                     >  2:528126/1‑101 (MQ=255)                                                                                aGTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCAc                   <  2:1188303/101‑1 (MQ=255)                                                                                 gTAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCAcc                  >  2:1136439/1‑101 (MQ=255)                                                                                  tAACAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCt                 >  2:1586441/1‑101 (MQ=255)                                                                                     cAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATc              <  2:123867/101‑1 (MQ=255)                                                                                     cAGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATc              <  1:296487/101‑1 (MQ=255)                                                                                      aGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCg             >  1:255317/1‑101 (MQ=255)                                                                                      aGAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCg             <  2:251888/101‑1 (MQ=255)                                                                                       gAAACGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGc            <  2:536137/101‑1 (MQ=255)                                                                                           cGCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCAt        <  1:1166005/101‑1 (MQ=255)                                                                                            gCCAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCAtt       >  1:1285953/1‑101 (MQ=255)                                                                                             ccAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTc      >  2:1457578/1‑101 (MQ=255)                                                                                             ccAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTc      >  2:292511/1‑101 (MQ=255)                                                                                              cAGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCg     <  1:1203638/101‑1 (MQ=255)                                                                                               aGAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGa    >  1:1672125/1‑101 (MQ=255)                                                                                                gAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAg   >  1:960800/1‑101 (MQ=255)                                                                                                gAATACTTTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCACAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAg   >  2:67495/1‑101 (MQ=255)                                                                                                 aaTACTTTTTTTTCATCATTTTCGTAATCCTTCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAGt  >  1:67747/1‑101 (MQ=255)                                                                                                      |                                                                                                CGCCCATTTGTCTTTACCCACTTCAAAACGGCTGACGGTAATCACCTGCGGTGCCGCCAGCGCCGCTGCAGAACTGGTGAGTAACAGAAACGCCAGAATAC‑TTTTTTTCATCATTTTCGTAATCCATCAAAAAGGTTCGAGAGTCGCAAGCAGGCTGATGATCACCAGTACTACCGCACCTATCGCCCATTCGAGT  >  REL606/1902169‑1902364

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