Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,744,362	G→A	100%	H104Y (CAC→TAC)	apt ←	adenine phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,744,362	0	G	A	100.0%	127.3 / NA	42	H104Y (CAC→TAC)	apt	adenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/20);  total (22/20)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTCAATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGC  >  CP000730/1744225‑1744496                                                                                                                                          |                                                                                                                                       ttcaattATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACAt                                                                                                                                     >  1:437141/1‑141 (MQ=255)     attATAAATGCAATACCTACTACGACACCGCCTAAGTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGtt                                                                                                                                 <  1:111138/141‑1 (MQ=255)          aaaTGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAAc                                                                                                                            >  1:246072/1‑141 (MQ=255)            aTGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATCAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACAt                                                                                                                          <  2:109941/141‑1 (MQ=255)            aTGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACAt                                                                                                                          <  2:3724/141‑1 (MQ=255)                aaTACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGt                                                                                                                      <  2:456219/141‑1 (MQ=255)                  tACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTAc                                                                                                                    <  1:416775/141‑1 (MQ=255)                   aCCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGtctc                                                                                                                   >  2:421396/1‑138 (MQ=255)                   aCCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTCAGACCGAACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACCTTGTTAAAACATTTGTAcc                                                                                                                   >  2:421477/1‑141 (MQ=255)                   aCCTACTACCATACCGCCTAATTTTTCAACTAACTCTATTGCTGCTTCAATCGTCCCACCAGTAGCTAATAAATCATCTGTAATTAACACACTTTGACCTGGTTTAATTGCATCTTTGTACAGTGTTAAAACATTTGTAcc                                                                                                                   <  1:181221/141‑1 (MQ=255)                     ctactaCGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTCTACATTGTTAAAACATTTGTACCaa                                                                                                                 >  2:131651/1‑140 (MQ=255)                      tactacGCTACCGCCTAATGTGTCAACCAATTTTATTGCTGCTTCAATCGTACCACCCGTAGCTCATAAATCATCTGTAACTAACACACTTTTACCTGGTTTAATTGCATCCTTGTACATTGTTAAAACCTTTGTACCata                                                                                                                >  2:180890/1‑141 (MQ=255)                        ctacGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATAtt                                                                                                              >  1:321423/1‑141 (MQ=255)                              tACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCATACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTTTACATTGTTAAAACATTTGTACCATATTCTAGGt                                                                                                        >  1:46483/1‑141 (MQ=255)                                          ttttttttttttttATTGCTGCTTCAACCGTACCACCAGTCGCTAATAAATCATCTGTAATTAACACCCCTTCCCCTGGTCCCATTGCATCTTTTTACATTTTTAAAACATTTTTACCATATTCTAGGTCTTACTCATAAc                                                                                            >  2:7606/11‑141 (MQ=255)                                              aaCTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAAt                                                                                        >  1:354519/1‑141 (MQ=255)                                                          tgctgcTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAgg                                                                            <  2:205394/141‑1 (MQ=255)                                                            ttgctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTa                                                                          <  2:448943/140‑1 (MQ=255)                                                            ctgctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTa                                                                          >  1:345683/1‑141 (MQ=255)                                                              gctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTCATTGCATCTTTGTACCTTGTTACACCATTTGTACCATATTCTAAGTCATTCTCATAACGACTGACTTGACGAGGGAAt                                                                        >  2:240693/1‑141 (MQ=255)                                                                        tccgccccgTAGCTAATAAATAATCTGTAATTAACACACGTTGGCCTGGTTTAATTGCATCTTTGTACATTGTTAACACATTTGTACTATATTATAGGTTATACTCATAACGTATGACTTCACGAGGTAATTTCCcttctt                                                              <  2:218717/133‑1 (MQ=255)                                                                           caccgGTAGCTAATAAATCAGCTGTAATTAACACTCGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATGCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTc                                                           <  1:437898/141‑1 (MQ=255)                                                                             ccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTa                                                         >  2:44497/1‑141 (MQ=255)                                                                                  aGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTACTTTTCTAACAgg                                                    <  2:329251/141‑1 (MQ=255)                                                                                    cTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATCACTTAACGAGGTAATAACCCTTCTTTTCTAACAggag                                                  >  2:360121/1‑139 (MQ=255)                                                                                       ataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTCCCCTTCTTTTCTAACAGGTCCaa                                               >  2:220579/1‑141 (MQ=255)                                                                                         aaaTCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAg                                             <  2:439686/141‑1 (MQ=255)                                                                                           atcatcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTGACAGGTGCAAAGcc                                           <  1:241377/141‑1 (MQ=255)                                                                                             gatcTGTAATTTACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCaa                                         <  2:150872/140‑1 (MQ=255)                                                                                             catcTGAAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCaa                                         >  2:255257/1‑141 (MQ=255)                                                                                                  gTAATTAACACACGTTGACCTGGGGTAAGTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATcccc                                    >  2:270037/1‑141 (MQ=255)                                                                                                   tAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATccccc                                   >  2:242581/1‑140 (MQ=255)                                                                                                      ttAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTg                                <  1:47112/141‑1 (MQ=255)                                                                                                             ccGTTGACCTGGTTCAATTGCATCTTTGTACATTGTTAAAACATTGGTACTATATTCTATGTCATACTCATAGCGAATGACTTCTCGGTCTAAGTTCCCGTCTTCTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGc                         <  2:2806/140‑1 (MQ=255)                                                                                                             acGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGc                         <  1:19833/141‑1 (MQ=255)                                                                                                              cGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCt                        <  2:437141/141‑1 (MQ=255)                                                                                                                   aCCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAgg                   <  1:190707/141‑1 (MQ=255)                                                                                                                      tGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACa                <  1:181234/141‑1 (MQ=255)                                                                                                                             aTTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTATGTCATACTCATAACGAATGACTTCACAAGGTAATATCCCTTCTTTTCTAACAGATGCAAAGCCAATCCCCAGTGACTAAGCTACAGGACAGCCAATg         >  2:353695/1‑141 (MQ=255)                                                                                                                             aTTGCATCTTTGTACATTGGTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCCCGAGGTAATTACCCTTCTTTTCTAACAGTTGCAAAGCCAAACCCCATTGAATAAGCTACAGGACAGCCAATg         >  2:353569/1‑141 (MQ=255)                                                                                                                                gCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATa      >  1:352487/1‑141 (MQ=255)                                                                                                                                    cTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGc  >  1:155234/1‑141 (MQ=255)                                                                                                                                          |                                                                                                                                       TTCAATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGC  >  CP000730/1744225‑1744496

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑