Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	3,815,825	C→G	intergenic (‑57/+9)	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_bme_pgi	3,815,825	0	C	G	97.8%	142.6 / ‑3.9	45	intergenic (‑57/+9)	pyrE/rph	orotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (23/21);  minor base . (1/0);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGC  >  NC_000913_3_bme_pgi/3815683‑3815965                                                                                                                                                |                                                                                                                                             cgcCCGGATTGCAGCGTAAACTCGCAAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTaaa                                                                                                                                        <  1:477529/149‑1 (MQ=255)      ggATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCa                                                                                                                                   <  2:637233/149‑1 (MQ=255)      ggATTTCAGCGTAAACTCCCCAAACTTTAACACCTGCTTGCCAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCa                                                                                                                                   <  1:428026/149‑1 (MQ=255)        aTTTCAGCGGAAATCCGCCAAACTTTAACACCTGCTTGCTAGGCTCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGt                                                                                                                                 <  2:90406/149‑1 (MQ=255)              gCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCa                                                                                                                           <  1:436400/149‑1 (MQ=255)                 tAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTAAAGACAAAAAAAAGGCGACCCACTCATTCG‑CTTAAAAATCAGTTTGCCAgcg                                                                                                                        >  2:509818/1‑149 (MQ=255)                     cTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTAA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGcgccgc                                                                                                                    >  1:275423/1‑149 (MQ=255)                          cAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCt                                                                                                               >  2:46761/1‑149 (MQ=255)                          cAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCt                                                                                                               <  1:310782/149‑1 (MQ=255)                          cAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCAACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCt                                                                                                               <  2:640756/149‑1 (MQ=255)                                   acacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTa                                                                                                      >  1:659167/1‑149 (MQ=255)                                   acacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTa                                                                                                      >  1:513158/1‑149 (MQ=255)                                       ctgctgGTTAAGGGTAAATTCAATAAACCGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAAt                                                                                                  <  2:52945/149‑1 (MQ=255)                                       cTGCTTGCCAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAAt                                                                                                  <  1:123934/149‑1 (MQ=255)                                         gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATgg                                                                                                >  1:589144/1‑149 (MQ=255)                                          cttgctAAGCGCAAATTCAATAAACTGGCGCTGATCTGAATTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGga                                                                                               >  1:571105/1‑149 (MQ=255)                                            tgctAAGCGCAAATTCAATAAACTGGCGCGGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAATTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                                                                                             <  1:740924/149‑1 (MQ=255)                                             gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc                                                                                            <  2:589144/149‑1 (MQ=255)                                              ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcg                                                                                           >  1:721917/1‑149 (MQ=255)                                                   cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTcc                                                                                      <  1:61784/149‑1 (MQ=255)                                                       aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTc                                                                                  >  1:44242/1‑149 (MQ=255)                                                       aaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTc                                                                                  >  1:42622/1‑149 (MQ=255)                                                        aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCg                                                                                 <  1:46761/149‑1 (MQ=255)                                                        aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCg                                                                                 <  2:65424/149‑1 (MQ=255)                                                                 aCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCCTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGc                                                                        >  2:58324/1‑149 (MQ=255)                                                                       gcTGATATGGTTTCATGCCTTCGCTCCTCCCGTTCTGTTTTGTCAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaacaa                                                                  <  1:403971/149‑1 (MQ=255)                                                                        cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaacaaa                                                                 >  2:143798/1‑148 (MQ=255)                                                                        cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCAAGGGCCAACAAg                                                                 >  2:141320/1‑149 (MQ=255)                                                                            tatTGTTTCATGCCTTCGCTCCCCAGCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGa                                                             <  1:97523/149‑1 (MQ=255)                                                                                gtgtCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                         >  2:204182/4‑149 (MQ=255)                                                                                gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TGAGCCGGCTTAAAAATCAGTTTGCCGGCCCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGTGCCAGCGCCAACAAGTTAAGTaa                                                         >  2:9513/1‑148 (MQ=255)                                                                                     aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCtt                                                    <  2:42622/149‑1 (MQ=255)                                                                                     aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCtt                                                    <  2:44242/149‑1 (MQ=255)                                                                                        acttggatcctcATCTTACTTTTCTACAGACAAAAAAAAGGCGAGTCA‑TCAGTCGGCTTAAAAATCAGTTTTCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCATCTCGGGCTTTAGCCTACAAGACGAGTAGCTCTTCAt                                                 <  2:745392/142‑1 (MQ=255)                                                                                         cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGCTTCCCCTCGGGCCAGAGCCAACAACATGAGTAGCTCTTCATg                                                >  1:289390/1‑149 (MQ=255)                                                                                           tCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGGCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATggg                                              <  1:70352/149‑1 (MQ=255)                                                                                                    aTCTTACTTTTCTACAGACAAAACAAAGGCGACTCA‑TCAATCGGCTCAAAAATCATTTCGCCAGCGCCCCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCCCTTCATCGCTGAACGGCt                                     >  2:371679/1‑149 (MQ=255)                                                                                                    aTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                                     >  1:42740/1‑149 (MQ=255)                                                                                                       ttACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCACCAAGATGAGTAGCTCTTCATGAGTGAACGGCTCGc                                  >  1:749798/1‑149 (MQ=255)                                                                                                                tACAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCg                         >  1:761364/1‑149 (MQ=255)                                                                                                                  cAGACAAAAAAAAGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTc                       >  1:206587/1‑149 (MQ=255)                                                                                                                         aaaaaaGGCGACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCAGCCTTCTGCGTCGCTACAATAGATTCGATTCCCCTCGCGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCa                >  1:695251/1‑149 (MQ=255)                                                                                                                                  gACTCA‑TCAGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAataa       >  1:230947/1‑147 (MQ=255)                                                                                                                                   aCTCA‑TCAGTCGGCTTAAAAATCAGGTTCCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgat      <  1:86497/149‑1 (MQ=255)                                                                                                                                       a‑tcaGTCGGCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTAAATGATgcgc  >  1:240736/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                             CGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCA‑TCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGC  >  NC_000913_3_bme_pgi/3815683‑3815965

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑