Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,333,815	T→G	pseudogene (174/189 nt)	yciX →	predicted protein; b4522(b4574)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	1,333,815	0	T	G	100.0%	120.9 / NA	31	pseudogene (174/189 nt)	yciX	predicted protein; b4522(b4574)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (17/14);  total (17/14)

CGTTATTCCAGACGACTGGCAACTAACATCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGATCAGCCGAAACAATAATTATCATCATTCTTATTACCCATTTTTAATGAATTAAAGGGCTTTTAATACA  >  REL606/1333746‑1333882                                                                      |                                                                    cgtTATTCCAGACGACTGGCAACTAACATCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCg                                                               >  1:42261/1‑76 (MQ=255)   ttATTCCAGACGACTGGCAACTAACATCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGaa                                                             >  2:1092850/1‑76 (MQ=255)   ttATTCCAGACGACTGGCAACTAACATCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGaa                                                             >  2:222434/1‑76 (MQ=255)     aTTCCAGACGACTGGCAACTAACATCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAAc                                                           <  1:669156/76‑1 (MQ=255)        ccAGACGACTGGCAACTAACATCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACaat                                                        >  1:1238535/1‑76 (MQ=255)         cAGACGACTGGCAACTAACATCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACaata                                                       <  1:760807/76‑1 (MQ=255)         cAGACGACTGGCAACTAACATCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACaata                                                       <  2:498211/76‑1 (MQ=255)             cgacTGGCAACTAACATCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTa                                                   >  1:1248638/1‑76 (MQ=255)                 tGGCAACTAACATCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTatcat                                               <  1:506827/76‑1 (MQ=255)                 tGGCAACTAACATCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTatcat                                               >  2:1340564/1‑76 (MQ=255)                  ggCAACTAACATCGCAGCAGCAAGCCTTTATAGAACGGTTTGCTGAAGATGAGCAGCCGAAACAATAATTatcatc                                              <  1:925885/76‑1 (MQ=255)                       cTAACATCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCt                                         >  2:660196/1‑76 (MQ=255)                         aaCATCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCtta                                       >  1:207823/1‑76 (MQ=255)                         aaCATCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCtta                                       >  1:112630/1‑76 (MQ=255)                            aTCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCttatta                                    >  1:679773/1‑76 (MQ=255)                                cagcagcaAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCa                                <  2:417652/76‑1 (MQ=255)                                cagcagcaAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCa                                >  1:971845/1‑76 (MQ=255)                                 agcagcaAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCAt                               <  1:204092/76‑1 (MQ=255)                                   cagcaAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCAttt                             >  1:1088487/1‑76 (MQ=255)                                   cagcaAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCAttt                             <  1:138067/76‑1 (MQ=255)                                   cagcaAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCAttt                             <  2:1356373/76‑1 (MQ=255)                                   cagcaAGCCTTTATAGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCAttt                             <  2:1052437/76‑1 (MQ=255)                                              tataGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCATTTTTAATGAATTa                  <  1:918244/76‑1 (MQ=255)                                               ataGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCATTTTTAATGAATTaa                 >  2:1350620/1‑76 (MQ=255)                                                taGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCATTTTTAATGAATTaaa                >  2:410580/1‑76 (MQ=255)                                                 aGAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCATTTTTAATGAATTAAAg               <  1:991892/76‑1 (MQ=255)                                                  gAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCATTTTTAATGAATTAAAgg              <  1:968450/76‑1 (MQ=255)                                                  gAACTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCATTTTTAATGAATTAAAgg              <  1:1218810/76‑1 (MQ=255)                                                    aCTGTTTGCTGAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCATTTTTAATGAATTAAAGGGc            >  1:474326/1‑76 (MQ=255)                                                              gAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCATTTTTAATGAATTAAAGGGCTTTTAATaca  >  1:492014/1‑76 (MQ=255)                                                              gAAGATGAGCAGCCGAAACAATAATTATCATCATTCTTATTACCCATTTTTAATGAATTAAAGGGCTTTTAATaca  >  2:896408/1‑76 (MQ=255)                                                                      |                                                                    CGTTATTCCAGACGACTGGCAACTAACATCGCAGCAGCAAGCCTTTATAGAACTGTTTGCTGAAGATGATCAGCCGAAACAATAATTATCATCATTCTTATTACCCATTTTTAATGAATTAAAGGGCTTTTAATACA  >  REL606/1333746‑1333882

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