Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	363,997	C→T	W770* (TGG→TAG)	lacZ ←	beta‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	363,997	0	C	T	100.0%	105.7 / NA	31	W770* (TGG→TAG)	lacZ	beta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/13);  total (18/13)

CCCGCCGCCTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCA  >  NC_000913/363862‑364120                                                                                                                                        |                                                                                                                            cccGCCGCCTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACAt                                                                                                                         <  1:29531/140‑1 (MQ=255)    gccgccTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCtg                                                                                                                      >  2:585829/1‑140 (MQ=255)    gccgccTTCCAGAGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCtg                                                                                                                      >  1:96208/1‑140 (MQ=255)          ttCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTgaaag                                                                                                                <  2:96208/140‑1 (MQ=255)                   tCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGa                                                                                                       >  2:726924/1‑140 (MQ=255)                    cGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGAc                                                                                                      <  1:640524/140‑1 (MQ=255)                       cccAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTgg                                                                                                   >  1:606419/1‑140 (MQ=255)                                 aGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGAc                                                                                                      <  1:684560/127‑1 (MQ=255)                                 aGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGAc                                                                                                      >  2:684560/1‑127 (MQ=255)                                 aGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAAtt                                                                                         >  1:618454/1‑140 (MQ=255)                                 aGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAAtt                                                                                         >  1:552797/1‑140 (MQ=255)                                      cAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCaa                                                                                    >  1:661896/1‑140 (MQ=255)                                      cAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCaa                                                                                    <  2:618454/140‑1 (MQ=255)                                       aaTGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAAtttcaaaa                                                                                   >  1:809138/1‑134 (MQ=255)                                       aaTGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCACAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAAc                                                                                   <  2:243584/140‑1 (MQ=255)                                            gggTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCtta                                                                              <  1:726924/140‑1 (MQ=255)                                              gTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCttatt                                                                            >  1:725333/1‑140 (MQ=255)                                                                 tGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAAATGCCAACGCTTATTACCCAGCTCGATGCaaaaa                                                         >  1:49730/1‑140 (MQ=255)                                                                     gTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCa                                                     <  1:265545/140‑1 (MQ=255)                                                                      ttATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCAt                                                    <  1:377885/140‑1 (MQ=255)                                                                             gCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCtg                                             <  2:552797/140‑1 (MQ=255)                                                                             gCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCtg                                             <  2:197451/140‑1 (MQ=255)                                                                              cGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAAATTCCAACGCTTATTATCCCGCTCGATGCAAAAAACCATTTTGCtgg                                            >  2:635884/1‑140 (MQ=255)                                                                              cGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTCTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTTCCAACGCTTATTACCCCGCTCCATGCAAAAATCCATTTCGCtgg                                            >  2:673492/1‑140 (MQ=255)                                                                                            aCTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCggg                              <  1:635884/140‑1 (MQ=255)                                                                                               gATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATg                           <  2:606419/140‑1 (MQ=255)                                                                                                aTCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTTAATTGCCAACGCTTATTACCCCGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATgg                          >  2:71239/1‑140 (MQ=255)                                                                                                aTCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTTCCAACCCTTATTACCCCCCCCCATGCAAAAAACCATTTTTCTGGTGGGCAGATGCggggggg                          >  2:188063/1‑136 (MQ=255)                                                                                                  cgcgCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTTCCCACGCTTTTTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCg                        >  2:686078/1‑140 (MQ=255)                                                                                                  cgcgCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAAATTCCAACGCTTATTACCCAGCCCCGTGCAAAAAAACATTTCGCTGGGGGTCAGATGCggggggggg                        >  2:207478/1‑134 (MQ=255)                                                                                                                        ttttttATCGCCAATCTACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAACTCGATGCAAAAAACCATTTCGCTGGTGGGCAGATGCGGGATGGCGGGGGACGCGGCggggggggtgc  >  2:51510/1‑133 (MQ=255)                                                                                                                                        |                                                                                                                            CCCGCCGCCTTCCAGCGTTCGACCCAGGCGTTAGGGTCAATGCGGGTCGCTTCACTTACGCCAATGTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCGTCA  >  NC_000913/363862‑364120

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