Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,521	A→G	V83A (GTG→GCG)	pyrE ←	orotate phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,521	0	A	G	100.0%	163.4 / NA	52	V83A (GTG→GCG)	pyrE	orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/29);  total (23/29)

GCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGC  >  NC_000913/3815374‑3815660                                                                                                                                                    |                                                                                                                                            gCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAg                                                                                                                                          <  2:208320/151‑1 (MQ=255)    gTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCgg                                                                                                                                       <  2:566974/151‑1 (MQ=255)    gTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCCCCCGTGCCGCAGCgg                                                                                                                                       >  1:678173/1‑151 (MQ=255)    gTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGGCCCGGGCGGGATGCCCCGCCCGTGCCGCCGCgg                                                                                                                                       >  2:124986/1‑151 (MQ=255)       ccGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGtt                                                                                                                                     <  1:490613/150‑1 (MQ=255)          gcggTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTgtgg                                                                                                                                 <  2:432396/151‑1 (MQ=255)            ggTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTgtggtg                                                                                                                               >  1:593291/1‑151 (MQ=255)                  cacaTCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATc                                                                                                                         <  1:259595/151‑1 (MQ=255)                       catcTACCAGCATTACGCGTCCTTGTAGCGCGCTACCAGCCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGa                                                                                                                     <  1:127228/150‑1 (MQ=255)                            aCCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATc                                                                                                                   <  2:240981/147‑1 (MQ=255)                            aCCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCt                                                                                                                <  2:423378/150‑1 (MQ=255)                             ccAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCttt                                                                                                              <  2:342239/151‑1 (MQ=255)                                       gcgTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGGGGCCATCGGGATCCCTTTGTAAGCAggg                                                                                                    >  1:435187/1‑151 (MQ=255)                                          tCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGGTCCCTTTGTAAGCAGGGcc                                                                                                  >  1:378977/1‑150 (MQ=255)                                              tGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCaaa                                                                                               <  1:357160/149‑1 (MQ=255)                                               gTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCCCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCaaa                                                                                               >  2:633392/1‑148 (MQ=255)                                               gTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGGCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGGCCAGGGCGTGAGGGTACGCCCCTGCCGCAGCGGGTGGGGTGGCCATCGGGATCCCTTTGTAAGCCGGGCCaa                                                                                                >  2:460164/1‑147 (MQ=255)                                                     cgcTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGGCGTGATGCTCCGCCAGTGCCGCAGCGGTTGGGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAAcagc                                                                                           >  1:121237/1‑146 (MQ=255)                                                             cccaGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCTGGTCGTGATGCTCCGCAATTGCCGCAGCGGTTGGGGTGGCAATCGGGATCCCTTTGTAAACAGGGCCCAACAGCAGATCGAAc                                                                                  >  2:526662/2‑147 (MQ=255)                                                              ccaGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCaa                                                                              >  1:229579/1‑150 (MQ=255)                                                               caGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATg                                                                            <  1:8194/151‑1 (MQ=255)                                                                 gATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGcc                                                                          >  1:9516/1‑151 (MQ=255)                                                                     gccgccTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGTAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGaa                                                                      <  1:372130/151‑1 (MQ=255)                                                                     gccgccTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGaa                                                                      <  2:284062/151‑1 (MQ=255)                                                                        gccTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGGCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATcc                                                                   >  2:367510/1‑151 (MQ=255)                                                                                        tttGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAg                                                      <  1:9396/148‑1 (MQ=255)                                                                                        tttGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCg                                                    >  1:592482/1‑150 (MQ=255)                                                                                             ttcttTGCGGTTAAAGCAGTACGGCAGGTCCAGGGCGTGATGCTCCTCCCGTGCCGCAGCGGTTGTGGTGGCCATCGGGGTCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAAACCAACAACGCTTCAGCGGaaaa                                               >  2:213355/1‑150 (MQ=255)                                                                                                 ttGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACgg                                           <  1:321172/150‑1 (MQ=255)                                                                                                     ggTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCCCCAACGCTTCAGCGTAAAAACGGCCTa                                       >  2:377356/1‑150 (MQ=255)                                                                                                      gTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAAc                                     <  1:441161/151‑1 (MQ=255)                                                                                                      gTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGGGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAAATCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTaa                                      >  2:261758/1‑150 (MQ=255)                                                                                                         aaaGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAg                                   >  1:191313/1‑150 (MQ=255)                                                                                                         aaaGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGt                                  <  2:671001/151‑1 (MQ=255)                                                                                                         aaaGCAGTACGGCAGGTCCAGGGCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGGGGGAAACGGGATCCCTTTGTAAGAAGGGCCAAAAAGCAGAACGAACTCAATGCCGGAAACCACCAACGCTTCAGCGGAAAAAAGGCCTAACa                                    >  2:349106/1‑149 (MQ=255)                                                                                                          aaGCAGTACGGCAGGTCCCGGTCGCGATGCCCCCCCAGTGCCGCACCGGTTGTGGGGGCCATCGGGATCCCTTTGTAAGCAGGGCCCAACAGCCGATCGCACTCAATGCCGGAATCCACCCACGCCTCCGCGTAAAAACGGccccac                                     >  2:364976/1‑143 (MQ=255)                                                                                                                 aCGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGa                            <  2:627319/149‑1 (MQ=255)                                                                                                                 aCGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGAt                           <  2:371961/150‑1 (MQ=255)                                                                                                                 aCGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGAt                           <  1:513587/150‑1 (MQ=255)                                                                                                                  cGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATcg                         <  2:642557/151‑1 (MQ=255)                                                                                                                    gCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATcgcg                       <  2:457006/151‑1 (MQ=255)                                                                                                                      aggtccaggtcGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGcc                     >  1:637260/1‑151 (MQ=255)                                                                                                                        gtccaggtcGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCg                   <  1:556036/151‑1 (MQ=255)                                                                                                                          ccaggtcGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATTCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGt                 >  1:366910/1‑151 (MQ=255)                                                                                                                          ccaggtcGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGt                 >  1:182104/1‑151 (MQ=255)                                                                                                                           caggtcGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTa                <  2:623899/151‑1 (MQ=255)                                                                                                                             ggtcGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTAt               <  2:291699/150‑1 (MQ=255)                                                                                                                              gtcGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTa             <  1:119047/151‑1 (MQ=255)                                                                                                                               tcGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTa             >  2:487759/1‑150 (MQ=255)                                                                                                                                cGTGATGCTCCGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTaaa           <  1:100751/151‑1 (MQ=255)                                                                                                                                         ccGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCgg   <  2:657198/150‑1 (MQ=255)                                                                                                                                         ccGCCAGTGCCGCAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGc  <  1:669204/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                            GCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAGGGCCAAACAGCAGATCGAACTCAATGCCGGAATCCACCAACGCTTCAGCGTAAAAACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGC  >  NC_000913/3815374‑3815660

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