Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	363,518	A→G	Y930H (TAC→CAC)	lacZ ←	beta‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	363,518	0	A	G	100.0%	168.0 / NA	47	Y930H (TAC→CAC)	lacZ	beta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/22);  total (25/22)

ATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTACGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCTAATCCGAGCCAGTTTACCCGCTCTGCTACCTGCGCCAGCTGGCAGT  >  NC_000913/363381‑363649                                                                                                                                          |                                                                                                                                    atggcgatggcTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCg                                                                                                                                   >  2:947746/1‑140 (MQ=255)   ggcgatggcTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCggg                                                                                                                                 >  1:85948/1‑140 (MQ=255)   ggcgatggcTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCggg                                                                                                                                 <  1:635355/140‑1 (MQ=255)   ggcgatggcTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCggg                                                                                                                                 <  1:530789/140‑1 (MQ=255)       atggcTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGtat                                                                                                                             <  1:117290/140‑1 (MQ=255)        tggcTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGtata                                                                                                                            >  2:980436/1‑140 (MQ=255)        tggcTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGtata                                                                                                                            <  1:439192/140‑1 (MQ=255)                   cATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATATATGTCTGACa                                                                                                                 <  1:980436/140‑1 (MQ=255)                       aGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATgg                                                                                                             <  1:14789/140‑1 (MQ=255)                                tGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCa                                                                                                    <  1:947746/140‑1 (MQ=255)                                       aGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCa                                                                                             >  2:960321/1‑140 (MQ=255)                                         cGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCaaa                                                                                           >  1:134597/1‑140 (MQ=255)                                           gCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAAc                                                                                         >  1:529908/1‑140 (MQ=255)                                           gCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAAc                                                                                         >  2:903103/1‑140 (MQ=255)                                           gCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAAc                                                                                         >  2:16859/1‑140 (MQ=255)                                           gCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAAc                                                                                         >  1:951521/1‑140 (MQ=255)                                              gATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAgg                                                                                      <  2:536667/140‑1 (MQ=255)                                              gATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAgg                                                                                      <  2:134597/140‑1 (MQ=255)                                              gATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAgg                                                                                      <  1:960321/140‑1 (MQ=255)                                              gATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAgg                                                                                      <  1:903103/140‑1 (MQ=255)                                                tGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAggcg                                                                                    >  1:69674/1‑140 (MQ=255)                                                         ggAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAgg                                                                           >  1:190588/1‑140 (MQ=255)                                                                gccgcGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCg                                                                         >  1:861616/1‑135 (MQ=255)                                                                gccgcGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCg                                                                         <  2:861616/135‑1 (MQ=255)                                                                gccgcGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCgg                                                                    >  2:815190/1‑140 (MQ=255)                                                                gccgcGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCgg                                                                    >  2:858732/1‑140 (MQ=255)                                                                  cgcgCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGa                                                                  <  2:327060/140‑1 (MQ=255)                                                                                 ggCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGc                                                   >  2:688236/1‑140 (MQ=255)                                                                                 ggCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGc                                                   >  1:112046/1‑140 (MQ=255)                                                                                 ggCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGc                                                   >  1:564828/1‑140 (MQ=255)                                                                                 ggCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGc                                                   >  1:654531/1‑140 (MQ=255)                                                                                 ggCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTATGACAATGCCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGc                                                   >  2:597143/1‑140 (MQ=255)                                                                                   ccATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGccc                                                 <  2:654531/140‑1 (MQ=255)                                                                                    cATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCt                                                >  1:343525/1‑140 (MQ=255)                                                                                      tAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCTaa                                              >  2:177482/1‑140 (MQ=255)                                                                                        attcaattcGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCTAATc                                            <  1:673772/140‑1 (MQ=255)                                                                                           caattcGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCTAATCCGa                                         <  2:85948/140‑1 (MQ=255)                                                                                              ttcGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCTAATCCGAGcc                                      >  2:917783/1‑140 (MQ=255)                                                                                                  cgcgTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCTAATCCGAGCCAGtt                                  <  2:69674/140‑1 (MQ=255)                                                                                                   gcgTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCTAATCCGAGCCAGttt                                 >  2:592482/1‑140 (MQ=255)                                                                                                     gTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCTAATCCGAGCCAGTTTAc                               <  1:592482/140‑1 (MQ=255)                                                                                                     gTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCTAATCCGAGCCAGTTTAc                               <  2:9566/140‑1 (MQ=255)                                                                                                      tctcgccgcgcagACCGTTTTCGCGCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAACACTGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGGGGCCCTAATCCGAGCCAGTTTAcc                              <  2:225840/137‑1 (MQ=255)                                                                                                      tCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCTAATCCGAGCCAGTTTAcc                              >  2:337900/1‑140 (MQ=255)                                                                                                       cccGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCTAATCCGAGCCAGTTTAccc                             <  1:177482/140‑1 (MQ=255)                                                                                                                 agACCGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCTAATCCGAGCCAGTTTACCCGCTCTGCTAc                   >  1:873522/1‑140 (MQ=255)                                                                                                                    ccGTTTTCGCTCGGGAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCTAATCCGAGCCAGTTTACCCGCTCTGCTACCTg                <  1:16859/140‑1 (MQ=255)                                                                                                                                 ggAAGACGTGCGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCTAATCCGAGCC‑GTTTACCCGCTCTGCTACCTGCGCCAGCTGGCAGt  <  2:343525/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                    ATGGCGATGGCTGGTTTCCATCAGTTGCTGTTGACTGTAGCGGCTGATGTTGAACTGGAAGTCGCCGCGCCACTGGTGTGGGCCATAATTCAATTCGCGCGTCCCGCAGCGCAGACCGTTTTCGCTCGGGAAGACGTACGGGGTATACATGTCTGACAATGGCAGATCCCAGCGGTCAAAACAGGCGGCAGTAAGGCGGTCGGGATAGTTTTCTTGCGGCCCTAATCCGAGCCAGTTTACCCGCTCTGCTACCTGCGCCAGCTGGCAGT  >  NC_000913/363381‑363649

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