Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,810	Δ1 bp	intergenic (‑42/+24)	pyrE ← / ← rph	orotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,809	0	G	.	100.0%	188.5 / NA	40	intergenic (‑41/+25)	pyrE/rph	orotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (18/22);  total (18/22)

GTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTC  >  NC_000913/3815661‑3815955                                                                                                                                                     |                                                                                                                                                   gTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAgc                                                                                                                                                   >  1:646264/1‑149 (MQ=255)     aagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACtc                                                                                                                                             <  1:9316/151‑1 (MQ=255)      agaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACtca                                                                                                                                            <  1:705611/151‑1 (MQ=255)        aaaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACtcatc                                                                                                                                          <  1:610293/151‑1 (MQ=255)                 cTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCtt                                                                                                                                 >  2:460835/1‑151 (MQ=255)                 cTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCtt                                                                                                                                 >  2:352965/1‑151 (MQ=255)                  ttttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCtt                                                                                                                                 >  2:787251/1‑150 (MQ=255)                  ttttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCtt                                                                                                                                 >  2:416314/1‑150 (MQ=255)                     tGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTaaa                                                                                                                              <  2:278989/150‑1 (MQ=255)                      gcgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTaaaa                                                                                                                             <  1:125038/150‑1 (MQ=255)                       cgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTaaaaa                                                                                                                            >  1:840188/1‑150 (MQ=255)                       cgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAAt                                                                                                                           <  1:250217/151‑1 (MQ=255)                                ttCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                                                                  >  1:629548/1‑151 (MQ=255)                                 tCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                                                                  <  2:161878/150‑1 (MQ=255)                                  cAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAg                                                                                                                <  1:610006/151‑1 (MQ=255)                                    gCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgcg                                                                                                              >  1:263945/1‑151 (MQ=255)                                     cGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgcgc                                                                                                             <  2:848847/151‑1 (MQ=255)                                           cTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCtt                                                                                                       >  1:725873/1‑151 (MQ=255)                                                         aaacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGATCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGACGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACa                                                                                          >  1:315839/3‑150 (MQ=255)                                                               gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgat                                                                                    <  1:199949/150‑1 (MQ=255)                                                                cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTGTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                                                                                   <  2:719678/150‑1 (MQ=255)                                                                  tgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcg                                                                                 >  1:93874/1‑150 (MQ=255)                                                                      aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgattc                                                                             <  2:93901/150‑1 (MQ=255)                                                                      aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTcc                                                                            >  1:458065/1‑151 (MQ=255)                                                                          gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTc                                                                        >  2:311175/1‑151 (MQ=255)                                                                                     aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGc                                                              <  2:204259/150‑1 (MQ=255)                                                                                      aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGcaa                                                            <  2:451998/151‑3 (MQ=255)                                                                                          ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTCGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaacaa                                                        <  1:202927/151‑1 (MQ=255)                                                                                            cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGa                                                      <  1:288263/151‑1 (MQ=255)                                                                                             gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaaca                                                         >  2:679580/1‑147 (MQ=255)                                                                                                   atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAACCAACAAGATGAGTa                                                <  1:238924/150‑1 (MQ=255)                                                                                                       tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGctct                                           >  2:598816/1‑151 (MQ=255)                                                                                                             gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATgg                                     <  1:6038/151‑1 (MQ=255)                                                                                                                 tCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTg                                  >  1:320259/1‑150 (MQ=255)                                                                                                                 tCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTg                                  <  2:438041/150‑1 (MQ=255)                                                                                                                         cATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                           <  2:294027/149‑1 (MQ=255)                                                                                                                         cATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTc                          <  1:652191/150‑1 (MQ=255)                                                                                                                             ttACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCtt                     >  1:744716/1‑151 (MQ=255)                                                                                                                               aCTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTc                    <  1:411527/150‑1 (MQ=255)                                                                                                                                  tttCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc                >  2:154725/1‑151 (MQ=255)                                                                                                                                       aCAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc           >  2:148157/1‑151 (MQ=255)                                                                                                                                                  aaaaGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTc  >  1:679107/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTC  >  NC_000913/3815661‑3815955

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