Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,591,144	T→C	N158S (AAT→AGT)	ydeH ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,591,144	0	T	C	100.0%	205.9 / NA	49	N158S (AAT→AGT)	ydeH	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/26);  total (23/26)

ATCGCCGGTTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGCCCCGGCAATCCCGTCA  >  REL606/1591071‑1591216                                                                          |                                                                         aTCGCCGGTTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTc                                                                        >  1:1303998/1‑76 (MQ=255)  tCGCCGGTTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCa                                                                       <  1:1088208/76‑1 (MQ=255)    gCCGGTTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCaga                                                                     <  1:1365825/76‑1 (MQ=255)     ccGGTTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCagag                                                                    <  2:650461/76‑1 (MQ=255)      cGGTTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAgg                                                                   >  2:1146538/1‑76 (MQ=255)       ggTTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGc                                                                  <  1:1268908/76‑1 (MQ=255)       ggTTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGc                                                                  <  1:505465/76‑1 (MQ=255)       ggTTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGc                                                                  >  1:158436/1‑76 (MQ=255)         ttAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTc                                                                <  2:699879/76‑1 (MQ=255)         ttAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTc                                                                <  1:530404/76‑1 (MQ=255)         ttAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTc                                                                <  2:1206962/76‑1 (MQ=255)         ttAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTc                                                                <  2:854693/76‑1 (MQ=255)          tAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCa                                                               >  2:1184114/1‑76 (MQ=255)          tAAATGCCCGTAGGTATCATTAACCAATTTAAATCGATCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCa                                                               <  1:329173/76‑1 (MQ=255)           aaaTGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAg                                                              <  2:834500/76‑1 (MQ=255)               gCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGtt                                                          <  1:1232200/76‑1 (MQ=255)               gCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGtt                                                          >  1:994058/1‑76 (MQ=255)                cccGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTg                                                         >  2:421762/1‑76 (MQ=255)                   gTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGt                                                      >  2:372150/1‑76 (MQ=255)                    tAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTa                                                     >  2:1635182/1‑76 (MQ=255)                     aGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTaa                                                    >  1:536447/1‑76 (MQ=255)                     aGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTaa                                                    <  2:716014/76‑1 (MQ=255)                     aGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTaa                                                    <  2:1320926/76‑1 (MQ=255)                        tATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACtg                                                 <  2:179645/76‑1 (MQ=255)                           cATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACtgatg                                              <  2:485930/76‑1 (MQ=255)                             ttAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACTgatgat                                            <  2:538245/76‑1 (MQ=255)                                aCCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCaa                                         <  1:400521/76‑1 (MQ=255)                                   aaTTTAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAgg                                      <  1:678669/76‑1 (MQ=255)                                     tttAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGAt                                    >  1:53695/1‑76 (MQ=255)                                      ttAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGAtt                                   <  1:1164051/76‑1 (MQ=255)                                      ttAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGAtt                                   >  1:1281601/1‑76 (MQ=255)                                       tAAATCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATtc                                  <  1:830694/76‑1 (MQ=255)                                        aaaTCGGTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATtca                                 <  1:120256/76‑1 (MQ=255)                                              gTCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGa                           <  1:595033/76‑1 (MQ=255)                                               tCAATATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGaa                          >  2:830607/1‑76 (MQ=255)                                                 aaTATCCAACAACATTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAAcc                        >  2:1324707/1‑76 (MQ=255)                                                          caacaTTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGcccc               >  1:612566/1‑76 (MQ=255)                                                          caacaTAAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGcccc               <  1:943219/76‑1 (MQ=255)                                                           aacaTTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGCCCCg              <  2:994986/76‑1 (MQ=255)                                                            acaTTAAATAAAGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGCCCCgg             >  1:1663165/1‑76 (MQ=255)                                                                 aaataaaGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGCCCCGGCAATc        >  2:964474/1‑76 (MQ=255)                                                                 aaataaaGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGCCCCGGCAATc        >  2:1432161/1‑76 (MQ=255)                                                                  aataaaGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGCCCCGGCAATcc       >  2:1359825/1‑76 (MQ=255)                                                                  aataaaGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGCCCCGGCAATcc       >  2:158654/1‑76 (MQ=255)                                                                      aaGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGCCCCGGCAATCCCGTc   >  2:692801/1‑76 (MQ=255)                                                                      aaGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGCCCCGGCAATCCCGTc   >  1:1342987/1‑76 (MQ=255)                                                                      aaGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGCCCCGGCAATCCCGTc   >  2:211170/1‑76 (MQ=255)                                                                      aaGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGCCCCGGCAATCCCGTc   >  2:974018/1‑76 (MQ=255)                                                                       aGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGCCCCGGCAATCCCGTCa  <  2:639348/76‑1 (MQ=255)                                                                       aGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGCCCCGGCAATCCCGTCa  >  1:984201/1‑76 (MQ=255)                                                                       aGACTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGCCCCGGCAATCCCGTCa  <  2:735449/76‑1 (MQ=255)                                                                          |                                                                         ATCGCCGGTTAAATGCCCGTAGGTATCATTAACCAATTTAAATCGGTCAATATCCAACAACATTAAATAAAGATTCAGAGGCTCAGCGTTGCGTAACTGATGATCAAAGGATTCATCAAGAACCCGACGCCCCGGCAATCCCGTCA  >  REL606/1591071‑1591216

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