Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,842,455	C→T	*445* (TAG→TAA)	adeQ ←	adenine permease, high affinity, adenine:H+ symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,842,455	0	C	T	100.0%	127.6 / NA	41	*445* (TAG→TAA)	adeQ	adenine permease, high affinity, adenine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/19);  total (22/19)

TCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCGGTGAT  >  NC_000913/3842342‑3842589                                                                                                                  |                                                                                                                                       tCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTg                                                                                                               >  2:606518/1‑139 (MQ=255)            aaaCTAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCt                                                                                                    <  1:580897/139‑1 (MQ=255)                 aaccaCAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGac                                                                                               <  2:402977/139‑1 (MQ=255)                  accaCAAACACTACCGCGAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGacc                                                                                              >  1:524865/1‑139 (MQ=255)                         acacTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGAccaccacca                                                                                       >  2:753053/1‑139 (MQ=255)                           acTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcaca                                                                                     <  1:753053/139‑1 (MQ=255)                             tACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                                   <  1:606518/139‑1 (MQ=255)                                  ccAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGt                                                                              >  1:607000/1‑139 (MQ=255)                                       gTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGaccacca                                                                                          <  2:525436/122‑1 (MQ=255)                                       gTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGaccacca                                                                                          >  1:525436/1‑122 (MQ=255)                                         ggCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCGACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATc                                                                       >  2:507641/1‑139 (MQ=255)                                         ggCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATc                                                                       >  1:9448/1‑139 (MQ=255)                                              gCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGcc                                                                  <  2:779910/139‑1 (MQ=255)                                                   aTCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAgcgc                                                             <  2:12809/139‑1 (MQ=255)                                                    tCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATATTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGcc                                                            >  2:349169/1‑139 (MQ=255)                                                        gATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGt                                                        >  1:674008/1‑139 (MQ=255)                                                            aTGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGcc                                                            <  1:273556/131‑1 (MQ=255)                                                            aTGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGcc                                                            >  2:273556/1‑131 (MQ=255)                                                              gAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTaa                                                                                                                                     <  1:101810/56‑1 (MQ=255)                                                              gAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTaa                                                                                                                                     >  2:101810/1‑56 (MQ=255)                                                              gAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATccacca                                                                                                                              >  1:712322/1‑63 (MQ=255)                                                              gAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATccacca                                                                                                                              <  2:712322/63‑1 (MQ=255)                                                              gAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATAc                                                  >  1:680593/1‑139 (MQ=255)                                                                   ctcCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCa                                             >  2:840775/1‑139 (MQ=255)                                                                         cGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGAt                                                                        <  2:148509/106‑1 (MQ=255)                                                                         cGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGAt                                                                        >  1:148509/1‑106 (MQ=255)                                                                         cGTAATCAAAACTCATGATATCTGGATACCCTCGCCAGTTATTCATCCACCAGATTAATCTTCAGTGCAAACGGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGc                                       >  2:178871/1‑139 (MQ=255)                                                                                       aTGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAg                         >  1:146118/1‑139 (MQ=255)                                                                                         gATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGcc                       <  2:607000/139‑1 (MQ=255)                                                                                            atCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGc                                                      <  1:106412/105‑1 (MQ=255)                                                                                            atCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGc                                                      >  2:106412/1‑105 (MQ=255)                                                                                                     cccTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATccc           >  2:390650/1‑139 (MQ=255)                                                                                                       cTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCtt         <  2:674008/139‑1 (MQ=255)                                                                                                         cGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGc                 >  1:184911/1‑129 (MQ=255)                                                                                                         cGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGc                 <  2:184911/129‑1 (MQ=255)                                                                                                           ccAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACa                                                                                                         <  2:473710/39‑1 (MQ=255)                                                                                                           ccAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACa                                                                                                         >  1:473710/1‑39 (MQ=255)                                                                                                           ccAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                                   <  1:570017/61‑1 (MQ=255)                                                                                                           ccAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                                   >  2:570017/1‑61 (MQ=255)                                                                                                              gTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCGGTGAt  <  1:390650/139‑1 (MQ=255)                                                                                                              gTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCGGTGAt  <  1:178871/139‑1 (MQ=255)                                                                                                               ttattaATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGccc                                                           <  2:591737/81‑1 (MQ=255)                                                                                                               ttattaATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGccc                                                           >  1:591737/1‑81 (MQ=255)                                                                                                                  |                                                                                                                                       TCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCGGTGAT  >  NC_000913/3842342‑3842589

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