Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	664,688	+T	intergenic (‑79/‑47)	gltL ← / → insJ‑2	ATP‑binding protein of glutamate/aspartate transport system; b0652_1/IS150 hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	664,684	1	.	T	100.0%	182.0 / NA	47	intergenic (‑75/‑51)	gltL/insJ‑2	ATP‑binding protein of glutamate/aspartate transport system; b0652_1/IS150 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (26/21);  total (26/21)

CACGCTGCTGCTGACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAG‑TTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACC  >  REL606/664609‑664753                                                                             |                                                                      cacGCTGCTGCTGACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCcagcag                                                                        <  1:1691147/76‑1 (MQ=255) cacGCTGCTGCTGACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCcagcag                                                                        >  1:1332290/1‑76 (MQ=255)     ctgctgctgACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAG‑tttt                                                                   >  1:1417045/1‑76 (MQ=255)        ctgctgACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAG‑tttttgt                                                                <  2:980906/76‑4 (MQ=255)        ctgctgACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAG‑tttttgt                                                                <  2:4066/76‑4 (MQ=255)             gACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGc                                                            >  1:870785/1‑76 (MQ=255)              gccgccgGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCa                                                           <  2:1777450/75‑1 (MQ=255)               ccgccgGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCAc                                                          >  2:1312380/1‑76 (MQ=255)                 gccaGAAAGCTGAGGGGGAAACTGATTGGCGTGAGCAGAAAGCCCGACACGGTCCAGCAGTTTTTGTACTGGAccc                                                        <  2:305211/76‑1 (MQ=255)                     gAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCAttt                                                    >  2:169082/1‑76 (MQ=255)                     gAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCAttt                                                    <  2:10057/76‑1 (MQ=255)                      aaaGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCAtttt                                                   >  2:512755/1‑76 (MQ=255)                      aaaGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCAtttt                                                   >  2:322338/1‑76 (MQ=255)                      aaaGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCAtttt                                                   <  1:976271/76‑1 (MQ=255)                       aaGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTttg                                                  >  2:243830/1‑76 (MQ=255)                            gcgcCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGa                                             >  2:8258/1‑76 (MQ=255)                            gcgcCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGa                                             >  2:830669/1‑76 (MQ=255)                            gcgcCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGa                                             <  1:292197/76‑1 (MQ=255)                                attaaaCTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGAt                                         >  2:110059/4‑76 (MQ=38)                                   aaaCTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGaa                                      >  2:398680/1‑76 (MQ=255)                                   aaaCTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGaa                                      >  2:398664/1‑76 (MQ=255)                                      cTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATg                                   <  1:1314794/76‑1 (MQ=255)                                        tATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGa                                 <  2:827870/76‑1 (MQ=255)                                        tATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGa                                 >  2:1395201/1‑76 (MQ=255)                                          ttGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt                               >  1:1362039/1‑76 (MQ=255)                                          ttGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAAt                               <  1:1230611/76‑1 (MQ=255)                                               gTGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGccc                          <  2:1373994/76‑1 (MQ=255)                                                tGAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGcccc                         >  2:788333/1‑76 (MQ=255)                                                 gAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt                        >  1:1186266/1‑76 (MQ=255)                                                 gAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt                        <  1:759372/76‑1 (MQ=255)                                                 gAGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCt                        <  2:1065723/76‑1 (MQ=255)                                                  aGCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTa                       >  2:1038019/1‑76 (MQ=255)                                                   gCAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTaa                      <  2:1851861/76‑1 (MQ=255)                                                    cAGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAat                     <  1:934812/76‑1 (MQ=255)                                                     aGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAata                    >  2:1116088/1‑76 (MQ=255)                                                     aGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAata                    >  1:246826/1‑76 (MQ=255)                                                     aGAAAGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAata                    >  1:76307/1‑76 (MQ=255)                                                       aaaGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAatatt                  >  2:15502/1‑75 (MQ=255)                                                       aaaGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATg                  <  1:1069681/76‑1 (MQ=255)                                                       aaaGCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATg                  <  1:1752546/76‑1 (MQ=255)                                                          gCCCGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCa               <  2:100406/76‑1 (MQ=255)                                                           cccGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCaa              >  2:1861713/1‑76 (MQ=255)                                                             cGACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAg            <  1:1530946/76‑1 (MQ=255)                                                              gACACGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGc           >  2:823844/1‑76 (MQ=255)                                                                cacGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCa         >  2:1848718/1‑76 (MQ=255)                                                                cacGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCa         <  1:1492230/76‑1 (MQ=255)                                                                cacGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCa         >  2:1848704/1‑76 (MQ=255)                                                                  cGCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaaa       >  1:312876/1‑76 (MQ=255)                                                                   gCTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCaaaa      <  2:452485/76‑1 (MQ=255)                                                                    cTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAAt     >  2:1442758/1‑76 (MQ=255)                                                                    cTCCAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAAt     <  1:8258/76‑1 (MQ=255)                                                                      ccAGCAGTTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAc   <  1:66985/76‑1 (MQ=255)                                                                       cagcagtTTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATAcc  >  2:1093517/1‑76 (MQ=255)                                                                             |                                                                      CACGCTGCTGCTGACCGCCGGAAAGCTGCGCCGGAAACTTATTGGCGTGAGCAGAAAGCCCGACACGCTCCAGCAG‑TTTTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACC  >  REL606/664609‑664753

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