Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,471,253	A→G	F65F (TTT→TTC)	pyrI ←	aspartate carbamoyltransferase, PyrI subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,471,253	0	A	G	100.0%	148.5 / NA	43	F65F (TTT→TTC)	pyrI	aspartate carbamoyltransferase, PyrI subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/21);  total (22/21)

GACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGAACAGACTCAACAGCTTAAAACCGATCTGGGCGGGGATA  >  NC_000913/4471124‑4471388                                                                                                                                  |                                                                                                                                        gACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTAtttt                                                                                                                               <  2:348182/140‑1 (MQ=255)            tGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATc                                                                                                                    >  2:653521/1‑140 (MQ=255)                gATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGAt                                                                                                                <  2:484407/140‑1 (MQ=255)                 aTGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATc                                                                                                               >  1:819165/1‑140 (MQ=255)                   gcgcTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCtt                                                                                                             <  1:754797/140‑1 (MQ=255)                    cgcTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCttt                                                                                                            <  1:803955/140‑1 (MQ=255)                        ccGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCgg                                                                                                        <  2:819165/140‑1 (MQ=255)                        ccGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCgg                                                                                                        >  1:289030/1‑140 (MQ=255)                            cAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCa                                                                                                    >  2:244021/1‑140 (MQ=255)                               aCTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCAtct                                                                                                 <  1:357283/140‑1 (MQ=255)                                 ttGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCg                                                                                               >  1:460345/1‑140 (MQ=255)                                   gggCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGcc                                                                                             >  2:765991/1‑140 (MQ=255)                                     gcgcGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAg                                                                                           <  2:186911/140‑1 (MQ=255)                                         gATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAgg                                                                                       >  2:189612/1‑140 (MQ=255)                                          aTTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGc                                                                                      <  2:289030/140‑1 (MQ=255)                                          aTTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGc                                                                                      <  1:623853/140‑1 (MQ=255)                                          aTTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGc                                                                                      <  1:244021/140‑1 (MQ=255)                                              aCCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGTCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGt                                                                                  >  1:529743/1‑140 (MQ=255)                                                   ccacTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGtt                                                                                 >  1:463705/1‑136 (MQ=255)                                                   ccacTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGtt                                                                                 <  2:463705/136‑1 (MQ=255)                                                      cTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCa                                                                          >  2:796957/1‑140 (MQ=255)                                                       ttCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCaa                                                                         >  2:815009/1‑140 (MQ=255)                                                         cATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATg                                                                       >  1:798062/1‑140 (MQ=255)                                                                gTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTACGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATgc                                                                >  2:769632/1‑140 (MQ=255)                                                                   gATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCt                                                             >  2:80380/1‑140 (MQ=255)                                                                   gATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCt                                                             >  2:252468/1‑140 (MQ=255)                                                                          ttAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGt                                                      <  1:653521/140‑1 (MQ=255)                                                                                gTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGt                                                <  1:612644/140‑1 (MQ=255)                                                                                        cggcgCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGa                                        <  1:774171/140‑1 (MQ=255)                                                                                        cggcgCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGa                                        >  1:551624/1‑140 (MQ=255)                                                                                             cATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCgg                                                 >  2:438408/1‑126 (MQ=255)                                                                                             cATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCgg                                                 <  1:438408/126‑1 (MQ=255)                                                                                             cATACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGAACAGa                                   >  2:460689/1‑140 (MQ=255)                                                                                               tACAATGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGAACAGAcc                                 <  2:793733/140‑2 (MQ=255)                                                                                                   aTGCCAGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGAACAGACTCAAc                             >  1:110376/1‑140 (MQ=255)                                                                                                        aGTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGAACAGACTCAACAGCtt                        <  1:252468/140‑1 (MQ=255)                                                                                                         gTTGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGAACAGACTCAACAGCTTa                       >  1:370310/1‑140 (MQ=255)                                                                                                          ttGATCTACTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGAACAGACTCAACAGCTTaa                      >  1:693713/1‑140 (MQ=255)                                                                                                                  cTTGATCTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGAACAGACTCAACAGCTTAAAACCGATc              <  1:57435/140‑1 (MQ=255)                                                                                                                        cTTCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGAACGGACTCAACAGCTTAAAACCGATCTGGGCg        <  1:420042/140‑1 (MQ=255)                                                                                                                          tCACTCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGAACAGACTCAACAGCTTAAAACCGATCTGGGCggg      >  1:11770/1‑140 (MQ=255)                                                                                                                              tCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGAACAGACTCAACAGCTTAAAACCGATCTGGGCGGGGata  <  2:370310/140‑1 (MQ=255)                                                                                                                              tCAAGAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGAACAGACTCAACAGCTTAAAACCGATCTGGGCGGGGata  <  1:460689/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                        GACCAGCACATTGTCGATGCGCTCCGGCAGACTTGGGCGCGATTTACCCACCACTTCATAGTTGTCGATACGGTTAACCGTGGCTTGCGGCGCATACAATGCCAGTTGATCTACTTGATCTTCACTCAAAAAGGTATTTTCGATTTTGATCAGATCTTTGCGGCCCATCTCGCCAGAAGGCAGGTTCAGACCAATGGTGATGCGCTGATCCGTTTCGGTCAGCTTGAACAGACTCAACAGCTTAAAACCGATCTGGGCGGGGATA  >  NC_000913/4471124‑4471388

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