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%	140.6 / 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 (21/24);  total (21/24)

AGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAAC  >  NC_000913/2169566‑2169808                                                                                                                                |                                                                                                                    aGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGtaa                                                                                                           <  2:540067/138‑1 (MQ=255) aGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGtaa                                                                                                           >  1:540067/1‑138 (MQ=255)   cgttgcgttgTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGtaatta                                                                                                        >  2:233994/1‑139 (MQ=255)              cGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCcgc                                                                                             <  2:733270/139‑1 (MQ=255)                  tACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGtct                                                                                         >  2:572380/1‑139 (MQ=255)                           aTTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATg                                                                                >  2:1242284/1‑139 (MQ=255)                           aTTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATg                                                                                <  2:1280370/139‑1 (MQ=255)                                     gTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATc                                                                      <  2:976584/139‑1 (MQ=255)                                     gTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATc                                                                      <  2:1184511/139‑1 (MQ=255)                                     gTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATc                                                                      <  1:572380/139‑1 (MQ=255)                                         cGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCAc                                                                   >  2:7460/1‑138 (MQ=255)                                         cGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCAc                                                                   <  1:7460/138‑1 (MQ=255)                                              aTCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATcc                                                                     <  2:487065/131‑1 (MQ=255)                                              aTCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATcc                                                                     >  1:487065/1‑131 (MQ=255)                                              aTCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATc                                                                      <  2:1273069/130‑1 (MQ=255)                                              aTCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATc                                                                      >  1:1273069/1‑130 (MQ=255)                                                 aaaCTATTCTCTCATGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTATACCGTAATTAATTTTAGCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGacac                                                          <  2:868383/139‑1 (MQ=255)                                                     tATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGca                                                      <  1:474354/139‑1 (MQ=255)                                                      aTTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGcat                                                     <  1:385690/139‑1 (MQ=255)                                                                 aaCCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGcatca                                                   <  2:1102001/130‑1 (MQ=255)                                                                 aaCCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGcatca                                                   >  1:1102001/1‑130 (MQ=255)                                                                 aaCCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGtaata                                          >  2:392638/1‑139 (MQ=255)                                                                      gggAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTg                                                                                     >  1:745346/1‑91 (MQ=255)                                                                      gggAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTaataataat                                     >  1:985011/1‑139 (MQ=255)                                                                      gggAGAATCCGTAATGCCGGATGAGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTg                                                                                     <  2:745346/91‑1 (MQ=255)                                                                             tCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGcc                                                                              >  2:74823/1‑91 (MQ=255)                                                                             tCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGcc                                                                              <  1:74823/91‑1 (MQ=255)                                                                             tCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAAtt                              >  1:228374/1‑139 (MQ=255)                                                                                 tAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGcc                                                                              <  1:497610/87‑1 (MQ=255)                                                                                 tAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGcc                                                                              >  2:497610/1‑87 (MQ=255)                                                                                 tAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGa                                                             <  2:506047/104‑1 (MQ=255)                                                                                 tAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGa                                                             >  1:506047/1‑104 (MQ=255)                                                                                    tGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGa                       <  1:1242284/139‑1 (MQ=255)                                                                                       cGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTa                                                                                                  >  2:974162/1‑61 (MQ=255)                                                                                       cGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTa                                                                                                  <  1:974162/61‑1 (MQ=255)                                                                                       cGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGcatca                                                   <  2:527980/108‑1 (MQ=255)                                                                                       cGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGcatca                                                   >  1:527980/1‑108 (MQ=255)                                                                                       cGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCg                        <  2:671529/135‑1 (MQ=255)                                                                                       cGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCg                        >  1:671529/1‑135 (MQ=255)                                                                                         gATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATcc                                                                     <  1:662290/88‑1 (MQ=255)                                                                                         gATGCGGTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATcc                                                                     >  2:662290/1‑88 (MQ=255)                                                                                              ggTGTAAACACCTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTAc             <  2:985011/139‑1 (MQ=255)                                                                                                         cTTATCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAAc  >  2:1068821/1‑139 (MQ=255)                                                                                                            aTCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGa                                                             <  2:1213598/77‑1 (MQ=255)                                                                                                            aTCCGGCATACAGAACAATATTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGa                                                             >  1:1213598/1‑77 (MQ=255)                                                                                                                                |                                                                                                                    AGCGTTGCGTTGTCGATTACCACCAGATTGTCCACTGTTGCGTTAATCAAACTATTCTCTCTTGAACCAGGGAGAATCCGTAATGCCGGATGCGGTGTAAACACCTTATCCGGCATACAGAACAATACTACACCGTAATTAATTTTACCCGCGTCTTCTGCAATGCCAGTTTATCCACCTCTGACACAGCATCATCGGTAATAATAATGTCAATTTTTTCGATTGGTAGTACCTGATTAAAAC  >  NC_000913/2169566‑2169808

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