Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,951,311	A→G	pseudogene (752/1645 nt)	ilvG →	acetolactate synthase II subunit IlvG, N‑terminal fragment (pseudogene)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,951,311	0	A	G	100.0%	181.9 / NA	50	pseudogene (752/1645 nt)	ilvG	acetolactate synthase II subunit IlvG, N‑terminal fragment (pseudogene)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/20);  total (30/20)

GCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATC  >  NC_000913/3951193‑3951444                                                                                                                       |                                                                                                                                      gcagTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACtt                                                                                                                  >  2:168262/1‑140 (MQ=255) gcagTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACtt                                                                                                                  >  1:665775/1‑140 (MQ=255)   agTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTcg                                                                                                                <  2:182950/140‑1 (MQ=255)        cGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTg                                                                                                           <  1:822184/140‑1 (MQ=255)           cTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAg                                                                                                        >  1:342490/1‑140 (MQ=255)             ttGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGa                                                                                                      >  1:188747/1‑140 (MQ=255)               gCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTGCGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTcg                                                                                                                >  1:466506/1‑128 (MQ=255)               gCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTcg                                                                                                                <  2:466506/128‑1 (MQ=255)               gCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAgt                                                                                                    >  2:152658/1‑140 (MQ=255)                                          gcctgccACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCAc                                                                         <  1:152658/140‑1 (MQ=255)                                          gcctgccACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCAc                                                                         >  1:776184/1‑140 (MQ=255)                                          gcctgccACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCAc                                                                         <  1:767479/140‑1 (MQ=255)                                                     gTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGAc                                                              >  1:296646/1‑140 (MQ=255)                                                         gCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCggg                                                          >  2:232135/1‑140 (MQ=255)                                                          cTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGt                                                         <  2:665775/140‑1 (MQ=255)                                                             aaaGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGAc                                                      >  2:399405/1‑140 (MQ=255)                                                             aaaGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGAc                                                      >  1:883952/1‑140 (MQ=255)                                                                gggCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCgg                                                   >  2:368452/1‑140 (MQ=255)                                                                gggCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCgg                                                   >  1:432975/1‑140 (MQ=255)                                                                  gCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCa                                                 <  2:342490/140‑1 (MQ=255)                                                                   cTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAGACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCaa                                                >  2:649466/1‑140 (MQ=255)                                                                            gTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAacac                                       >  1:800489/1‑140 (MQ=255)                                                                                aaGCAGCTTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTc                                   >  2:573240/1‑140 (MQ=255)                                                                                     gATTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGcc                              >  1:424518/1‑140 (MQ=255)                                                                                      aTTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCca                             <  1:649466/140‑1 (MQ=255)                                                                                      aTTATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCca                             <  2:188747/140‑1 (MQ=255)                                                                                       ttATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCcac                            >  1:398043/1‑140 (MQ=255)                                                                                       ttATCCGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCcac                            >  1:451480/1‑140 (MQ=255)                                                                                           gcGTGCTATTTGTGCATGCTGCGGACGCGCGGATGCAGAGCGGCAAACTTCGCGGTGCAGGGGTGTGGCCTGCTGATCGCCGTGGGCCATCGTTTTGACGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCGCACGc                        <  2:690700/139‑1 (MQ=255)                                                                                           ccGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGc                        <  2:432975/140‑1 (MQ=255)                                                                                           ccGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGc                        <  2:776184/140‑1 (MQ=255)                                                                                           ccGTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGc                        >  1:586760/1‑140 (MQ=255)                                                                                           ccGTACTATCTGGGCATGCGGGGGATGCGCGGCCCCAAAGCGGCAAACTGCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGc                        <  2:424518/140‑1 (MQ=255)                                                                                             gTACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCa                      <  1:368452/140‑1 (MQ=255)                                                                                              tACTATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAg                     <  1:399405/140‑1 (MQ=255)                                                                                                 tATCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAgtgt                  <  2:800489/140‑1 (MQ=255)                                                                                                  aTCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGtt                 <  1:518299/140‑1 (MQ=255)                                                                                                  aTCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGtt                 >  1:287300/1‑140 (MQ=255)                                                                                                   tCTGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTa                <  2:451480/140‑1 (MQ=255)                                                                                                     tGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTa                <  2:795680/138‑1 (MQ=255)                                                                                                     tGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTa                >  1:795680/1‑138 (MQ=255)                                                                                                     tGGGCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATc              <  2:287300/140‑1 (MQ=255)                                                                                                      gggCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATcc             >  1:620563/1‑140 (MQ=255)                                                                                                       ggCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATccc            >  1:690700/1‑139 (MQ=255)                                                                                                        gCATGCTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCat           >  2:154290/1‑140 (MQ=255)                                                                                                             cTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAgtg                   <  2:122615/127‑1 (MQ=255)                                                                                                             cTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAgtg                   >  1:122615/1‑127 (MQ=255)                                                                                                             cTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGa      >  1:1021159/1‑140 (MQ=255)                                                                                                             cTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGa      >  2:941755/1‑140 (MQ=255)                                                                                                             cTGGGGATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGa      >  2:707466/1‑140 (MQ=255)                                                                                                                gggATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGatat   <  2:1021159/140‑1 (MQ=255)                                                                                                                 ggATGCGCGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATc  >  2:427493/1‑140 (MQ=255)                                                                                                                       |                                                                                                                                      GCAGTTCCGGCTTTGCGTGAATTTCTCGCTGCCACAAAAATGCCTGCCACCTGTACGCTGAAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATC  >  NC_000913/3951193‑3951444

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