Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,000,899	A→G	E185E (GAA→GAG)	xdhA →	xanthine dehydrogenase, molybdenum binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,000,899	0	A	G	100.0%	130.9 / NA	38	E185E (GAA→GAG)	xdhA	xanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/17);  total (21/17)

AAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACGAGTCATCAAACCATTTGTCGG  >  NC_000913/3000754‑3001045                                                                                                                                                  |                                                                                                                                                   aagaCGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGt                                                                                                                                               >  1:649817/1‑151 (MQ=255)      gCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACGTGGAGAGCGTAACAt                                                                                                                                          >  1:811619/1‑151 (MQ=255)      gCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACAt                                                                                                                                          <  2:646126/151‑1 (MQ=255)         gCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGc                                                                                                                                       <  2:767031/151‑1 (MQ=255)            ccAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTg                                                                                                                                     <  2:600865/150‑1 (MQ=255)                                   aCTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTcgcg                                                                                                              <  1:822548/150‑1 (MQ=255)                                   aCTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTcgcg                                                                                                              >  2:435581/1‑150 (MQ=255)                                   aCTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTcgcg                                                                                                              <  2:308028/150‑1 (MQ=255)                                      gAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAAt                                                                                                           <  2:47231/150‑1 (MQ=255)                                       aaaCAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTa                                                                                                         <  1:155229/151‑1 (MQ=255)                                                aCGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGttt                                                                                                >  2:102346/1‑151 (MQ=255)                                                      tCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGc                                                                                           >  2:571614/1‑150 (MQ=255)                                                          cGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCAccc                                                                                       >  2:472388/1‑150 (MQ=255)                                                           gggTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCAccc                                                                                       >  1:850767/1‑149 (MQ=255)                                                           gggTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAg                                                                                     <  1:783947/151‑1 (MQ=255)                                                                           acaaTCGATGCCGCCGACTACCTAGTACAGGGGCACTATCAGACCCTCGATATTCAACATTGTCACATGGAGAGCCTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTc                                                                     <  1:47609/151‑1 (MQ=255)                                                                           acaaTCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGTTGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACAtt                                                                         >  2:678034/1‑147 (MQ=255)                                                                             aaTCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTcgc                                                                   >  1:107024/1‑151 (MQ=255)                                                                              aTCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCa                                                                                                      >  1:266349/1‑115 (MQ=255)                                                                              aTCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCa                                                                                                      <  2:266374/115‑1 (MQ=255)                                                                              aTCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTAGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCATCACATTGTTcgcc                                                                  >  1:665062/1‑151 (MQ=255)                                                                                cGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTcgccgc                                                                >  1:645814/1‑151 (MQ=255)                                                                                  aTGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGca                                                                             >  1:443495/1‑136 (MQ=255)                                                                                  aTGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGca                                                                             <  2:443529/136‑1 (MQ=255)                                                                                    gccgccGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGtgg                                                            <  1:422347/151‑1 (MQ=255)                                                                                       gccgACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGtggttg                                                         <  2:443685/151‑1 (MQ=255)                                                                                                cAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGgcgc                                                >  1:340121/1‑151 (MQ=255)                                                                                                  ggTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTg                                              >  2:310270/1‑151 (MQ=255)                                                                                                        ggggCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTg                                              >  2:136999/1‑145 (MQ=255)                                                                                                        ggggCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGa                                            <  1:136976/147‑1 (MQ=255)                                                                                                        ggggCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATAtt                                        <  1:653885/151‑1 (MQ=255)                                                                                                                      cccccGTTATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGt                          <  1:199721/151‑1 (MQ=255)                                                                                                                            ttATTCAACATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACGAGt                    >  1:118507/1‑151 (MQ=255)                                                                                                                                   aCATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACGAGTCATCaaa             <  1:784321/151‑1 (MQ=255)                                                                                                                                    cATTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACGAGTCATCAAAc            >  2:646503/1‑151 (MQ=255)                                                                                                                                     aTTGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACGAGTCATCAAAcc           >  1:690447/1‑151 (MQ=255)                                                                                                                                      ttGTCACATGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACGAGTCATCAAAcc           >  1:249550/1‑150 (MQ=255)                                                                                                                                              tGGAGAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACGAGTCATCAAACCATTTGTCgg  >  2:439935/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                   AAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGCGCTGGATATTCCCTGGTCATGCGTACGAGTCATCAAACCATTTGTCGG  >  NC_000913/3000754‑3001045

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