Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	239,002	A→C	A4A (GCA→GCC)	dnaQ →	DNA polymerase III subunit epsilon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	239,002	0	A	C	100.0%	299.5 / NA	74	A4A (GCA→GCC)	dnaQ	DNA polymerase III subunit epsilon
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (38/36);  total (38/36)

AATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTT  >  REL606/238906‑239098                                                                                                 |                                                                                                 aaTTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTa                                                                                              >  2:21261/1‑101 (MQ=255) aaTTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTa                                                                                              <  1:504135/101‑1 (MQ=255)   tttCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTaca                                                                                            <  1:627885/101‑1 (MQ=255)     tCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACg                                                                                          <  1:222927/101‑1 (MQ=255)        aCCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCa                                                                                       <  2:1144091/101‑1 (MQ=255)         ccTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACGGCTATGAGCACTGCCATTACACGCCAg                                                                                      >  1:509466/1‑101 (MQ=255)         ccTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAg                                                                                      <  2:190341/101‑1 (MQ=255)          cTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGa                                                                                     >  1:1532792/1‑101 (MQ=255)           tGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGAt                                                                                    <  2:1651625/101‑1 (MQ=255)            gTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATc                                                                                   <  1:43602/101‑1 (MQ=255)              ttAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGt                                                                                 >  1:1568212/1‑101 (MQ=255)                aaGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTtc                                                                               >  1:728900/1‑101 (MQ=255)                 aGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTtct                                                                              >  1:1414557/1‑101 (MQ=255)                 aGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTtct                                                                              >  1:675211/1‑101 (MQ=255)                 aGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTtct                                                                              >  2:162486/1‑101 (MQ=255)                 aGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTtct                                                                              >  2:647158/1‑101 (MQ=255)                      ctctGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATa                                                                         >  1:620457/1‑101 (MQ=255)                         tGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCg                                                                      <  2:1491934/101‑1 (MQ=255)                             aGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAac                                                                  >  2:1430899/1‑101 (MQ=255)                             aGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAac                                                                  <  1:280554/101‑1 (MQ=255)                              gACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAacc                                                                 <  1:143435/101‑1 (MQ=255)                              gACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAacc                                                                 >  1:801777/1‑101 (MQ=255)                                  tCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCg                                                             >  2:1395891/1‑101 (MQ=255)                                    cTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                           <  1:804689/101‑1 (MQ=255)                                      gTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                         <  1:742610/101‑1 (MQ=255)                                      gTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                         >  1:325623/1‑101 (MQ=255)                                       tAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                        <  1:1590382/101‑1 (MQ=255)                                        aaTTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                       >  1:44553/1‑101 (MQ=255)                                         aTTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                      >  1:1386918/1‑101 (MQ=255)                                         aTTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACAGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                      <  1:369840/101‑1 (MQ=255)                                          ttGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                     >  1:1252788/1‑101 (MQ=255)                                           tGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAcc                                                    <  2:1674688/101‑1 (MQ=255)                                            gAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                   >  1:131433/1‑101 (MQ=255)                                                cGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                               <  2:1330939/101‑1 (MQ=255)                                                cGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                               >  1:995655/1‑101 (MQ=255)                                                   aCTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGt                                            <  1:607678/101‑1 (MQ=255)                                                    cTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTg                                           >  1:1649122/1‑101 (MQ=255)                                                     tGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgc                                          <  1:1432611/101‑1 (MQ=255)                                                      gTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                         <  2:1024265/101‑1 (MQ=255)                                                      gTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                         >  1:39311/1‑101 (MQ=255)                                                            cGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                   <  1:471039/101‑1 (MQ=255)                                                            cGATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                   <  2:1260248/101‑1 (MQ=255)                                                             gATAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                  >  2:195612/1‑101 (MQ=255)                                                               tAAGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                <  2:153039/101‑1 (MQ=255)                                                                aaGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                               >  1:917067/1‑101 (MQ=255)                                                                aaGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                               <  1:647163/101‑1 (MQ=255)                                                                 aGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATAGCGAAACCACCGGTATGAACCAGATTAGTGCGCACTATGGAgg                              >  1:87981/1‑101 (MQ=255)                                                                 aGTCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                              >  2:420962/1‑101 (MQ=255)                                                                   tCTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                            <  1:319228/101‑1 (MQ=255)                                                                    cTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                           >  2:947853/1‑101 (MQ=255)                                                                    cTGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                           >  2:1497790/1‑101 (MQ=255)                                                                     tGACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                          >  1:36906/1‑101 (MQ=255)                                                                      gACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcaca                         >  1:1147752/1‑101 (MQ=255)                                                                      gACATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcaca                         <  1:423616/101‑1 (MQ=255)                                                                        cATAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                       <  2:728895/101‑1 (MQ=255)                                                                         aTAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGTGCACTATGAAGGCCACAAGa                      >  1:1448522/1‑101 (MQ=255)                                                                         aTAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                      >  2:1529439/1‑101 (MQ=255)                                                                         aTAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                      >  2:1671780/1‑101 (MQ=255)                                                                          tAAATGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                     <  2:916997/101‑1 (MQ=255)                                                                           aaaTGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGGTc                    <  1:1430905/101‑1 (MQ=255)                                                                           aaaTGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc                    <  2:620452/101‑1 (MQ=255)                                                                             aTGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt                  >  2:639893/1‑101 (MQ=255)                                                                              tGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                 <  2:130408/101‑1 (MQ=255)                                                                              tGACCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                 <  2:1528446/101‑1 (MQ=255)                                                                                aCCGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga               >  2:1562389/1‑101 (MQ=255)                                                                                 ccGCTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag              <  1:886230/101‑1 (MQ=255)                                                                                    cTATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt           >  1:90375/1‑101 (MQ=255)                                                                                     tATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg          <  1:1410718/101‑1 (MQ=255)                                                                                     tATGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg          >  1:448128/1‑101 (MQ=255)                                                                                      aTGAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg         >  2:1251918/1‑101 (MQ=255)                                                                                        gAGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTg       <  1:1677930/101‑1 (MQ=255)                                                                                         aGCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGc      >  2:1208105/1‑101 (MQ=255)                                                                                          gCACTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGcc     <  2:995650/101‑1 (MQ=255)                                                                                             cTGCCATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt  <  2:44553/101‑1 (MQ=255)                                                                                                 |                                                                                                 AATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGATAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTT  >  REL606/238906‑239098

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