Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	142,598	A→G	F25L (TTT→CTT)	can ←	carbonic anhydrase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	142,598	0	A	G	100.0%	191.5 / NA	55	F25L (TTT→CTT)	can	carbonic anhydrase 2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/37);  total (18/37)

CAACATTACGGTGAACAAAGAGTTCGCCCGGCTCAAGACCGGTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAACCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGC  >  NC_000913/142472‑142708                                                                                                                               |                                                                                                               cAACATTACGGTGAACAAAGAGTTCGCCCGGCTCAAGACCGGTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTCCGCTTGTGCCACTTTCTCAAAAAGCCCggg                                                                                                          >  1:185441/1‑133 (MQ=255)   aCATTACGGTGAACAAAGAGTTCGCCCGGCTCAAGACCGGTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGAt                                                                                                        <  2:34074/133‑1 (MQ=255)   aCATTACGGTGAACAAAGAGTTCGCCCGACTCAAGACCGGTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGAt                                                                                                        <  2:177262/133‑1 (MQ=255)      ttACGGTGAACAAAGAGTTCGCCCGGCTCAAGACCGGTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCct                                                                                                     >  2:318575/1‑133 (MQ=255)       tACGGTGAACAAAGAGTTCGCCCGGCTCAAGACCGGTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCctc                                                                                                    >  2:465406/1‑133 (MQ=255)              aaCAAAGAGTTCGCCCGGCTCAAGACCGGTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGATTCTCAACAAGACCGGGATCCTCTTccacc                                                                                             >  1:224831/1‑133 (MQ=255)                      gTTCGCCCGGCTCAAGACCGGTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCAtttt                                                                                     >  1:23053/1‑133 (MQ=255)                        tCGCCCGGCTCAAGACCGGTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGGCCGAGATCCTCTACCACCAGCATTTTTg                                                                                   >  2:202070/1‑133 (MQ=255)                         cGCCCGGCTCTAGACCGGTTAAACCTTCCGCACGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGa                                                                                  <  2:130114/133‑1 (MQ=255)                         cGCCCGGCTCAAGGCCGGTTAATCGTTCTGTAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGa                                                                                  <  1:212653/133‑1 (MQ=255)                         cGCCCGGCTCAAGACCGGTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGa                                                                                  <  1:252484/133‑1 (MQ=255)                            ccGGCTCAAGACCGGTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCa                                                                               <  1:318575/133‑1 (MQ=255)                               gCTCAAGACCGGTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCCAAAAGCCCGGGACCTTATTCCACCAACATTGTTGACCATAg                                                                            >  2:252484/1‑133 (MQ=255)                               gCTCAAGACCGGTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAg                                                                            >  2:20819/1‑133 (MQ=255)                                   aaGACCGGTTCAACGATCTGCAGGAACGCGACTGTCGGAACTTCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCa                                                                        <  2:207462/133‑1 (MQ=255)                                      aCCGGTTAAAGGTTCTGCAGGAACTCGACTGTAGGAACATCGAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCAttg                                                                     <  2:185441/133‑1 (MQ=255)                                       ccGGTTAACCGTTCTGCACGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCAttgt                                                                    <  1:165974/133‑1 (MQ=255)                                         ggTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCAttgttg                                                                  >  1:275412/1‑133 (MQ=255)                                          gTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGTTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGc                                                                 >  2:369133/1‑133 (MQ=255)                                          gTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGc                                                                 >  1:85433/1‑133 (MQ=255)                                            tAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCtg                                                               <  2:23053/133‑1 (MQ=255)                                            tAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCtg                                                               <  2:292526/133‑1 (MQ=255)                                            tAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCtg                                                               <  2:337690/133‑1 (MQ=255)                                            tAAACGTTCTGCAGGAACGCGACTATCGAAACTTCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCtg                                                               <  2:174964/133‑1 (MQ=255)                                              aaCGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCtgat                                                             <  2:308734/133‑1 (MQ=255)                                                cGTTCTGCAGGAATGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCtgatga                                                           <  2:44115/133‑1 (MQ=255)                                                cGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCtgatga                                                           <  1:480121/133‑1 (MQ=255)                                                 gTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAg                                                          <  1:426786/133‑1 (MQ=255)                                                 gTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAg                                                          <  1:255324/133‑1 (MQ=255)                                                  ttCTGCAGAAACGCGCCTTTCGCAACATCCCATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAgt                                                         <  1:202070/133‑1 (MQ=255)                                                     tGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTa                                                      >  1:90379/1‑133 (MQ=255)                                                       cAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATc                                                    <  2:117475/133‑1 (MQ=255)                                                          gAACGCGACTGTCGGTACATCCAACCCATAGATAGCGCGGTCTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTAt                                                 <  2:181458/133‑1 (MQ=255)                                                          gAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTAt                                                 <  1:465406/133‑1 (MQ=255)                                                                    gTCGGAACCTCCAATACACAGACAGCGCGGCTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATg                                       <  1:135332/133‑1 (MQ=255)                                                                    gTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATg                                       <  2:340201/133‑1 (MQ=255)                                                                     tCGGAACATCCAATCAATAGAATGCGCGGTTTTTGCGCTTGTGTCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATgg                                      <  1:174050/133‑1 (MQ=255)                                                                            atccaatccaTAGAAAGCGCGGTTTTTGCGTTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTaa                               <  1:369133/133‑1 (MQ=255)                                                                            atccaatccaTAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTaa                               >  2:82695/1‑133 (MQ=255)                                                                            atccaatccaGAGAAGTCGCGGTTTTTGTGCTTGTGCCAGTCTCTCAAAAAGCTCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTaa                               <  2:224831/133‑1 (MQ=255)                                                                             tccaatccaTAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAAc                              <  1:20819/133‑1 (MQ=255)                                                                                 atccaTAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGAcc                          <  2:275412/133‑1 (MQ=255)                                                                                      tAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTaa                     >  2:402211/1‑133 (MQ=255)                                                                                      tAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTaa                     >  2:443689/1‑133 (MQ=255)                                                                                         aaaGCGCGGTCTTTGCGCTTGTGCGAGTTTCTCAAAAAGCCTGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTAACCa                  <  2:155486/133‑1 (MQ=255)                                                                                           aGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTAACCaaa                >  1:507018/1‑133 (MQ=255)                                                                                             cgcgGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTAACCAaata              >  2:342206/1‑133 (MQ=255)                                                                                              gcgGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGCTTAACGGCCTGTAACCAaataa             >  1:226599/1‑133 (MQ=255)                                                                                               cgGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTAACCAaataat            >  1:397284/1‑133 (MQ=255)                                                                                                    tttGCGCTTGTGCCATTTTCTCAACAAGTCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGt       <  2:226599/133‑1 (MQ=255)                                                                                                    tttGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGt       <  1:342206/133‑1 (MQ=255)                                                                                                    tttGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGt       <  2:266717/133‑1 (MQ=255)                                                                                                    tttGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGt       <  2:85433/133‑1 (MQ=255)                                                                                                    tttGCGCTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGt       <  2:89398/133‑1 (MQ=255)                                                                                                         gcTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGc  <  1:402211/133‑1 (MQ=255)                                                                                                         gcTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGc  <  1:443689/133‑1 (MQ=255)                                                                                                         gcTTGTGCCAGTTTCTCAAAAAGCCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGc  <  2:507018/133‑1 (MQ=255)                                                                                                                               |                                                                                                               CAACATTACGGTGAACAAAGAGTTCGCCCGGCTCAAGACCGGTTAAACGTTCTGCAGGAACGCGACTGTCGGAACATCCAATCCATAGAAAGCGCGGTTTTTGCGCTTGTGCCAGTTTCTCAAAAAACCCGGGATCCTCTTCCACCAGCATTTTTGACCATAGTGCATTGTTGCTGATGAGTGTATCTATGTCTTTCATGGAGGTTAACGACCTGTAACCAAATAATTACGTTTGGC  >  NC_000913/142472‑142708

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