Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,332,839	G→A	D598N (GAC→AAC)	topA →	DNA topoisomerase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,332,839	0	G	A	100.0%	115.9 / NA	34	D598N (GAC→AAC)	topA	DNA topoisomerase 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/15);  total (19/15)

GGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGC  >  NC_000913/1332706‑1332974                                                                                                                                      |                                                                                                                                        ggCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGCAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGcc                                                                                                                                   <  1:794418/140‑1 (MQ=255)   cAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCg                                                                                                                                 >  2:98535/1‑140 (MQ=255)                       aaaGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATggg                                                                                                             >  2:783741/1‑140 (MQ=255)                       aaaGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATggg                                                                                                             <  1:555169/140‑1 (MQ=255)                       aaaGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGCTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACT‑CCCGACTTGTGGTCGCAAAATgggg                                                                                                            >  1:166495/1‑140 (MQ=255)                                     tCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGa                                                                                               >  1:599594/1‑140 (MQ=255)                                        cttcttCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAg                                                                                              <  2:136213/138‑1 (MQ=255)                                        cttcttCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAg                                                                                              >  1:136213/1‑138 (MQ=255)                                        cttcttCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCa                                                                                            >  2:708599/1‑140 (MQ=255)                                        cttcttCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCa                                                                                            >  2:373706/1‑140 (MQ=255)                                        cttcttCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCa                                                                                            <  1:98535/140‑1 (MQ=255)                                             tctcGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCggg                                                                                       >  1:122099/1‑140 (MQ=255)                                               tcGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGt                                                                                     <  2:599594/140‑1 (MQ=255)                                                 ggATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCgggg                                                                                      <  1:468658/137‑1 (MQ=255)                                                 ggATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCgggg                                                                                      >  2:468658/1‑137 (MQ=255)                                                             gcagTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGtt                                                                       <  1:373706/140‑1 (MQ=255)                                                                         aGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATgcgc                                                           >  2:5733/1‑140 (MQ=255)                                                                         aGCTGAAAAAGATCCGGAAGAGGGAGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATgcgc                                                           >  1:84908/1‑140 (MQ=255)                                                                               aaaaGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTgccgc                                                     >  2:75243/1‑140 (MQ=255)                                                                                      ccGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACT‑CCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAagag                                             <  2:166495/140‑1 (MQ=255)                                                                                        tccggagGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGc                                            <  1:260415/136‑1 (MQ=255)                                                                                              gggTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCa                                      <  2:84908/140‑1 (MQ=255)                                                                                               ggtggtATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCaa                                     >  2:44470/1‑140 (MQ=255)                                                                                                  ggtATGCGCCCGAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAac                                  <  2:122099/140‑1 (MQ=255)                                                                                                             gAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAAcc                       <  1:75243/140‑1 (MQ=255)                                                                                                             gAACCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATGAAcc                       >  2:739173/1‑140 (MQ=255)                                                                                                              aaCCAGATGGTTCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCt                      >  2:63599/1‑140 (MQ=255)                                                                                                                         tCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAAcc                       <  1:9642/128‑1 (MQ=255)                                                                                                                         tCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAAcc                       >  2:9642/1‑128 (MQ=255)                                                                                                                         tCTGACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGaaa           >  2:593272/1‑140 (MQ=255)                                                                                                                            gACCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACg        <  1:5733/140‑1 (MQ=255)                                                                                                                             aCCAGCATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGa       <  1:708599/140‑1 (MQ=255)                                                                                                                                  cATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGc  >  2:441684/1‑140 (MQ=255)                                                                                                                                  cATTAACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGc  >  1:418944/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                        GGCAAATCACGAAGCAGAGTGGAAAGCTGTACTGGATCACTTCTTCTCGGATTTCACCCAGCAGTTAGATAAAGCTGAAAAAGATCCGGAAGAGGGTGGTATGCGCCCGAACCAGATGGTTCTGACCAGCATTGACTGCCCGACTTGTGGTCGCAAAATGGGGATTCGCACAGCGAGCACCGGGGTATTCCTTGGCTGTTCTGGCTATGCGCTGCCGCCGAAAGAGCGTTGCAAAACCACCATTAACCTGGTGCCGGAAAACGAAGTGC  >  NC_000913/1332706‑1332974

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