Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,207,861	T→G	N47T (AAC→ACC)	ycfC ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,207,861	0	T	G	100.0%	167.8 / NA	43	N47T (AAC→ACC)	ycfC	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (25/18);  total (25/18)

GCGACTGCTGGCATTGAGCACGCCGAGCAAGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAAT  >  REL606/1207767‑1207958                                                                                               |                                                                                                  gCGACTGCTGGCATTGAGCACGCCGAGCAAGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCatat                                                                                             <  2:1335237/101‑1 (MQ=255)   gACTGCTGGCATTGAGCACGCCGAGCAAGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCa                                                                                           >  2:1541401/1‑101 (MQ=255)     ctgctgGCATTGAGCACGCCGAGCAAGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCaat                                                                                         <  1:972262/101‑1 (MQ=255)        ctgGCATTGAGCACGCCGAGCAAGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCaataat                                                                                      <  2:100007/101‑1 (MQ=255)         tgGCATTGAGCACGCCGAGCAAGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAataata                                                                                     <  2:1337338/101‑1 (MQ=255)           gCATTGAGCACGCCGAGCAAGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACt                                                                                   <  1:1237880/101‑1 (MQ=255)            cATTGAGCACGCCGAGCAAGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTg                                                                                  <  2:45468/101‑1 (MQ=255)               tGAGCACGCCGAGCAAGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTtg                                                                               <  2:395150/101‑1 (MQ=255)               tGAGCACGCCGAGCAAGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTtg                                                                               <  2:843862/101‑1 (MQ=255)                  gcccggggAGCAAGGTGTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCAAAAACCCGCCAGCGTCGAGCTCGGGGTCATATCAATAATACTGTtgagt                                                                            <  2:953041/94‑1 (MQ=255)                   cACGCCGAGCAAGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTtgcgtg                                                                           >  1:69762/1‑101 (MQ=255)                     cGCCGAGCAAGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTtgagtgag                                                                         <  1:109013/101‑1 (MQ=255)                     cGCCGAGCAAGGTTTCCAGCCCGACCCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCCTCGAGCTTGGGGTCATATCAATAATAATGTtgagttag                                                                         >  2:1231147/1‑101 (MQ=255)                             aaGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAg                                                                 >  1:466693/1‑101 (MQ=255)                             aaGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAg                                                                 >  1:780558/1‑101 (MQ=255)                               ggTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAgcg                                                               >  2:1489112/1‑101 (MQ=255)                                     cAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCgg                                                         >  2:159792/1‑101 (MQ=255)                                     cAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCgg                                                         >  2:348224/1‑101 (MQ=255)                                      aGCCAGACGCGCAGGTTGGCTTCGCCACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCAGc                                                        >  2:1699905/1‑101 (MQ=255)                                            aCGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATcaca                                                  >  1:1606986/1‑101 (MQ=255)                                                     ttGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTg                                         >  2:993495/1‑101 (MQ=255)                                                     ttGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTg                                         >  1:179129/1‑101 (MQ=255)                                                      tGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGGCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTgg                                        >  1:279748/1‑101 (MQ=255)                                                       ggCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTggg                                       >  1:154957/1‑100 (MQ=255)                                                       ggCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGt                                       >  1:894810/1‑101 (MQ=255)                                                           tCGCTACCGCCAAAAACCGCCAGCGTGGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGc                                   <  2:187642/101‑1 (MQ=255)                                                             gCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGa                                 <  1:427395/101‑1 (MQ=255)                                                                 ccGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAgttg                             >  1:288621/1‑101 (MQ=255)                                                                  cGCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAgttgt                            >  1:1434309/1‑101 (MQ=255)                                                                   gCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGTCCCTGGTGAGCGAgttgtt                           <  1:1488946/101‑1 (MQ=255)                                                                   gCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAgttgtt                           <  1:532944/101‑1 (MQ=255)                                                                   gCCAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAgttgtt                           <  2:1538332/101‑1 (MQ=255)                                                                     cAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGc                         >  1:977917/1‑101 (MQ=255)                                                                     cAAAAACCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGc                         >  1:93469/1‑101 (MQ=255)                                                                      aaaaaCCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGGTGTTGCa                        >  2:132450/1‑101 (MQ=255)                                                                       aaaaCCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCAc                       >  2:1126107/1‑101 (MQ=255)                                                                        aaaCCGCCAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCAcc                      >  2:65571/1‑101 (MQ=255)                                                                         acccGCCAGCGTCGAGCTGGGGGTCATATCCATAATTCTGTTGAGTGGGACGTGTAGCGCATCAGCATCACAATGCCCCTCGTGACCGAGTTGTTGCACCa                     <  2:1408786/99‑1 (MQ=255)                                                                              ccAGCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGt                >  2:955541/1‑101 (MQ=255)                                                                                 gCGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGcc             <  2:1155732/101‑1 (MQ=255)                                                                                  cGTCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCg            >  1:25778/1‑101 (MQ=255)                                                                                    tCGAGCTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGAc          <  2:997572/101‑1 (MQ=255)                                                                                         cTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTATCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACa     <  1:249054/101‑1 (MQ=255)                                                                                         cTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACa     >  1:891134/1‑101 (MQ=255)                                                                                         cTGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACa     >  2:125985/1‑101 (MQ=255)                                                                                          tGGGGGTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACaa    >  1:1477393/1‑101 (MQ=255)                                                                                            ggggTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAAt  <  2:977912/101‑1 (MQ=255)                                                                                               |                                                                                                  GCGACTGCTGGCATTGAGCACGCCGAGCAAGGTTTCCAGCCCGACGCGCAGGTTGGCTTCGCTACCGCCAAAAACCGCCAGCGTCGAGCTGGGGTTCATATCAATAATACTGTTGAGTGAGACGTGTAGCGCATCGGCATCACAATGCCCCTGGTGAGCGAGTTGTTGCACCAGGCGTGCCGACTGACAAAT  >  REL606/1207767‑1207958

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