Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,930,575	C→A	F137L (TTC→TTA)	yecE →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,930,575	0	C	A	97.1%	244.1 / 0.4	69	F137L (TTC→TTA)	yecE	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (37/30);  total (38/31)
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.64e-02

ACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTCGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTC  >  REL606/1930477‑1930672                                                                                                   |                                                                                                   acaTTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGc                                                                                                  >  2:1147361/1‑101 (MQ=255)         cccACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGgag                                                                                          <  2:939137/101‑1 (MQ=255)         cccACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGgag                                                                                          <  1:1321206/101‑1 (MQ=255)          ccACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGgaga                                                                                         >  1:420246/1‑101 (MQ=255)              gggAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAagag                                                                                     >  2:1224493/1‑101 (MQ=255)               ggAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAgg                                                                                    >  2:1483406/1‑101 (MQ=255)                gAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGa                                                                                   <  1:258945/101‑1 (MQ=255)                     gcctgcGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTCGCCAAAGGAGAAGAGGAACaaa                                                                              <  2:622664/101‑1 (MQ=255)                     gcctgcGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACaaa                                                                              <  2:310385/101‑1 (MQ=255)                       ctgcGCTTTGGCATTTTCTCGATTTTCTTCCTGGTGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACg                                                                            <  1:91210/101‑1 (MQ=255)                       ctgcGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACg                                                                            >  2:435172/1‑101 (MQ=255)                        tgcGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGc                                                                           >  1:1228478/1‑101 (MQ=255)                          cgcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCtt                                                                         <  1:384102/101‑1 (MQ=255)                           gcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTa                                                                        <  2:1569673/101‑1 (MQ=255)                           gcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCCCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTa                                                                        <  2:28146/101‑1 (MQ=255)                            cTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTaa                                                                       >  1:541857/1‑101 (MQ=255)                             tttGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAAt                                                                      >  2:835952/1‑101 (MQ=255)                             tttGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAAt                                                                      <  1:612272/101‑1 (MQ=255)                              ttGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATc                                                                     >  2:1209219/1‑101 (MQ=255)                               tGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATcg                                                                    >  2:418289/1‑101 (MQ=255)                                   aTTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGt                                                                >  2:505121/1‑101 (MQ=255)                                   aTTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGt                                                                >  2:990047/1‑101 (MQ=255)                                      ttCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTa                                                             >  2:252946/1‑101 (MQ=255)                                      ttCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTa                                                             <  1:176525/101‑1 (MQ=255)                                       tctcGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTAc                                                            <  1:207754/101‑1 (MQ=255)                                        ctcGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACa                                                           >  1:528383/1‑101 (MQ=255)                                         tcGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACAt                                                          >  2:397843/1‑101 (MQ=255)                                         tcGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACAt                                                          <  2:889103/101‑1 (MQ=255)                                         tcGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACAt                                                          <  1:941108/101‑1 (MQ=255)                                          cGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATc                                                         <  1:1240288/101‑1 (MQ=255)                                           gATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCa                                                        >  1:1560809/1‑101 (MQ=255)                                           gATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCa                                                        >  2:1487873/1‑101 (MQ=255)                                           gATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCa                                                        >  2:981814/1‑101 (MQ=255)                                            aTTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAg                                                       <  2:1159287/101‑1 (MQ=255)                                            aTTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAg                                                       <  2:855354/101‑1 (MQ=255)                                              tctcTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAgcg                                                     >  1:151015/1‑101 (MQ=255)                                               ctctTCCTGGCGAATTTAATTATGGCGTGGAAGGCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAgcgc                                                    >  1:1445532/1‑101 (MQ=255)                                                tctcCCTGGCGAATTTAATTATGGCGTGGAAGGCTGCCAGCCACAGGTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGcgcg                                                   <  2:585659/101‑1 (MQ=255)                                                   tCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCgcggcg                                                <  1:1451652/101‑1 (MQ=255)                                                     tcttccgatctAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGtt                                              <  1:498741/91‑1 (MQ=255)                                                      tGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTa                                             >  1:1016301/1‑101 (MQ=255)                                                         cGAATTTAATTAGGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATc                                         <  1:39671/101‑1 (MQ=255)                                                          gAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCg                                        <  2:503771/101‑1 (MQ=255)                                                          gAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCg                                        >  2:1364903/1‑101 (MQ=255)                                                            aTTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCggg                                      >  1:291807/1‑101 (MQ=255)                                                                aaTTATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTAGATCggaaga                                   >  2:498736/1‑97 (MQ=255)                                                                  ttATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGAttt                                >  2:521050/1‑101 (MQ=255)                                                                   tATGGCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGAtttt                               >  1:1709410/1‑101 (MQ=255)                                                                       gCGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGAc                           <  1:1487879/101‑1 (MQ=255)                                                                        cGTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACa                          <  1:844515/101‑1 (MQ=255)                                                                         gTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTG‑ATCGGGTGATTTTAGACac                         >  2:445773/1‑100 (MQ=255)                                                                         gTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAg                         >  2:589709/1‑101 (MQ=255)                                                                         gTGGAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAg                         >  1:1222821/1‑101 (MQ=255)                                                                           ggtaGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAgcc                       >  1:280037/4‑101 (MQ=255)                                                                           gaaaGCCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAgcc                       <  2:1660011/99‑1 (MQ=255)                                                                            gAAGTCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAgccg                      <  1:1147366/101‑1 (MQ=255)                                                                             aaGTCCGCCATCCACAGTTTTTCGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAgccgc                     >  1:584744/1‑101 (MQ=255)                                                                                tCCGCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCg                  >  2:742750/1‑101 (MQ=255)                                                                                   gCCATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCGGtt               >  2:32715/1‑101 (MQ=255)                                                                                    ccATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCGGTTc              <  2:1016296/101‑1 (MQ=255)                                                                                     cATCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCGGTTCa             >  2:861562/1‑101 (MQ=255)                                                                                      aTCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCGGTTCAt            <  1:639765/101‑1 (MQ=255)                                                                                      aTCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCGGTTCAt            <  2:1323347/101‑1 (MQ=255)                                                                                      aTCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCGGTTCAt            >  2:1355594/1‑101 (MQ=255)                                                                                      aTCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCGGTTCAt            >  1:1196118/1‑101 (MQ=255)                                                                                      aTCCACAGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCGGTTCAt            >  2:229975/1‑101 (MQ=255)                                                                                         cacaGTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCGGTTCATgca         <  1:1278378/101‑1 (MQ=255)                                                                                             gTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCAc     <  1:981819/101‑1 (MQ=255)                                                                                             gTTTTTAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCAc     >  2:1445105/1‑101 (MQ=255)                                                                                               ttttAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGt   <  1:435177/101‑1 (MQ=255)                                                                                                tttAGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTc  >  2:320147/1‑101 (MQ=255)                                                                                                   |                                                                                                   ACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTTTCGCCAAAGGAGAAGAGGAACAAACGCTTAATCGCGGTTTACATCAGCGCGGCGTTA‑ATCGGGTGATTTTAGACAGCCGCCCGGTTCATGCAGCACGTC  >  REL606/1930477‑1930672

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