Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,333,620	G→A	G858D (GGC→GAC)	topA →	DNA topoisomerase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,333,620	0	G	A	100.0%	164.7 / NA	48	G858D (GGC→GAC)	topA	DNA topoisomerase 1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/26);  total (22/26)

TTATCTGGCCGATGCGCCACAGCAGGATCCGGAAGGTAATAAGACCATGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGGCAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAG  >  NC_000913/1333483‑1333753                                                                                                                                          |                                                                                                                                      ttATCTGGCCGATGCGCCACAGCAGGATCCGGAAGGTAATAAGACCATGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACa                                                                                                                                     >  1:776473/1‑140 (MQ=255)   aTCTGGCCGATGCGCCACAGCAGGATCCGGAAGGTAATAAGACCATGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACaaa                                                                                                                                   <  1:292070/140‑1 (MQ=255)         ccGATGCGCCACAGCAGGATCCGGAAGGTAATAAGACCATGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGtt                                                                                                                             >  1:829518/1‑140 (MQ=255)           gATGCGCCACAGCAGGATCCGGAAGGTAATAAGACCATGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGa                                                                                                                           <  1:542368/140‑1 (MQ=255)            aTGCGCCACAGCAGGATCCGGAAGGTAATAAGACCATGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACaaa                                                                                                                                   <  2:358568/131‑1 (MQ=255)            aTGCGCCACAGCAGGATCCGGAAGGTAATAAGACCATGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACaaa                                                                                                                                   >  1:358568/1‑131 (MQ=255)                gcCACAGCAGGATCCGGAAGGTAATAAGACCATGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGaa                                                                                                                      <  2:829518/140‑1 (MQ=255)                gcCACAGCAGGATCCGGAAGGTAATAAGACCATGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGaa                                                                                                                      <  2:70773/140‑1 (MQ=255)                                gAAGGTAATAAGACCATGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAAt                                                                                                      <  2:908439/140‑1 (MQ=255)                                 aaGGTAATAAGACCATGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAAtt                                                                                                     >  2:126381/1‑140 (MQ=255)                                        taaGACCATGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAg                                                                                              >  2:160236/1‑140 (MQ=255)                                              cATGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGttttt                                                                                        >  1:300046/1‑140 (MQ=255)                                              cATGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGttttt                                                                                        >  2:112153/1‑140 (MQ=255)                                               aTGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAAtt                                                                                                     >  2:718654/1‑126 (MQ=255)                                               aTGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAAtt                                                                                                     <  1:718654/126‑1 (MQ=255)                                                tGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTtat                                                                                      <  1:112153/140‑1 (MQ=255)                                                   ttCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATaaa                                                                                   >  2:560918/1‑140 (MQ=255)                                                    tCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAAc                                                                                  <  1:126381/140‑1 (MQ=255)                                                     cGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGt                                                                                            >  1:878223/1‑129 (MQ=255)                                                     cGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGt                                                                                            <  2:878223/129‑1 (MQ=255)                                                           aGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACa                                                                                 <  2:516724/134‑1 (MQ=255)                                                           aGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACa                                                                                 >  1:516724/1‑134 (MQ=255)                                                           aGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGt                                                                           >  1:901592/1‑140 (MQ=255)                                                               gTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTcgcg                                                                       <  2:776473/140‑1 (MQ=255)                                                                        aaCAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGcc                                                              <  2:714982/140‑1 (MQ=255)                                                                           aGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCtt                                                           >  1:858566/1‑140 (MQ=255)                                                                            gCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCttt                                                          >  1:816570/1‑140 (MQ=255)                                                                                      tctTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCat                                                <  1:137836/140‑1 (MQ=255)                                                                                          cGGAAAAAGACGGAAAGGCGCCTGCCTGGCCAGCATTTTATGTTGATGAAAAATGGGTGGAAGGAAAAAAAAAACCTTTATTTCTGTAAGGTTTTTATAAACAAAGGGCCGCAAAACGGGCCCTTTTTTTTTGCAtattt                                            >  2:235605/1‑139 (MQ=255)                                                                                            gAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATattatt                                          >  2:636299/1‑140 (MQ=255)                                                                                              aaaaGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTAtttt                                        <  2:858566/140‑1 (MQ=255)                                                                                                     ggAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAGCCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTca                                 <  2:17560/140‑1 (MQ=255)                                                                                                     ggAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTca                                 <  1:636299/140‑1 (MQ=255)                                                                                                     ggAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTca                                 >  2:535394/1‑140 (MQ=255)                                                                                                             gACTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTCACACCtata                         >  2:259626/1‑140 (MQ=255)                                                                                                              aCTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTca                                 >  1:616939/1‑131 (MQ=255)                                                                                                              aCTGGCTGGTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTca                                 <  2:616939/131‑1 (MQ=255)                                                                                                                 ggctggTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTCACACCTATACACt                     <  1:259626/140‑1 (MQ=255)                                                                                                                 ggctggTCAGCATTTTATGTTGATGACAAATGGGTTGAAAGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTCACACCTATACACt                     <  1:714620/140‑1 (MQ=255)                                                                                                                     ggTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAgg                 >  2:565214/1‑140 (MQ=255)                                                                                                                     ggTCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAgg                 >  1:897901/1‑140 (MQ=255)                                                                                                                       tCAGCATTTTATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCt               <  1:535394/140‑1 (MQ=255)                                                                                                                                tATGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGat      >  1:790331/1‑140 (MQ=255)                                                                                                                                  tGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGatat    <  2:300046/140‑1 (MQ=255)                                                                                                                                  tGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGatat    <  2:790331/140‑1 (MQ=255)                                                                                                                                  tGTTGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGatat    <  1:560918/140‑1 (MQ=255)                                                                                                                                    ttGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAg  <  2:685842/140‑1 (MQ=255)                                                                                                                                    ttGATGACAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAg  <  2:901592/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                      TTATCTGGCCGATGCGCCACAGCAGGATCCGGAAGGTAATAAGACCATGGTTCGCTTTAGCCGTAAAACCAAACAGCAATATGTCTCTTCGGAAAAAGACGGAAAGGCGACTGGCTGGTCAGCATTTTATGTTGATGGCAAATGGGTTGAAGGAAAAAAATAACCTTTAATTCTGTCAGGTTTTTATAAACAAAGGGTCGCGAAAGCGGCCCTTTTTTATTGCATATTATTTTTTCTTCACACCTATACACTAAGGCTATAAATGATATAG  >  NC_000913/1333483‑1333753

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