Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,000,775	A→G	N144S (AAC→AGC)	xdhA →	putative xanthine dehydrogenase molybdenum‑binding subunit XdhA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,000,775	0	A	G	100.0%	198.0 / NA	54	N144S (AAC→AGC)	xdhA	putative xanthine dehydrogenase molybdenum‑binding subunit XdhA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/28);  total (26/28)

GCCGTTGCCATCGTCGTGGCCCGCGATGAACTCACGGCAGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGA  >  NC_000913/3000639‑3000898                                                                                                                                         |                                                                                                                            gCCGTTGCCATCGTCGTGGCCCGCGATGAACTCACGGCAGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCgg                                                                                                                          >  2:162383/1‑140 (MQ=255)  ccGTTGCCATCGTCGTGGCCCGCGATGAACTCACGGCAGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGt                                                                                                                         >  2:755750/1‑140 (MQ=255)     ttGCCATCGTCGTGGCCCGCGATGAACTCACGGCAGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGc                                                                                                                      >  2:468171/1‑140 (MQ=255)       gCCATCGTCGTGGCCCGCGATGAACTCACGGCAGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCaa                                                                                                                    <  2:157318/140‑1 (MQ=255)       gCCATCGTCGTGGCCCGCGATGAACTCACGGCAGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCaa                                                                                                                    <  2:525056/140‑1 (MQ=255)         cATCGTCGTGGCGCGCGATGAACTCACGGCAGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTAGCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACGAATCCATAGCGGTGGCAAtt                                                                                                                  <  1:573818/140‑1 (MQ=255)                   gCCCGCGATGAACTCACGGCAGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACa                                                                                                        <  1:755517/140‑1 (MQ=255)                       gcgATGAACTCACGGCAGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGc                                                                                                    >  2:109267/1‑140 (MQ=255)                        cgATGAACTCACGGCAGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCa                                                                                                   >  1:594763/1‑140 (MQ=255)                        cgATGAACTCACGGCAGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAAACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCa                                                                                                   >  1:136196/1‑140 (MQ=255)                             aaCTCACGGCAGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATg                                                                                              <  1:162383/140‑1 (MQ=255)                             aaCTCACGGCAGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATg                                                                                              <  1:528784/140‑1 (MQ=255)                                tCACGGCAGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGTCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCg                                                                                           <  1:755750/140‑1 (MQ=255)                                       aGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTa                                                                                    >  1:121819/1‑140 (MQ=255)                                         aaaaaGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAAt                                                                                  <  1:109267/140‑1 (MQ=255)                                          aaaaGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATg                                                                                 <  2:121819/140‑1 (MQ=255)                                            aaGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTa                                                                                    <  2:282972/135‑1 (MQ=255)                                            aaGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTa                                                                                    >  1:282972/1‑135 (MQ=255)                                               cggcgCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCcaa                                                                            <  1:850122/140‑1 (MQ=255)                                                ggcgCAATTGGTCAGCATTCAGTGGCAAGAATTACCCGTTATCGCCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCcaac                                                                           >  2:806868/1‑140 (MQ=255)                                                         ggTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCg                                                                  >  1:113816/1‑140 (MQ=255)                                                            cAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATg                                                               >  1:360356/1‑140 (MQ=255)                                                                 ttGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATgccgcc                                                          <  1:644956/140‑1 (MQ=255)                                                                   gAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGa                                                        <  1:292956/140‑1 (MQ=255)                                                                    aGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTCCTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGAc                                                       >  1:64058/1‑140 (MQ=255)                                                                    aGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGAc                                                       >  1:244546/1‑140 (MQ=255)                                                                    aGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGAc                                                       >  2:444463/1‑140 (MQ=255)                                                                        gCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTAcc                                                   <  1:444463/140‑1 (MQ=255)                                                                         cAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCa                                                  >  2:454232/1‑140 (MQ=255)                                                                         cAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCa                                                  >  1:853014/1‑140 (MQ=255)                                                                         cAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCa                                                  >  2:171668/1‑140 (MQ=255)                                                                              aTTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTAc                                             <  2:853014/140‑1 (MQ=255)                                                                               ttACCCGTTATCGCCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACa                                            <  1:806868/140‑1 (MQ=255)                                                                                  cccGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAggg                                         >  2:199202/1‑140 (MQ=255)                                                                                  cccGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAggg                                         >  2:473553/1‑140 (MQ=255)                                                                                   ccGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAgggg                                        <  1:171668/140‑1 (MQ=255)                                                                                    cGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCTATGCCGCCGACTACCAGGTACAGGGGc                                       <  2:594763/140‑1 (MQ=255)                                                                                    cGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGc                                       <  2:244546/140‑1 (MQ=255)                                                                                       tATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACt                                    >  1:576216/1‑140 (MQ=255)                                                                                            ccacGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAg                               >  1:446230/1‑140 (MQ=255)                                                                                              acGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGAc                             <  1:95642/140‑1 (MQ=255)                                                                                               cGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGAcc                            >  2:520075/1‑140 (MQ=255)                                                                                                gCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGAccc                           >  1:581072/1‑140 (MQ=255)                                                                                                 ccAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGAcccc                          >  2:93876/1‑140 (MQ=255)                                                                                                    gAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGAc                             >  2:569337/1‑134 (MQ=255)                                                                                                    gAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGAc                             <  1:569337/134‑1 (MQ=255)                                                                                                     aaGCGGCGCTGGCAGAAGACGCTGCACC‑ATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCa                                      >  2:396276/1‑123 (MQ=255)                                                                                                     aaGCGGCGCTGGCAGAAGACGCTGCACC‑ATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCa                                      <  1:396276/123‑1 (MQ=255)                                                                                                      aGCGGCGCTGGCAGAAGACGCTGCATCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTa                     <  1:113409/140‑1 (MQ=255)                                                                                                      aGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTa                     <  2:587554/140‑1 (MQ=255)                                                                                                               ggCAGAAGACGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACAtt            <  1:454232/140‑1 (MQ=255)                                                                                                                   gaaaaCGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTca        <  2:136196/140‑1 (MQ=255)                                                                                                                    aagaCGCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTcac       <  1:473553/140‑1 (MQ=255)                                                                                                                         gCTGCACCAATCCATAGCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGa  <  2:113816/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                            GCCGTTGCCATCGTCGTGGCCCGCGATGAACTCACGGCAGAAAAAGCGGCGCAATTGGTCAGCATTGAGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGA  >  NC_000913/3000639‑3000898

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