Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,744,362	G→A	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%	124.0 / NA	37	H104Y (CAC→TAC)	apt	adenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/22);  total (15/22)

TTTCAATTCAATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTT  >  CP000730/1744219‑1744505                                                                                                                                                |                                                                                                                                                tttCAATTCAATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACAtt                                                                                                                                             <  1:174389/148‑1 (MQ=255)        ttaattATAAGTGCAATACCGACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCCCCTGGATCTAATAATTCGTCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTaaag                                                                                                                                      <  2:395511/146‑2 (MQ=255)        tcaattATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTaaaa                                                                                                                                      <  2:350231/148‑1 (MQ=255)         caattATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAAc                                                                                                                                     <  2:393002/148‑1 (MQ=255)           attATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACAt                                                                                                                                   <  1:15922/148‑1 (MQ=255)           attATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACAt                                                                                                                                   <  1:297577/148‑1 (MQ=255)           attATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACAt                                                                                                                                   <  1:132577/148‑1 (MQ=255)           attATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACAt                                                                                                                                   <  1:398592/148‑1 (MQ=255)            ttATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACAtt                                                                                                                                  <  2:375290/148‑1 (MQ=255)              ataAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTg                                                                                                                                <  1:318275/148‑1 (MQ=255)                aaaTGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTa                                                                                                                              >  2:71559/1‑148 (MQ=255)                      aaTACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCatat                                                                                                                        >  2:132093/1‑148 (MQ=255)                      aaTACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCatat                                                                                                                        <  1:50009/148‑1 (MQ=255)                      aaTACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCatat                                                                                                                        >  1:237599/1‑148 (MQ=255)                      aaTACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCatat                                                                                                                        >  1:148954/1‑148 (MQ=255)                         aCCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCCATCGTACCACCACTCGCTCATCAACCATAAGTAAGTAACCCACGTCGCCCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCt                                                                                                                     <  1:83054/148‑1 (MQ=255)                                   aTACCGCCTAATTTTTTTTTTTTTTTTTTTGCTGCTTCAATCGTACCACCATTTGCTAATAAATCATGTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTTCCATATTCTAGATCATACt                                                                                                           >  2:290912/1‑148 (MQ=255)                                          cTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACg                                                                                                    <  2:325628/148‑1 (MQ=255)                                                   cAACTAATTTTATTGCTGCTTCAATCGTACCAGCAGTAGCTAATAAATCATCTGTAATTAAAACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTc                                                                                           >  2:272614/1‑148 (MQ=255)                                                        aaTTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAg                                                                                      >  1:196031/1‑148 (MQ=255)                                                                     ctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCcttc                                                                         <  1:411853/148‑1 (MQ=255)                                                                         aaTCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATGGTTAAAGCATTTGTTCCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCtttt                                                                     <  2:404194/148‑1 (MQ=255)                                                                                 caccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAgg                                                             <  1:206679/148‑1 (MQ=255)                                                                                                  tcatcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCa                                            <  1:252170/148‑1 (MQ=255)                                                                                                    atcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCAtt                                          >  2:350362/1‑148 (MQ=255)                                                                                                    atcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCAtt                                          >  2:182077/1‑148 (MQ=255)                                                                                                          aaTTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATaa                                    >  2:137356/1‑148 (MQ=255)                                                                                                          aaTTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATaa                                    >  2:105365/1‑148 (MQ=255)                                                                                                           aTTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAg                                   <  1:105820/148‑1 (MQ=255)                                                                                                           aTTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAg                                   <  2:140862/148‑1 (MQ=255)                                                                                                           aTTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATCTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAg                                   <  1:220775/148‑1 (MQ=255)                                                                                                               acacacGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAc                               >  1:269222/1‑148 (MQ=255)                                                                                                               acacacGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTATTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATGAGCTAc                               >  2:397226/1‑148 (MQ=255)                                                                                                                   acGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGa                           <  1:109553/148‑1 (MQ=255)                                                                                                                              gTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGAt                >  2:354031/1‑148 (MQ=255)                                                                                                                               tttAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATAATCATACCGAATGACTTCCCGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATa               <  1:290912/148‑1 (MQ=255)                                                                                                                                            ttGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCtt  >  2:365122/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                TTTCAATTCAATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTT  >  CP000730/1744219‑1744505

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