Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,169,693	C→T	pseudogene (475/475 nt)	gatR ←	pseudogene, repressor for gat operon, interrupted by IS3, split galactitol utilization operon repressor, fragment 2, split galactitol utilization operon repressor, interrupted

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,169,693	0	C	T	100.0%	139.9 / NA	45	pseudogene (475/475 nt)	gatR	pseudogene, repressor for gat operon, interrupted by IS3, split galactitol utilization operon repressor, fragment 2, split galactitol utilization operon repressor, interrupted
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/22);  total (23/22)

CGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGAT  >  NC_000913/2169573‑2169828                                                                                                                         |                                                                                                                                        cgttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAAAACAATATTACACCGTAATTAAtttt                                                                                                                       >  1:158508/1‑139 (MQ=255)                     ttGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGc                                                                                                   <  2:175441/139‑1 (MQ=255)                        tCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAAt                                                                                                     <  2:473473/134‑1 (MQ=255)                        tCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAAt                                                                                                     >  1:473473/1‑134 (MQ=255)                        tCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAg                                                                                                >  1:14479/1‑139 (MQ=255)                              gTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTaca                                                                                                                                     <  2:415546/96‑1 (MQ=255)                              gTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTaca                                                                                                                                     >  1:415546/1‑96 (MQ=255)                                  cGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGtc                                                                                                              <  2:188378/115‑1 (MQ=255)                                  cGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGtc                                                                                                              >  1:188378/1‑115 (MQ=255)                                  cGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCAcc                                                                                      >  2:518103/1‑139 (MQ=255)                                     tAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACctct                                                                                   <  2:605478/139‑1 (MQ=255)                                       aTCAAACTATTCTCTCTTGAACCAGGGAGCATCCGTAATGCCGGATGCGGTGTAAACACCATATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGc                                                                                                        >  2:277881/1‑116 (MQ=255)                                       aTCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGc                                                                                                        <  1:277881/116‑1 (MQ=255)                                       aTCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGa                                                                                 >  2:683629/1‑139 (MQ=255)                                            aCTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAg                                                                            <  1:175069/139‑1 (MQ=255)                                                ttCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATcc                                                                                         <  2:701610/122‑1 (MQ=255)                                                ttCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATcc                                                                                         >  1:701610/1‑122 (MQ=255)                                                 tctctcTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGtc                                                                                                              >  1:116012/1‑100 (MQ=255)                                                 tctctcTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGtc                                                                                                              <  2:116012/100‑1 (MQ=255)                                                 tctctcTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGcatca                                                                       <  2:14479/139‑1 (MQ=255)                                                  ctctctTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGcatcat                                                                      <  1:683629/139‑1 (MQ=255)                                                  ctctctTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGcatcat                                                                      <  1:518103/139‑1 (MQ=255)                                                               gggAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTaataataat                                                         >  2:3237/1‑139 (MQ=255)                                                               gggAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTaataataat                                                         >  2:621089/1‑139 (MQ=255)                                                               gggAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTaataataat                                                         >  1:490362/1‑139 (MQ=255)                                                               gggAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTaataataat                                                         >  1:561758/1‑139 (MQ=255)                                                                           aaTGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACctc                                                                                    <  2:288641/100‑1 (MQ=255)                                                                           aaTGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACctc                                                                                    >  1:288641/1‑100 (MQ=255)                                                                             tGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGa                                           <  2:158508/139‑1 (MQ=255)                                                                             tGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGa                                           <  1:3237/139‑1 (MQ=255)                                                                             tGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGa                                           <  2:490362/139‑1 (MQ=255)                                                                             tGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGa                                           <  2:561758/139‑1 (MQ=255)                                                                                   aTGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGgta                                     <  2:694930/139‑1 (MQ=255)                                                                                         tgtAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCt                               >  2:266454/1‑139 (MQ=255)                                                                                         tgtAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTACAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGa                                           >  2:65277/1‑127 (MQ=255)                                                                                         tgtAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTACAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGa                                           <  1:65277/127‑1 (MQ=255)                                                                                                   ttATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAAcc                     >  2:678371/1‑139 (MQ=255)                                                                                                   ttATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAAcc                     >  2:459957/1‑139 (MQ=255)                                                                                                    tATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGATACAGcatcat                                                                      >  1:134277/1‑89 (MQ=255)                                                                                                    tATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGATACAGcatcat                                                                      <  2:134277/89‑1 (MQ=255)                                                                                                                aGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGa                                           <  2:652163/104‑1 (MQ=255)                                                                                                                aGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGa                                           >  1:724415/1‑104 (MQ=255)                                                                                                                aGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGa                                           >  1:652163/1‑104 (MQ=255)                                                                                                                aGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGa                                           <  2:724415/104‑1 (MQ=255)                                                                                                                aGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAAt        >  2:472783/1‑139 (MQ=255)                                                                                                                      atatTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCAACATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGTGACGGTTAAATTTCGAt  <  1:459957/139‑1 (MQ=255)                                                                                                                         |                                                                                                                                        CGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAACCGCGACGGTTAAATTTCGAT  >  NC_000913/2169573‑2169828

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