Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,830,265	T→C	S441P (TCC→CCC)	xanP →	xanthine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,830,265	0	T	C	100.0%	127.1 / NA	36	S441P (TCC→CCC)	xanP	xanthine permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/20);  total (16/20)

CTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCC  >  NC_000913/3830143‑3830392                                                                                                                           |                                                                                                                                ctcGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATcgccgc                                                                                                                 <  2:990221/139‑1 (MQ=255)    gTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCgcggg                                                                                                              <  1:1032832/139‑1 (MQ=255)    gTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCgcggg                                                                                                              >  1:191479/1‑139 (MQ=255)    gTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCgcggg                                                                                                              >  2:139950/1‑139 (MQ=255)        gCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGt                                                                                                          <  1:266026/139‑1 (MQ=255)         ccGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTa                                                                                                         >  2:894862/1‑139 (MQ=255)          cGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTAt                                                                                                        <  1:402994/139‑1 (MQ=255)            cTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTa                                                                                                      >  1:1044776/1‑139 (MQ=255)                  cGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCa                                                                                                <  2:141289/139‑1 (MQ=255)                       ggCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTg                                                                                           <  1:930204/139‑1 (MQ=255)                       ggCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTg                                                                                           <  1:894862/139‑1 (MQ=255)                           aTTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCa                                                                                                <  2:738670/130‑1 (MQ=255)                           aTTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCa                                                                                                >  1:738670/1‑130 (MQ=255)                             tCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTATGCCCCTGAATGGCTGAAAAAACTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAAt                                                                                     <  1:647270/139‑1 (MQ=255)                                 aTTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACCGCCATCGTGCTGAATCTGa                                                                                 <  1:139950/139‑1 (MQ=255)                                 aTTATCGCGCTGTCGCTGGCGCTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTACCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGa                                                                                 >  2:635189/1‑139 (MQ=255)                                        cgcTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTccc                                                                          >  1:133848/1‑139 (MQ=255)                                          cTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCccac                                                                        <  2:15984/139‑1 (MQ=255)                                                                     tctcAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGgcggcg                                             >  2:469718/1‑139 (MQ=255)                                                                         agcagcCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTaa                                         >  2:684229/1‑139 (MQ=255)                                                                          gcagcCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTAAg                                        <  2:1044776/139‑1 (MQ=255)                                                                              ccGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTAAGATaa                                    >  1:967608/1‑139 (MQ=255)                                                                              ccGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTAAGATaa                                    >  1:344951/1‑139 (MQ=255)                                                                              ccGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTAAGATaa                                    >  2:402254/1‑139 (MQ=255)                                                                                      tttGCAGTTTGACCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTAAGATAATTAccgcc                            <  1:684229/139‑1 (MQ=255)                                                                                              ttGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCCGAAAAACAGTAATTCACTGGCGGCGGTAAGATAATTACCGCCGCAGtctc                    >  2:553660/1‑139 (MQ=255)                                                                                               tGCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCa                   <  2:750033/139‑1 (MQ=255)                                                                                                gCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCAc                  <  1:613985/139‑1 (MQ=255)                                                                                                gCCCCTGAATGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCAc                  <  1:553660/139‑1 (MQ=255)                                                                                                       aaTGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCAc                  <  1:777090/132‑1 (MQ=255)                                                                                                       aaTGGCTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCAc                  >  2:777090/1‑132 (MQ=255)                                                                                                            cTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCAc                  >  2:484151/1‑127 (MQ=255)                                                                                                            cTGAAAAACCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCAc                  <  1:484151/127‑1 (MQ=255)                                                                                                                aaaaCCTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTAAGATAATCACCGCCGCAGTCTCACCTCTTTAACCATTccc  <  1:469718/139‑1 (MQ=255)                                                                                                                     cTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTccc                                                                          >  2:364046/1‑62 (MQ=255)                                                                                                                     cTGCTCCCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTccc                                                                          <  1:364046/62‑1 (MQ=255)                                                                                                                           |                                                                                                                                CTCGTGAGCCGCTGAACCGTCGGGCGATTCTGATTATCGCGCTGTCGCTGGCGGTTGGTCTGGGCGTGTCTCAGCAGCCGCTGATTTTGCAGTTTGCCCCTGAATGGCTGAAAAACCTGCTCTCCTCCGGGATCGCCGCGGGCGGTATTACTGCCATCGTGCTGAATCTGATTTTCCCACCAGAAAAACAGTAATTCACTGGCGGCGGTAAGATAATTACCGCCGCAGTCTCACCTCTTTAACCATTCCC  >  NC_000913/3830143‑3830392

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: ‑