Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	666,133	+T	pseudogene (355/387 nt)	gltL ←	ATP‑binding protein of glutamate/aspartate transport system; b0652_1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	666,130	1	.	T	100.0%	254.0 / NA	67	pseudogene (358/387 nt)	gltL	ATP‑binding protein of glutamate/aspartate transport system; b0652_1
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (40/27);  total (40/27)

TGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA‑TTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACA  >  REL606/666043‑666225                                                                                         |                                                                                                tGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttt                                                                                     >  2:775918/1‑101 (MQ=255) tGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttt                                                                                     >  2:1150422/1‑101 (MQ=255) tGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttt                                                                                     <  1:872176/101‑1 (MQ=255)  gAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCtttt                                                                                    <  2:1321153/101‑1 (MQ=255)   aaaTTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                                                                   <  2:54763/101‑1 (MQ=255)   aaaTTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                                                                   <  1:1486536/101‑1 (MQ=255)   aaaTTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCttttt                                                                                   <  2:45971/101‑1 (MQ=255)     aTTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCa                                                                                 >  1:534251/1‑101 (MQ=255)         aaaGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                                                                             <  1:1021004/101‑1 (MQ=255)         aaaGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAg                                                                             >  2:473770/1‑101 (MQ=255)           aGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGcc                                                                           >  1:386256/1‑101 (MQ=255)           aGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGcc                                                                           <  1:514019/101‑1 (MQ=255)                 gACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCt                                                                     >  2:754383/1‑101 (MQ=255)                      cAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc                                                                >  2:495947/1‑101 (MQ=255)                      cAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATc                                                                >  2:798676/1‑101 (MQ=255)                       aaTTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg                                                               >  1:1227492/1‑101 (MQ=255)                       aaTTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcg                                                               >  1:1235210/1‑101 (MQ=255)                        aTTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATcgc                                                              >  1:206869/1‑101 (MQ=255)                          tGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGt                                                            >  1:798510/1‑101 (MQ=255)                           gAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGtt                                                           >  1:1229718/1‑101 (MQ=255)                            aaTATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGttt                                                          <  1:150360/101‑1 (MQ=255)                              tatCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTaa                                                        >  1:1126712/1‑101 (MQ=255)                                  ggAATCAGTCCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCAc                                                    <  1:495952/101‑1 (MQ=255)                                   gAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACt                                                   <  1:183604/101‑1 (MQ=255)                                    aaTCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACtt                                                  >  1:1560540/1‑101 (MQ=255)                                     aTCAGACCTATATGGCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACttt                                                 <  2:1143716/101‑1 (MQ=255)                                     aTCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACttt                                                 >  2:68144/1‑101 (MQ=255)                                     aTCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACttt                                                 >  2:532109/1‑101 (MQ=255)                                        aGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCAc                                              >  1:1219144/1‑101 (MQ=255)                                         gACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCAcc                                             <  1:1246411/101‑1 (MQ=255)                                          aCCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCt                                            >  2:1273793/1‑101 (MQ=255)                                           ccTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTg                                           >  1:164424/1‑101 (MQ=255)                                            cTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTgc                                          >  2:532788/1‑101 (MQ=255)                                               tatGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGcc                                       <  1:775923/101‑1 (MQ=255)                                                atGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCa                                      >  2:367352/1‑101 (MQ=255)                                                 tGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAg                                     <  2:248679/101‑1 (MQ=255)                                                   ccTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAggg                                   >  1:1658054/1‑101 (MQ=255)                                                    cTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGt                                  >  2:910633/1‑101 (MQ=255)                                                     tCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGATc                                 <  2:1071544/101‑1 (MQ=255)                                                       gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAg                               >  2:944486/1‑101 (MQ=255)                                                       gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAg                               >  1:595948/1‑101 (MQ=255)                                                        tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAgg                              >  2:1210726/1‑101 (MQ=255)                                                        tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAgg                              >  1:106661/1‑101 (MQ=255)                                                        tgtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAgg                              <  2:265284/101‑1 (MQ=255)                                                         gtTTAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGt                             >  1:146578/1‑101 (MQ=255)                                                          tttAACTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGTCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGtt                            <  2:206923/101‑1 (MQ=255)                                                             aaCTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTc                         >  1:736755/1‑101 (MQ=255)                                                              aCTGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCg                        <  1:798681/101‑1 (MQ=255)                                                                tGTCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGAt                      >  2:523231/1‑101 (MQ=255)                                                                tGTCCAACTTTTTGGGGGCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGAt                      >  2:1428502/1‑101 (MQ=255)                                                                  tCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATaa                    <  2:1210062/101‑1 (MQ=255)                                                                  tCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATaa                    <  2:456252/101‑1 (MQ=255)                                                                  tCCAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATaa                    >  2:1503475/1‑101 (MQ=255)                                                                   ccAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAAt                   <  2:274111/101‑1 (MQ=255)                                                                   ccAACTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAAt                   >  1:1257312/1‑101 (MQ=255)                                                                     aaCTTTTTGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCg                 <  1:331095/101‑1 (MQ=255)                                                                         ttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAg             <  1:532793/101‑1 (MQ=255)                                                                         ttttGGGGTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAg             <  2:1265962/101‑1 (MQ=255)                                                                              gggTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAg        <  1:532114/101‑1 (MQ=255)                                                                               ggTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAgg       >  1:418779/1‑101 (MQ=255)                                                                                gTCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAggg      >  1:1144591/1‑101 (MQ=255)                                                                                 tCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGa     >  1:1543606/1‑101 (MQ=255)                                                                                 tCAGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGa     >  2:898719/1‑101 (MQ=255)                                                                                   aGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAAc   <  2:1126707/101‑1 (MQ=255)                                                                                   aGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAAc   >  1:347107/1‑101 (MQ=255)                                                                                   aGTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAAc   >  2:1637490/1‑101 (MQ=255)                                                                                    gTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACa  <  1:569397/101‑1 (MQ=255)                                                                                    gTACATTTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACa  <  1:1579741/101‑1 (MQ=255)                                                                                         |                                                                                                TGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACA‑TTTCAGGGCTTTTTCACGAGCCGGCGCTTTATCGCGTTTAAGCACTTTCACCTGCGCCAGGGTCAGGTTTTCGATAATCGACAGATGAGGGAACA  >  REL606/666043‑666225

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: ‑