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%	155.5 / NA	48	*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/26);  total (22/26)

CTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCG  >  NC_000913/3842329‑3842584                                                                                                                               |                                                                                                                                  cTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGa                                                                                                                       <  1:76485/139‑1 (MQ=255)         ggATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCtt                                                                                                               >  1:849862/1‑139 (MQ=255)                 gCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCaa                                                                                                       <  2:1661688/139‑1 (MQ=255)                               accaCAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGacc                                                                                         >  2:1518110/1‑139 (MQ=255)                               accaCAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGacc                                                                                         >  2:844221/1‑139 (MQ=255)                                          tACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                              <  1:844221/139‑1 (MQ=255)                                           aCCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacacac                                                                             <  1:1562510/139‑1 (MQ=255)                                                    gTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGa                                                                    >  1:243470/1‑139 (MQ=255)                                                        cGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATcgc                                                                >  2:474661/1‑139 (MQ=255)                                                            cATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCa                                                            <  1:474661/139‑1 (MQ=255)                                                                aTCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCCCCACCACACACAGGTTCAGATCGCGCCAgcgc                                                        <  1:1518110/139‑1 (MQ=255)                                                                 tCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGACACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                              <  1:1052829/116‑1 (MQ=255)                                                                 tCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                              >  2:1052829/1‑116 (MQ=255)                                                                      aTCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                              <  2:734769/111‑1 (MQ=255)                                                                      aTCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                              >  1:734769/1‑111 (MQ=255)                                                                           gAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACagag                                                                                                 <  2:408474/87‑1 (MQ=255)                                                                           gAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACagag                                                                                                 >  1:408474/1‑87 (MQ=255)                                                                           gAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATAc                                             >  1:1474939/1‑139 (MQ=255)                                                                                ctcCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCt                                                                                               >  1:1655122/1‑84 (MQ=255)                                                                                ctcCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCt                                                                                               <  2:1655122/84‑1 (MQ=255)                                                                                  ccTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTc                                                                                                              >  2:1364455/1‑67 (MQ=255)                                                                                  ccTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTc                                                                                                              <  1:1364455/67‑1 (MQ=255)                                                                                   cTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCa                                                >  2:305039/1‑128 (MQ=255)                                                                                   cTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCa                                                <  1:305039/128‑1 (MQ=255)                                                                                   cTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCatga                                     <  2:849862/139‑1 (MQ=255)                                                                                     gCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCatgatg                                   <  2:1474939/139‑1 (MQ=255)                                                                                           tCAGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTAc                             >  2:894799/1‑139 (MQ=255)                                                                                             aGAACTCATGATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGa                           >  2:447550/1‑139 (MQ=255)                                                                                                      gATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCa                                                >  2:1635808/1‑109 (MQ=255)                                                                                                      gATATCTGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCa                                                <  1:1635808/109‑1 (MQ=255)                                                                                                            tGGAAACCCTCGCCAGTTATTAATCCCCCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGc            <  2:1594951/139‑1 (MQ=255)                                                                                                            tGGAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGc            <  1:447550/139‑1 (MQ=255)                                                                                                             ggAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCa           <  1:894799/139‑1 (MQ=255)                                                                                                             ggAAACCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCa           <  2:243470/139‑1 (MQ=255)                                                                                                                 aCCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCa                                                >  2:1640970/1‑98 (MQ=255)                                                                                                                 aCCCTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCa                                                <  1:1640970/98‑1 (MQ=255)                                                                                                                    cTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                              <  2:239525/65‑1 (MQ=255)                                                                                                                    cTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACcacaca                                                                              >  1:239525/1‑65 (MQ=255)                                                                                                                    cTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATAc                                             <  2:39829/98‑1 (MQ=255)                                                                                                                    cTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATAc                                             >  1:39829/1‑98 (MQ=255)                                                                                                                    cTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCtt    >  1:1480520/1‑139 (MQ=255)                                                                                                                    cTCGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCtt    <  1:1621218/139‑1 (MQ=255)                                                                                                                      cGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACCGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCg  >  2:1037502/1‑139 (MQ=255)                                                                                                                      cGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGcc                  <  2:575626/123‑1 (MQ=255)                                                                                                                      cGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGcc                  >  1:575626/1‑123 (MQ=255)                                                                                                                      cGCCAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCg  <  2:1480520/139‑1 (MQ=255)                                                                                                                         cAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATcgcg                                                               >  2:610451/1‑75 (MQ=255)                                                                                                                         cAGTTATTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATcgcg                                                               <  1:610451/75‑1 (MQ=255)                                                                                                                               |                                                                                                                                  CTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTCAGTGCAAACAGAGCTGCGACCACCACCACACACAGGTTCAGATCGCGCCAGCGCCCGGTGCATACTTTCATGATGCAGTACGACATAAAGCCGAGTGCAATCCCTTCG  >  NC_000913/3842329‑3842584

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