Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	364,333	G→A	A658V (GCG→GTG)	lacZ ←	beta‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	364,333	0	G	A	100.0%	130.9 / NA	39	A658V (GCG→GTG)	lacZ	beta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/20);  total (19/20)

TTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCG  >  NC_000913/364197‑364456                                                                                                                                         |                                                                                                                            ttGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCAcca                                                                                                                          <  1:627524/140‑1 (MQ=255)      cTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGCGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATc                                                                                                                     <  1:423280/140‑1 (MQ=255)         cgcgTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAg                                                                                                                  >  2:517849/1‑140 (MQ=255)         cgcgTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAg                                                                                                                  >  2:664364/1‑140 (MQ=255)           cgTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTg                                                                                                                >  1:668835/1‑140 (MQ=255)                          agTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAg                                                                                                          <  2:650540/131‑1 (MQ=255)                          agTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAg                                                                                                          >  1:650540/1‑131 (MQ=255)                          agTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATc                                                                                                 >  1:146654/1‑140 (MQ=255)                               cccGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTAt                                                                                            >  2:62229/1‑140 (MQ=255)                                 gacagcTCTCCGGCTGCTGTAGCTCAGGCAGTTCAATCAACTGTTTACTTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGa                                                                                          <  2:334173/136‑1 (MQ=255)                                 cGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGa                                                                                          <  1:590268/140‑1 (MQ=255)                                      ctctCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGaa                                                                                     <  1:62229/140‑1 (MQ=255)                                      ctctCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGaa                                                                                     <  1:517849/140‑1 (MQ=255)                                      ctctCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGaa                                                                                     <  1:664364/140‑1 (MQ=255)                                      ctctCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGaa                                                                                     <  2:685714/140‑1 (MQ=255)                                      ctctCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGaa                                                                                     <  2:668835/140‑1 (MQ=255)                                         tcCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAg                                                                                  >  1:703215/1‑140 (MQ=255)                                         tcCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAg                                                                                  >  2:109501/1‑140 (MQ=255)                                                 cGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGAAGGAGCTAGTTATCACTATGAGGGAACAGGTactgtc                                                                          >  1:334173/1‑135 (MQ=255)                                                     aGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGt                                                                      <  1:569215/140‑1 (MQ=255)                                                           ggCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTg                                                                        <  1:614532/132‑1 (MQ=255)                                                           ggCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTg                                                                        >  2:614532/1‑132 (MQ=255)                                                                ttCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATgg                                                           >  1:586197/1‑140 (MQ=255)                                                                ttCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACAGAATAGATATTCGCTGGTCACTTCGATgg                                                           >  1:373624/1‑140 (MQ=255)                                                                     tcaaCTGTTTACCTTGTGGAGCGACATCCAGAGGCACATCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGc                                                      <  1:307202/140‑1 (MQ=255)                                                                          tGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGcc                                                     >  2:618920/1‑136 (MQ=255)                                                                          tGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGcc                                                     <  1:618920/136‑1 (MQ=255)                                                                          tGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGa                                                 >  2:709963/1‑140 (MQ=255)                                                                             ttACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGACCAGTTATTCGCTGGTCACTTCCATGGTTTGCCCGGATaa                                              >  2:307202/1‑140 (MQ=255)                                                                              tACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATaaa                                             <  1:109501/140‑1 (MQ=255)                                                                                ccTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACg                                           <  1:709963/140‑1 (MQ=255)                                                                                      ggAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTGGCTTGTCACTTCGATGATTTGCCCGGATAAACGGAACTg                                     >  2:227708/1‑140 (MQ=255)                                                                                      ggAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTg                                     >  2:47120/1‑140 (MQ=255)                                                                                           gACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGaaaa                                <  2:703215/140‑1 (MQ=255)                                                                                             cATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAAc                              <  2:146654/140‑1 (MQ=255)                                                                                             cATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAAc                              <  1:47120/140‑1 (MQ=255)                                                                                               tCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAActg                            >  2:749313/1‑140 (MQ=255)                                                                                                      gCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTgg                     >  1:815445/1‑140 (MQ=255)                                                                                                                         cGGCTTACCATCCAGCACCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAgcg  <  2:586197/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                            TTGCACTACGCGTACTGTGAGCCAGAGTTGCCCGGCGCTCTCCGGCTGCGGTAGTTCAGGCAGTTCAATCAACTGTTTACCTTGTGGAGCGACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTACCATCCAGCGCCACCATCCAGTGCAGGAGCTCGTTATCGCTATGACGGAACAGGTATTCGCTGGTCACTTCGATGGTTTGCCCGGATAAACGGAACTGGAAAAACTGCTGCTGGTGTTTTGCTTCCGTCAGCG  >  NC_000913/364197‑364456

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