Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	545,039	A→G	V85A (GTG→GCG)	allC ←	allantoate amidohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	545,039	0	A	G	100.0%	141.7 / NA	41	V85A (GTG→GCG)	allC	allantoate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/16);  total (25/16)

ACATACGGGAAGCGGCTGCCTTCTTCTTCTGCCATCGCCACCACTTCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCACGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACGTGTTTCCAGCCCGCTTGCTGCCATTCTTTTTTTAAATTGCTGCTGGGTTTCCAGCC  >  NC_000913/544892‑545186                                                                                                                                                    |                                                                                                                                                    aCATACGGGAAGCGGCTGCCTTCTTCTTCTGCCATCGCCACCACTTCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCgg                                                                                                                                                  >  2:32656/1‑151 (MQ=255)   aTACGGGAAGCGGCTGCCTTCTTCTTCTGCCATCGCCACCACTTCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGt                                                                                                                                                 >  2:567126/1‑150 (MQ=255)        ggAAGCGGCTGCCTTCTTCTTCTGCCATCGCCACCACTTCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGat                                                                                                                                           >  1:25228/1‑151 (MQ=255)        ggAAGCGGCTGCCTTCTTCTTCTGCACCCGCCACCACTCCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCACCCAGACAATTGCCAGCAACCCAGCAAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGat                                                                                                                                           <  1:103158/151‑1 (MQ=255)          acgCGGCTGCCTTCTTCTTCTGCCATCGCCACCACTTCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGatat                                                                                                                                         >  1:187213/3‑151 (MQ=255)              ggCTGCCTTCTTCTTCTGCCATCGCCACCACTTCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCg                                                                                                                                      >  2:642338/1‑150 (MQ=255)                 ggCCTTCTTCTTCTGCCATCGCCACCCCTTCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGGATCGATAGGCGaaac                                                                                                                                  <  1:115438/150‑3 (MQ=255)                   ccTTCTTCTTCTGCCATCGCCACCACTTCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGc                                                                                                                                >  1:209311/1‑151 (MQ=255)                      tcttcttctGCCATCGCCACCACTTCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCt                                                                                                                               <  1:608678/149‑1 (MQ=255)                                ccATCGCCCCCACTTCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACtt                                                                                                                    <  1:144112/150‑1 (MQ=255)                                 cATCGCCACCACTTCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTcc                                                                                                                  >  1:646496/1‑151 (MQ=255)                                      ccaccaCTTCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTgg                                                                                                             >  2:562367/1‑151 (MQ=255)                                          cacTTCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGatat                                                                                                         >  1:55140/1‑151 (MQ=255)                                            cTTCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTc                                                                                                       >  1:459404/1‑151 (MQ=255)                                                  ccGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATGCCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTCCCTGTGGATATTCGGTGc                                                                                                  <  1:567143/150‑1 (MQ=255)                                                   cgtacgtaGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGcc                                                                                                 >  2:86130/1‑150 (MQ=255)                                                       cgtaGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTc                                                                                             >  1:669805/1‑150 (MQ=255)                                                        gtggcggcgCGCCGTATTGCGTTTCCACCCAGTCAATTCCCGGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAg                                                                                           <  1:815/148‑1 (MQ=255)                                                              gcgcgcCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGAcc                                                                                     >  2:408101/1‑151 (MQ=255)                                                                 cgcCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGt                                                                                   >  1:440177/1‑150 (MQ=255)                                                                       aTTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAAtt                                                                            >  1:433877/1‑151 (MQ=255)                                                                         tGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTcc                                                                          >  1:137743/1‑151 (MQ=255)                                                                            gTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTcccc                                                                        >  1:145526/1‑150 (MQ=255)                                                                             ttttCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCAc                                                                      >  1:233832/1‑151 (MQ=255)                                                                             ttttCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCAc                                                                      >  2:324297/1‑151 (MQ=255)                                                                             ttttCAGCCAGTCAATTGCCACCCCCGCCGCCAGCGCGCCGATTTCCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTcc                                                                          <  2:611824/147‑1 (MQ=255)                                                                                    ccAGTCAATTGCCACCCCCGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATc                                                                <  2:335121/150‑1 (MQ=255)                                                                                           aTTGCCAGCCACCCCCCCAGCCCGCCAAATTGCCGGTCAAGGTAGCCGCCGTTAACCGCGGTATCGATATGCTAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACgt                                                         <  2:370186/150‑1 (MQ=255)                                                                                             tGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACGTGtt                                                      >  2:37119/1‑151 (MQ=255)                                                                                                         cgccAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACGTGTTTCCAGCCCGCtt                                          <  2:145132/151‑1 (MQ=255)                                                                                                           ccAGCGCGCCAAATTCCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACGTGTTTCCAGCCCGCTtg                                         <  2:462288/150‑1 (MQ=255)                                                                                                                 cgcCGAATTGCCCGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACGTGTTTCCAGCCCGCTTGCTGCCAt                                  >  1:631783/1‑151 (MQ=255)                                                                                                                           cccGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACGTGTTTCCAGCCCGCTTGCTGCCATTCttttttt                         >  1:54517/1‑150 (MQ=255)                                                                                                                           cccGTCAAGGTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACGTGTTTCCAGCCCGCTTGCTGCCATTCttttttt                         <  2:55138/150‑1 (MQ=255)                                                                                                                             cGTCAAGTTTACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACGTGTTTCCAGCCCGCTTGCTGCCATTCTTTTTTTaaa                      <  1:660599/151‑1 (MQ=255)                                                                                                                                     ttACCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACGTGTTTCCAGCCCGCTTGCTGCCATTCTTTTTTTAAATtgctgc               <  1:362377/150‑1 (MQ=255)                                                                                                                                       aCCGCCGTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACGTGTTTCCAGCCCGCTTGCTGCCATTCTTTTTTTAAATTGCTGCTgg            >  2:311632/1‑151 (MQ=255)                                                                                                                                        ccgccgTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACGTGTTTCCAGCCCGCTTGCTGCCATTCTTTTTTTAAATTGCTGCTggg           >  2:337166/1‑151 (MQ=255)                                                                                                                                           ccgTTAACCGCGGTATCGATATGCAAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACGTGTTTCCAGCCCGCTTGCTGCCATTCTTTTTTTAAATTGCTGCTGGGttt        <  2:406988/151‑1 (MQ=255)                                                                                                                                            cgTTAACCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACGTGTTTCCAGCCCGCTTGCTGCCATTCTTTTTTTAAATTGCTGCTGGGTTTc       >  2:672966/1‑151 (MQ=255)                                                                                                                                                 aCCGCGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACGTGTTTCCAGCCCGCTTGCTGCCATTCTTTTTTTAAATTGCTGCTGGGTTTCCAGcc  <  2:584149/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ACATACGGGAAGCGGCTGCCTTCTTCTTCTGCCATCGCCACCACTTCGACCGTACGTAGCGGCGCGCCGTATTGCGTTTTCAGCCAGTCAATTGCCAGCCACGCCGCCAGCGCGCCGAATTGCCCGTCAAGGTTACCGCCGTTAACCACGGTATCGATATGCGAACCGCTCAGAACCACTTCCTGTGGATATTCGGTGCCATTCAGGCGACCGTATAAATTCCCCACTTCATCGAAACGTGTTTCCAGCCCGCTTGCTGCCATTCTTTTTTTAAATTGCTGCTGGGTTTCCAGCC  >  NC_000913/544892‑545186

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