Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,002,255	T→C	N637N (AAT→AAC)	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,002,255	0	T	C	100.0%	193.3 / NA	53	N637N (AAT→AAC)	xdhA	putative xanthine dehydrogenase molybdenum‑binding subunit XdhA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/27);  total (26/27)

CTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGAC  >  NC_000913/3002119‑3002389                                                                                                                                         |                                                                                                                                       ctgAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCa                                                                                                                                     <  2:947127/140‑1 (MQ=255) ctgAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCa                                                                                                                                     <  2:258492/140‑1 (MQ=255) ctgAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCa                                                                                                                                     <  2:567081/140‑1 (MQ=255)   gAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCAct                                                                                                                                   <  1:198184/140‑1 (MQ=255)       gCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCActgctg                                                                                                                               >  1:804415/1‑140 (MQ=255)         tCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGc                                                                                                                             >  1:1013551/1‑140 (MQ=255)         tCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGc                                                                                                                             >  2:769819/1‑140 (MQ=255)         tCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATAGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGc                                                                                                                             >  1:911156/1‑140 (MQ=255)           cATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAg                                                                                                                           <  2:657735/140‑1 (MQ=255)            aTCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGa                                                                                                                          <  1:769819/140‑1 (MQ=255)                 aaCCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTc                                                                                                                     >  2:203067/1‑140 (MQ=255)                  accaccACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCa                                                                                                                    >  1:707090/1‑140 (MQ=255)                         cTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTacac                                                                                                             <  2:804415/140‑1 (MQ=255)                           aaCCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACAcgg                                                                                                           <  1:541128/140‑1 (MQ=255)                           aaCCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACAcgg                                                                                                           <  1:203067/140‑1 (MQ=255)                            aCCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACAcggc                                                                                                          <  2:885026/140‑1 (MQ=255)                                aCCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCggaa                                                                                                      <  1:159090/140‑1 (MQ=255)                                 ccGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGTCAGAAGGTCAGGTACACGGCggaat                                                                                                     >  2:566455/1‑140 (MQ=255)                                  cGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACAcgg                                                                                                           <  1:49747/133‑1 (MQ=255)                                  cGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACAcgg                                                                                                           >  2:49747/1‑133 (MQ=255)                                      gTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCggaatgggaa                                                                                                <  2:663932/140‑1 (MQ=255)                                        ttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCggaatgggaatg                                                                                              >  1:470371/1‑140 (MQ=255)                                               tACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTg                                                                                       >  2:535378/1‑140 (MQ=255)                                                aCCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTgg                                                                                      >  2:114543/1‑140 (MQ=255)                                                    ttgttgATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTgg                                                                                  >  2:433666/1‑140 (MQ=255)                                                      gttgATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGgc                                                                                <  2:1013551/140‑1 (MQ=255)                                                       ttgATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGgcg                                                                               <  1:433666/140‑1 (MQ=255)                                                        tgATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGgcgc                                                                              <  1:874190/140‑1 (MQ=255)                                                        tgATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGgcgc                                                                              <  1:717356/140‑1 (MQ=255)                                                           tCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTAt                                                                           >  1:661015/1‑140 (MQ=255)                                                                cGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGaa                                                                      <  2:470371/140‑1 (MQ=255)                                                                                 tgCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGatcat                                                           <  1:362830/134‑1 (MQ=255)                                                                                 tgCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGatcat                                                           >  2:362830/1‑134 (MQ=255)                                                                                 tgCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGc                                                     >  2:170937/1‑140 (MQ=255)                                                                                 tgCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGc                                                     >  1:598434/1‑140 (MQ=255)                                                                                 tgCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCGGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGc                                                     >  1:903872/1‑140 (MQ=255)                                                                                         cACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAgcg                                             <  2:704891/140‑1 (MQ=255)                                                                                           ccATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAgcggc                                           >  1:949964/1‑140 (MQ=255)                                                                                           ccATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAgcggc                                           >  1:143546/1‑140 (MQ=255)                                                                                                  ccGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAgcggcg                                          <  1:706856/134‑1 (MQ=255)                                                                                                  ccGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAgcggcg                                          >  2:706856/1‑134 (MQ=255)                                                                                                      aTCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTg                                        >  1:645260/1‑132 (MQ=255)                                                                                                      aTCCTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTg                                        <  2:645260/132‑1 (MQ=255)                                                                                                        ccTCAACGTTCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAAcc                              >  2:895007/1‑140 (MQ=255)                                                                                                                tttaTGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCTTTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATctg                      <  2:496231/137‑1 (MQ=255)                                                                                                                ttCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATctg                      <  2:707090/140‑1 (MQ=255)                                                                                                                ttCATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATctg                      <  1:170937/140‑1 (MQ=255)                                                                                                                  cATGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGCCAGGTACACGGCGGAATGGGCATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATctgct                    >  1:845305/1‑140 (MQ=255)                                                                                                                    tGATTCAGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTgg                  >  2:31446/1‑140 (MQ=255)                                                                                                                          aGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATctg                      >  2:360805/1‑130 (MQ=255)                                                                                                                          aGGGCATATTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATctg                      <  1:360805/130‑1 (MQ=255)                                                                                                                                 atTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGcc     >  1:795238/1‑140 (MQ=255)                                                                                                                                 atTCTTAACCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGcc     <  2:143546/140‑1 (MQ=255)                                                                                                                                    cTTAACCCACTGCGGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGAc  <  2:903872/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       CTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGAC  >  NC_000913/3002119‑3002389

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