Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	2,478,071	C→T	M173I (ATG→ATA)	yfeY ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	2,478,071	0	C	T	93.0%	137.4 / 3.4	43	M173I (ATG→ATA)	yfeY	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base T (23/17);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 9.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-02

TGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTTGACTGCCTTCGGCTTTAC  >  REL606/2478000‑2478141                                                                        |                                                                       ttttCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTa                                                                    >  2:197952/3‑76 (MQ=255) tGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTa                                                                    >  1:406579/1‑76 (MQ=255)          gCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAg                                                           <  2:1332510/76‑1 (MQ=255)           cAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAgg                                                          <  1:726496/76‑1 (MQ=255)           cAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAgg                                                          <  1:1522083/76‑1 (MQ=255)           cAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAgg                                                          <  2:820642/76‑1 (MQ=255)            aaaTTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGa                                                         >  1:963433/1‑76 (MQ=255)            aaaTTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGa                                                         >  1:292221/1‑76 (MQ=255)            aaaTTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGa                                                         <  2:929026/76‑1 (MQ=255)                tAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCAc                                                     <  1:1437239/76‑1 (MQ=255)                tAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCAc                                                     >  1:666589/1‑76 (MQ=255)                 aaCGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACg                                                    >  2:1349256/1‑76 (MQ=255)                 aaCGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACg                                                    >  2:845344/1‑76 (MQ=255)                          aGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTcccc                                           >  1:481662/1‑76 (MQ=255)                          aGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTcccc                                           <  2:1399823/76‑1 (MQ=255)                           gATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCg                                          <  2:777718/76‑1 (MQ=255)                             tAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCt                                        <  2:1142958/76‑1 (MQ=255)                              aaTTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTg                                       >  2:764102/1‑76 (MQ=255)                              aaTTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTg                                       <  2:1286626/76‑1 (MQ=255)                                ttttACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGaa                                     <  1:456682/76‑1 (MQ=255)                                  ttACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACt                                   <  2:1190782/76‑1 (MQ=255)                                   tACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTg                                  >  1:1217983/1‑76 (MQ=255)                                      tGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTa                               >  2:1183300/1‑76 (MQ=255)                                      tGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTa                               >  2:927848/1‑76 (MQ=255)                                      tGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTa                               <  2:821679/76‑1 (MQ=255)                                      tGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTa                               <  2:821675/76‑1 (MQ=255)                                           ntCCAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTaa                          >  1:1129378/2‑76 (MQ=255)                                              cAGTTTTTCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAatat                       <  1:605781/76‑1 (MQ=255)                                                  ttttCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGACCTGGAAGCTAATATGTTg                   >  1:661405/1‑76 (MQ=255)                                                  ttttCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTTg                   >  1:661408/1‑76 (MQ=255)                                                  ttttCAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTTg                   <  2:13503/76‑1 (MQ=255)                                                   tttCGGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTTGa                  >  1:844743/1‑76 (MQ=255)                                                      cAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTTGACTg               >  1:976466/1‑76 (MQ=255)                                                      cAGGGTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTTGACTg               <  1:1427018/76‑1 (MQ=255)                                                          gTATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTTGACTGCCtt           <  2:406579/76‑1 (MQ=255)                                                           tATCGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTTGACTGCCTTc          >  1:940688/1‑76 (MQ=255)                                                              cGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGGTGAACTGGTAGCTAATATGTTGACTGCCTTCGGc       >  1:182770/1‑76 (MQ=255)                                                              cGTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTTGACTGCCTTCGGc       >  1:176892/1‑76 (MQ=255)                                                               gTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGCACTGGTAGCTAATATGTTGACTGCCTTCGGCt      >  2:97015/1‑76 (MQ=255)                                                               gTCCGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTTGACTGCCTTCGGCt      >  1:500268/1‑76 (MQ=255)                                                                 ccGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTTGACTGCCTTCGGCttt    <  2:658817/76‑1 (MQ=255)                                                                  cGAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTTGACTGCCTTCGGCTTTa   >  1:169262/1‑76 (MQ=255)                                                                   gAAGGTATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTTGACTGCCTTCGGCTTTAc  <  1:502163/76‑1 (MQ=255)                                                                        |                                                                       TGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTCCCCGCTGAACTGGTAGCTAATATGTTGACTGCCTTCGGCTTTAC  >  REL606/2478000‑2478141

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