Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,951,372	A→G	pseudogene (813/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,372	0	A	G	100.0%	133.7 / NA	39	pseudogene (813/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 (23/16);  total (23/16)

AAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCT  >  NC_000913/3951253‑3951507                                                                                                                        |                                                                                                                                        aaaGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGAc                                                                                                                     >  2:440737/1‑140 (MQ=255)   aGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCg                                                                                                                   <  2:18856/140‑1 (MQ=255)   aGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGCGGGCGCGCGTTTTGATGACCGGGTGACCg                                                                                                                   <  2:548287/140‑1 (MQ=255)     ggCTGGGCGCAGTAGGAGCAGAGTACCCGCACTATCTTGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGc                                                                                                                 <  1:322790/140‑1 (MQ=255)                         gATTATCCGTACTATCTGGGCATTCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGcc                                                                                             >  1:344935/1‑140 (MQ=255)                          aTTATCCGTACTATCTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCca                                                                                            <  1:440737/140‑1 (MQ=255)                           ttATCCGTACTATCTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCcac                                                                                           >  1:383721/1‑140 (MQ=255)                             aTCCGTACTATCTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCacac                                                                                         <  1:144704/140‑1 (MQ=255)                               ccGTACTATCTGGGCGTGCTGGTGAGGCACGGCACCATAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGc                                                                                       <  1:318701/140‑1 (MQ=255)                               ccGTACTATCTGGGCATGCTGGGGGTGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGc                                                                                       <  1:514092/140‑1 (MQ=255)                               ccGTACTATCTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCCTTTTGATGTCCGGGCGACCGGCAAAATGAACACCTTCGCGCCACACGc                                                                                       >  2:298159/1‑140 (MQ=255)                                 gTACTATCTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCa                                                                                     >  1:333120/1‑140 (MQ=255)                                    cTATCTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAgtg                                                                                  >  1:394505/1‑140 (MQ=255)                                        cTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTAt                                                                              >  1:140292/1‑140 (MQ=255)                                        cTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTAt                                                                              >  1:67162/1‑140 (MQ=255)                                         tGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGc                                                                                       >  2:654633/1‑130 (MQ=255)                                         tGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGc                                                                                       <  1:654633/130‑1 (MQ=255)                                                 cTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGa                                                                     >  1:583731/1‑140 (MQ=255)                                                          cACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCCCGCCACACGCCAGTATTATCCCTATGGATATCGAccc                                                            >  1:252243/1‑140 (MQ=255)                                                              gCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCa                                                        >  2:662184/1‑140 (MQ=255)                                                                  cAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGaaa                                                    <  1:543897/140‑1 (MQ=255)                                                                     aGCGGCAAACTTAGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCAGATAGATATCGACCCGGCAGAAATGa                                                 >  1:252220/1‑140 (MQ=255)                                                                      gCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGGGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGaa                                                <  2:344935/140‑1 (MQ=255)                                                                             aCTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTg                                         <  2:167460/140‑1 (MQ=255)                                                                               ttCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCg                                       <  2:333120/140‑1 (MQ=255)                                                                                tCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGt                                      >  1:366714/1‑140 (MQ=255)                                                                                      tGCAGGAGTGTGACCCGCTGATTGCGGTGTGCGCGCGGTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGca                                <  1:298159/140‑1 (MQ=255)                                                                                            aGTGTGACCTGCTGATGGCCGTCGGCGCGCGTTTTGATGACCGGGTGATCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACAtgtg                          <  2:291856/140‑1 (MQ=255)                                                                                                tgACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGGAGAGATGAACAAACTGCGTCAGGCACATGTGGCAt                      >  2:279783/1‑140 (MQ=255)                                                                                                tgACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCAt                      >  2:72749/1‑140 (MQ=255)                                                                                                tgACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCAt                      >  2:492001/1‑140 (MQ=255)                                                                                                tgACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCAt                      >  1:344570/1‑140 (MQ=255)                                                                                                tgACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCAt                      <  1:662184/140‑1 (MQ=255)                                                                                                tgACCTGCTGATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCAt                      >  1:502651/1‑140 (MQ=255)                                                                                                     tgctgATCGCCGTGGTCGCGCGCTTCGATCACCGGGTTACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACaa                 <  1:205458/140‑1 (MQ=255)                                                                                                       ctgATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAgg               >  1:108379/1‑140 (MQ=255)                                                                                                       ctgATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGAGTGAGGCACGTGTGGCATTACAAgg               >  1:168397/1‑140 (MQ=255)                                                                                                       ctgATCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAACCTGCGTAAGGCACGTGTGGCATTACAAgg               >  1:278183/1‑140 (MQ=255)                                                                                                          aTCGCCGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCAGGCAGAAAGGAAGAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGa            >  2:205458/1‑140 (MQ=255)                                                                                                               cGTGGGCGCGCGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTaa       <  2:383721/140‑1 (MQ=255)                                                                                                                    gcgcgcGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGct  <  2:344570/140‑1 (MQ=255)                                                                                                                        |                                                                                                                                        AAAGGGCTGGGCGCAGTAGAAGCAGATTATCCGTACTATCTGGGCATGCTGGGGATGCACGGCACCAAAGCGGCAAACTTCGCGGTGCAGGAGTGTGACCTGCTGATCGCCGTGGGCGCACGTTTTGATGACCGGGTGACCGGCAAACTGAACACCTTCGCGCCACACGCCAGTGTTATCCATATGGATATCGACCCGGCAGAAATGAACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCT  >  NC_000913/3951253‑3951507

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