Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,332,930	A→G	E628G (GAG→GGG)	topA →	DNA topoisomerase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,332,930	0	A	G	100.0%	125.9 / NA	36	E628G (GAG→GGG)	topA	DNA topoisomerase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/16);  total (20/16)

GGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTA  >  NC_000913/1332803‑1333062                                                                                                                                |                                                                                                                                     ggtATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAac                                                                                                                          <  1:430970/140‑1 (MQ=255)         cccGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAac                                                                                                                          >  2:852198/1‑132 (MQ=255)         cccGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAac                                                                                                                          <  1:852198/132‑1 (MQ=255)         cccGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTa                                                                                                                  >  1:961032/1‑140 (MQ=255)            gAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAAcc                                                                                                               <  2:149249/140‑1 (MQ=255)             aaCCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCt                                                                                                              >  2:976515/1‑140 (MQ=255)                      gTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGa                                                                                                     >  1:352848/1‑140 (MQ=255)                                 cATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGc                                                                                          >  2:739571/1‑140 (MQ=255)                                 cATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGc                                                                                          >  2:354099/1‑140 (MQ=255)                                 cATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACC‑GGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCt                                                                                         >  2:24126/1‑140 (MQ=255)                                  aTTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCt                                                                                         <  2:762013/140‑1 (MQ=255)                                   ttGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACC‑GGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGa                                                                                       <  1:24126/140‑1 (MQ=255)                                    tGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGa                                                                                       <  1:354099/140‑1 (MQ=255)                                    tGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGa                                                                                       <  2:961032/140‑1 (MQ=255)                                    tGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGa                                                                                       <  1:976515/140‑1 (MQ=255)                                      aCTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGa                                                                                       >  2:341509/1‑138 (MQ=255)                                      aCTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGa                                                                                       <  1:341509/138‑1 (MQ=255)                                           ccGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCt                                                                                >  1:906696/1‑140 (MQ=255)                                            cGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTg                                                                               <  2:906696/140‑1 (MQ=255)                                              aCTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGa                                                                             >  2:847451/1‑140 (MQ=255)                                                                        cacaGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGACGCTGAAACCAACgcgc                                                   >  2:412091/1‑140 (MQ=255)                                                                                         tATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAAcgtcgt                                  <  1:787980/140‑1 (MQ=255)                                                                                             ccTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGACGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGcc                              <  1:412091/140‑1 (MQ=255)                                                                                               ttGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCg                            <  1:42505/140‑1 (MQ=255)                                                                                                 ggCTGTTCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGaa                          <  1:847451/140‑1 (MQ=255)                                                                                                      ttCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCg                     >  2:395907/1‑140 (MQ=255)                                                                                                      ttCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCg                     >  2:895709/1‑140 (MQ=255)                                                                                                      ttCTGGCTATGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCg                     >  1:73339/1‑140 (MQ=255)                                                                                                              aTGCGCTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCg             >  2:779043/1‑140 (MQ=255)                                                                                                                gcgcTGCCGCCGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGAt           >  2:855732/1‑140 (MQ=255)                                                                                                                     gccgccGAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACa      >  2:70212/1‑140 (MQ=255)                                                                                                                      ccgccgAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATg          <  1:690130/135‑1 (MQ=255)                                                                                                                      ccgccgAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATg          >  2:690130/1‑135 (MQ=255)                                                                                                                      ccgccgAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGa        >  2:938180/1‑137 (MQ=255)                                                                                                                      ccgccgAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGa        <  1:938180/137‑1 (MQ=255)                                                                                                                         ccgAAAGGGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTa  >  2:40981/1‑140 (MQ=255)                                                                                                                                |                                                                                                                                     GGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGCTGAACGTGCTGGAAGGCGAAGATGCTGAAACCAACGCGCTGCGCGCAAAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTA  >  NC_000913/1332803‑1333062

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