Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,843,762	T→G	S10R (AGT→CGT)	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,843,762	0	T	G	100.0%	214.1 / NA	58	S10R (AGT→CGT)	adeQ	adenine permease, high affinity; adenine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (24/34);  total (24/34)

TTTGTGCCGCGCCGAGGATTTGCGGGTTCACAAAAACGATGTACACCATGGTTAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGC  >  NC_000913/3843615‑3843908                                                                                                                                                    |                                                                                                                                                   tttGTGCCGCGCCGAGGATTTGCGGGTTCACAAAAACGATGTACACCATGGTTAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCAc                                                                                                                                                 >  2:258731/1‑151 (MQ=255)    gtgCCGCGCCGAGGATTTGCGGGTTCACAAAAACGATGTACACCATGGTTAAACAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGATTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGt                                                                                                                                               >  2:602539/1‑150 (MQ=255)     tgCCGCGCCGAGGATTTGCGGGTTCACAAAAACGATGTACACCATGGTTAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTa                                                                                                                                              >  1:425712/1‑150 (MQ=255)      gCCGCGCCGAGGATTTGCGGGTTCACAAAAACGATGTACACCATGGTTAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAAt                                                                                                                                            >  1:65832/1‑151 (MQ=255)          cgcCGAGGATTTGCGGGTTCACAAAAACGATGTACACCATGGTTAAAAAAGTGTTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGt                                                                                                                                        <  2:77328/151‑1 (MQ=255)          cgcCGAGGATTTGCGGGTTCACAAAAACGATGTACACCATGGTTAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGt                                                                                                                                        <  1:583662/151‑1 (MQ=255)               aGGATTTGCGGGTTCACAAAAACGATGTACACCATGGTTAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAATTAATCGCGAAAGCGCCCCTGATTCTTTACGAACGAATTCggatt                                                                                                                                      >  1:302862/1‑145 (MQ=255)                     ttcGGGTTCACAAAAACGATGTACACCATGGTTAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCatta                                                                                                                              <  2:862964/148‑1 (MQ=255)                       cGGGTTCACAAAAACGATGTACACCATGGTTAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTc                                                                                                                           <  1:477188/151‑1 (MQ=255)                             cacaAAAACGATGTACACCATGGTTAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCAttttt                                                                                                                     >  1:805414/1‑151 (MQ=255)                                aaaaaCGATGTACACCATGGTTAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTg                                                                                                                  <  2:780808/151‑1 (MQ=255)                                aaaaaCGATGTACACCATGGTTAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTg                                                                                                                  >  2:753375/1‑151 (MQ=255)                                      gATGTACACCATGGTTAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCAtt                                                                                                            >  2:897336/1‑151 (MQ=255)                                         gTACACCATGGTTAAAAAAGTGGTCAGCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTg                                                                                                           <  1:605381/149‑1 (MQ=255)                                           acacCATGGTTAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCGTTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTggag                                                                                                        <  1:377294/150‑1 (MQ=255)                                                   gTTAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTaa                                                                                               <  1:698118/151‑1 (MQ=255)                                                    ttAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAAt                                                                                              >  2:407309/1‑151 (MQ=255)                                                        aaaaGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCt                                                                                          <  1:409228/151‑1 (MQ=255)                                                         aaaGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCtt                                                                                         <  1:590786/151‑1 (MQ=255)                                                              ggTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCttt                                                                                    <  2:142133/151‑1 (MQ=255)                                                                  aTCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGt                                                                                 >  2:612434/1‑150 (MQ=255)                                                                    ccccGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTaa                                                                               >  2:467505/1‑150 (MQ=255)                                                                     cccGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTaaaa                                                                             <  1:94439/151‑1 (MQ=255)                                                                          aatcaatTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAAt                                                                        <  2:763239/151‑1 (MQ=255)                                                                           atcaatTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATg                                                                       <  2:118220/151‑1 (MQ=255)                                                                             caatTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGac                                                                     <  1:902895/151‑1 (MQ=255)                                                                                    gtgCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGacaaca                                                                 <  2:457741/148‑1 (MQ=255)                                                                                     tgCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACgg                                                             >  1:296634/1‑151 (MQ=255)                                                                                       cGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGtt                                                           <  2:226524/151‑1 (MQ=255)                                                                                       cGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGtt                                                           >  2:740863/1‑151 (MQ=255)                                                                                               ggtCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGAt                                                    <  1:239748/150‑1 (MQ=255)                                                                                               ggtCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGAt                                                    <  1:321371/150‑1 (MQ=255)                                                                                                gtCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATa                                                   >  1:85639/1‑150 (MQ=255)                                                                                                  cccATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAAtt                                                <  1:92084/151‑1 (MQ=255)                                                                                                                gTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAAtt                                  <  2:561225/151‑1 (MQ=255)                                                                                                                gTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAAtt                                  >  2:434723/1‑151 (MQ=255)                                                                                                                gTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAAtt                                  >  1:58948/1‑151 (MQ=255)                                                                                                                   tAAATAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGc                                >  1:96292/1‑150 (MQ=255)                                                                                                                    aaaTAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGcac                              >  2:538809/1‑151 (MQ=255)                                                                                                                    aaaTAATCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGca                               <  1:302657/150‑1 (MQ=255)                                                                                                                        aatCGCGAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACAtt                           <  2:333628/150‑1 (MQ=255)                                                                                                                            gcgAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATgg                      <  2:287083/151‑1 (MQ=255)                                                                                                                            gcgAAAGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATgg                      <  1:872094/151‑1 (MQ=255)                                                                                                                                aaGCGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGatat                  <  1:281190/151‑1 (MQ=255)                                                                                                                                   cGTCCCTTATTCATTACGAACGTAATCGGTATTGTCATTATTCATTTTTATGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGat                    <  2:483796/146‑1 (MQ=255)                                                                                                                                   cGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCt               <  2:578001/151‑1 (MQ=255)                                                                                                                                   cGTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCt               >  1:136753/1‑151 (MQ=255)                                                                                                                                    gTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCtt              <  2:904272/151‑1 (MQ=255)                                                                                                                                    gTCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCtt              <  2:203398/151‑1 (MQ=255)                                                                                                                                     tCCCTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCttt             >  1:217574/1‑151 (MQ=255)                                                                                                                                       ccTGATTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTta            >  1:708527/1‑150 (MQ=255)                                                                                                                                           aTTCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTaataat       >  1:550816/1‑151 (MQ=255)                                                                                                                                            ttCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAAtc      >  1:147358/1‑151 (MQ=255)                                                                                                                                             tCATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAAtct     >  2:572622/1‑151 (MQ=255)                                                                                                                                              cATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAAtct     <  2:760776/150‑1 (MQ=255)                                                                                                                                              cATTACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAAtct     <  2:479836/150‑1 (MQ=255)                                                                                                                                                 tACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTcgc  <  2:739637/150‑1 (MQ=255)                                                                                                                                                 tACGCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTcgc  <  2:821042/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TTTGTGCCGCGCCGAGGATTTGCGGGTTCACAAAAACGATGTACACCATGGTTAAAAAAGTGGTCATCCCCGCAATCAATTCTGTGCGGACGGTGGTCCCATGCTGAGGTAGTTTAAATAATCGCGAAAGCGTCCCTGATTCATTACTCACGTAATCGGTATTGTCATTATTCATTTTTTTGTCCATTGGAGGAGATTTAATCCCTTGCTTTTGTAAAACAAATGACAACAACGGTTCAGTGATAATTATCACATTTCAATTGCACATTAATGGATATTCTTTAATAATCTCGC  >  NC_000913/3843615‑3843908

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