Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,815,808	(A)8→7	intergenic (‑40/+55)	pyrE ← / ← rph	orotate phosphoribosyltransferase/truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,815,801	0	A	.	100.0%	174.2 / NA	41	intergenic (‑33/+62)	pyrE/rph	orotate phosphoribosyltransferase/truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (16/25);  total (16/25)

GCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC  >  NC_000913/3815723‑3815875                                                                               |                                                                           gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaa                                                                           >  1:683654/1‑80 (MQ=255)  cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaa                                                                          <  2:683598/80‑1 (MQ=255)    tgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaa                                                                        >  1:683655/1‑80 (MQ=255)     gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaa                                                                        <  2:683599/79‑1 (MQ=255)         aGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaaaaggc                                                                   >  2:683662/1‑77 (MQ=255)         aGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaaaaggc                                                                   >  2:683641/1‑77 (MQ=255)          gcgcAAATTCAATAAACTGGCGCTGATATGGATTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGa                                                                <  2:683600/81‑1 (MQ=255)            gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACt                                                              >  2:683675/1‑81 (MQ=255)            gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACt                                                              >  2:683646/1‑81 (MQ=255)            gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACt                                                              >  1:683658/1‑81 (MQ=255)               aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACtcat                                                           >  1:683640/1‑81 (MQ=255)                aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACtcatc                                                          <  2:683601/81‑1 (MQ=255)                   cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGt                                                       >  2:683679/1‑81 (MQ=255)                    aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTc                                                      <  1:683603/81‑1 (MQ=255)                     aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCg                                                     <  1:683604/81‑1 (MQ=255)                        aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCt                                                  <  1:683606/81‑1 (MQ=255)                        aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCt                                                  <  2:683604/81‑1 (MQ=255)                         aCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCtt                                                 <  1:683607/81‑1 (MQ=255)                          cTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTa                                                <  1:683608/81‑1 (MQ=255)                           tGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTa                                                >  1:683652/1‑80 (MQ=255)                             gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTaaa                                              <  2:683609/80‑1 (MQ=255)                             gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAGGCGACTCATCAGTCGCCTTaaaa                                             <  1:683609/80‑1 (MQ=255)                                      tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTg                                    <  2:683610/81‑1 (MQ=255)                                        gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGc                                   <  2:683611/80‑1 (MQ=255)                                           tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgc                               >  2:683670/1‑81 (MQ=255)                                            cATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgc                               <  2:683612/80‑1 (MQ=255)                                             aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgcgc                             <  1:683614/81‑1 (MQ=255)                                                 cTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCgccgcc                         >  2:683665/1‑81 (MQ=255)                                                 cTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCgccgcc                         >  2:683688/1‑81 (MQ=255)                                                   tCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCtt                       >  2:683656/1‑81 (MQ=255)                                                    cGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCtt                       <  2:683614/80‑1 (MQ=255)                                                    cGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCt                        <  1:683616/79‑1 (MQ=255)                                                    cGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTc                      >  2:683678/1‑81 (MQ=255)                                                     gCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCt                     <  1:683617/81‑1 (MQ=255)                                                        cctcATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCg                  >  1:683710/1‑81 (MQ=255)                                                        cctcATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCg                  >  2:683655/1‑81 (MQ=255)                                                            aTCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTc                <  2:683618/79‑1 (MQ=255)                                                            aTCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCg               <  1:683618/80‑1 (MQ=255)                                                              cTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTa            >  1:683670/1‑81 (MQ=255)                                                                 aCTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAc           <  2:683619/79‑1 (MQ=255)                                                                 aCTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAc           <  1:683621/79‑1 (MQ=255)                                                                 aCTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACaa         >  2:683684/1‑81 (MQ=255)                                                                 aCTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACa          <  2:683621/80‑1 (MQ=255)                                                                  cTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAAt        <  1:683623/81‑1 (MQ=255)                                                                         aCAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc  >  2:683667/1‑80 (MQ=255)                                                                          cAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc  <  2:683624/79‑1 (MQ=255)                                                                          cAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc  <  1:683624/79‑1 (MQ=255)                                                                               |                                                                           GCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTC  >  NC_000913/3815723‑3815875

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