Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	40,676	G→A	T419I (ACC→ATC)	caiT ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	40,676	0	G	A	100.0%	149.0 / NA	44	T419I (ACC→ATC)	caiT	putative transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/25);  total (19/25)

TACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGGTGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTTGATGATGTTTTTATCT  >  NC_000913/40551‑40810                                                                                                                              |                                                                                                                                       tACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCa                                                                                                                           <  1:284257/139‑1 (MQ=255) tACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCa                                                                                                                           <  1:980105/139‑1 (MQ=255)             cAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATgaagaa                                                                                                               <  1:437077/139‑1 (MQ=255)              aaCCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGTAATAAAg                                                                                                                             <  2:993927/124‑1 (MQ=255)              aaCCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGTAATAAAg                                                                                                                             >  1:993927/1‑124 (MQ=255)               aCCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGAGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGc                                                                                                             >  1:693107/1‑139 (MQ=255)                ccAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTCGACATCGCCAGGGTATAAGAGCACGCGTTAGCCAGCGTAACGATGGCAATCAAGCAGAGGATGAAGAAGcc                                                                                                            >  1:706093/1‑139 (MQ=255)                         aCCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACAtggt                                                                                                   >  1:478240/1‑139 (MQ=255)                              ccGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCgg                                                                                              >  1:1038619/1‑139 (MQ=255)                              ccGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCgg                                                                                              <  2:189023/139‑1 (MQ=255)                                gATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTg                                                                                            <  2:693107/139‑1 (MQ=255)                                 aTACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGc                                                                                           >  2:362960/1‑139 (MQ=255)                                  tACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCt                                                                                          <  2:751717/139‑1 (MQ=255)                                         cagcagAGGTGGTTCTTCACCATCGCGTACTTCGCGCCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGc                                                                                   >  1:593015/1‑139 (MQ=255)                                                      tcttcACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAgg                                                                      <  1:1070740/139‑1 (MQ=255)                                                       cttcACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGt                                                                     <  2:1013640/139‑1 (MQ=255)                                                          cACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCa                                                                                                                           >  1:853759/1‑82 (MQ=255)                                                          cACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCa                                                                                                                           <  2:853759/82‑1 (MQ=255)                                                              aTCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAatg                                                              >  1:761450/1‑139 (MQ=255)                                                                   gTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGgc                                                         <  2:203780/139‑1 (MQ=255)                                                                      cTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGgcgcg                                                      <  1:1071969/139‑1 (MQ=255)                                                                           cgGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCa                                                                                                                           >  2:470023/1‑65 (MQ=255)                                                                           cgGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCa                                                                                                                           <  1:470023/65‑1 (MQ=255)                                                                             gCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCcaca                                               <  2:593015/139‑1 (MQ=255)                                                                                 gTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGt                                           <  2:478240/139‑1 (MQ=255)                                                                                         cGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCg                                   >  2:847169/1‑139 (MQ=255)                                                                                           ccAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGAt                                 <  2:365475/139‑1 (MQ=255)                                                                                           ccAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGAt                                 <  2:814477/139‑1 (MQ=255)                                                                                           ccAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGAt                                 <  2:761450/139‑1 (MQ=255)                                                                                            cAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATc                                >  2:17656/1‑139 (MQ=255)                                                                                              gggTATAAGAGCAGGCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAg                              >  1:724568/1‑139 (MQ=255)                                                                                                         cAGGCGTTAACCAGCGTAACGATGGCAATAAAGCa                                                                                                                           <  2:968307/35‑1 (MQ=25)                                                                                                         cAGGCGTTAACCAGCGTAACGATGGCAATAAAGCa                                                                                                                           >  1:968307/1‑35 (MQ=25)                                                                                                           ggCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACAtggtgg                                                                                                 >  1:529156/1‑59 (MQ=255)                                                                                                           ggCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACAtggtgg                                                                                                 <  2:529156/59‑1 (MQ=255)                                                                                                           ggCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCgg                                                                                              <  1:926664/62‑1 (MQ=255)                                                                                                           ggCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCgg                                                                                              >  2:926664/1‑62 (MQ=255)                                                                                                           ggCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAg                                                                                       >  2:626064/1‑69 (MQ=255)                                                                                                           ggCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAg                                                                                       <  1:626064/69‑1 (MQ=255)                                                                                                           ggCGTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTtg                 <  1:550940/139‑1 (MQ=255)                                                                                                              gTTAACCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTtgatg              >  2:341498/1‑139 (MQ=255)                                                                                                                 aaCCAGCGTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGACGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTTGatgatg           >  2:1025866/1‑139 (MQ=255)                                                                                                                        gTAACGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTTGATGATGTTTTtat    <  1:341498/139‑1 (MQ=255)                                                                                                                          aaCGATGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTTGATGATGTTTTtatct  <  2:1038619/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                       TACCGATAATGCCAACCAGAATTGACCAACCGATACGCACCAGCAGAGGTGGTTCTTCACCATCGCGTACTTCGCGGCAAGTGGACATCGCCAGGGTATAAGAGCAGGCGTTAACCAGCGTAACGGTGGCAATAAAGCAGAGGATGAAGAAGCCCCACATGGTGGCGGTGCTGAGTGGCAGAGCGGCCCAGGTTTCAATGATGGCGCGCGCCACACCGTACTGTTCGATCAGATTTGGAATGTTGATGATGTTTTTATCT  >  NC_000913/40551‑40810

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