Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,251,588	A→G	D250G (GAC→GGC)	nfo →	endonuclease IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,251,588	0	A	G	100.0%	181.3 / NA	55	D250G (GAC→GGC)	nfo	endonuclease IV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/25);  total (30/25)

GGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATAACCGGGAAAAGGAGA  >  NC_000913/2251461‑2251721                                                                                                                                |                                                                                                                                      ggCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGa                                                                                                                           >  1:1050329/1‑140 (MQ=255)      tAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCa                                                                                                                      <  2:1050329/140‑1 (MQ=255)         gTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGTCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGa                                                                                                                           >  1:1016212/1‑132 (MQ=255)         gTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGTCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGa                                                                                                                           <  2:1016212/132‑1 (MQ=255)           aTCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCCTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCg                                                                                                                 <  2:1231945/140‑1 (MQ=255)             cTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCt                                                                                                               >  1:1242527/1‑140 (MQ=255)                cgcgGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGAt                                                                                                            >  2:1351030/1‑140 (MQ=255)                    ggATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTc                                                                                                        >  2:1145061/1‑140 (MQ=255)                          aCCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAAcc                                                                                                  >  1:996930/1‑140 (MQ=255)                           ccTTAACGATGCGAAAAGCGCCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAAGCa                                                                                                 <  2:915959/140‑1 (MQ=255)                              tAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCa                                                                                              <  2:996930/140‑1 (MQ=255)                              tAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCa                                                                                              <  2:572248/140‑1 (MQ=255)                               aaCGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGGTGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCaa                                                                                             >  2:907583/1‑140 (MQ=255)                               aaCGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCaa                                                                                             >  2:99967/1‑140 (MQ=255)                               aaCGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCaa                                                                                             >  2:1359065/1‑140 (MQ=255)                               aaCGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCaa                                                                                             >  2:1083905/1‑140 (MQ=255)                                  gATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAAccc                                                                                          <  1:992428/140‑1 (MQ=255)                                   aTGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                                                                                         <  2:613397/140‑1 (MQ=255)                                   aTGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCg                                                                                         <  2:1041767/140‑1 (MQ=255)                                     gCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGa                                                                                       <  1:1083905/140‑1 (MQ=255)                                      cGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGACAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGat                                                                                      >  1:848531/1‑140 (MQ=255)                                        aaaaGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGatat                                                                                    <  1:1145061/140‑1 (MQ=255)                                             cACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTggg                                                                               <  1:1439264/140‑1 (MQ=255)                                                cTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCag                                                                            <  1:907583/140‑1 (MQ=255)                                                cTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCag                                                                            <  1:99967/140‑1 (MQ=255)                                                       aGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGAt                                                                     >  1:440710/1‑140 (MQ=255)                                                          cgcgTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGc                                                                  >  2:681719/1‑140 (MQ=255)                                                          cgcgTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGc                                                                  >  1:917143/1‑140 (MQ=255)                                                            cgTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCt                                                                >  2:233672/1‑140 (MQ=255)                                                            cgTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCt                                                                >  1:493137/1‑140 (MQ=255)                                                                gACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCt                                                            <  1:155595/140‑1 (MQ=255)                                                                  ccGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGa                                                          <  1:233672/140‑1 (MQ=255)                                                                      catcatAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAgc                                                      >  2:1376939/1‑140 (MQ=255)                                                                           tAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGcaac                                                 >  2:1325751/1‑140 (MQ=255)                                                                                  ggTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGa                                          >  1:1070010/1‑140 (MQ=255)                                                                                   gTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGaa                                         >  2:1066497/1‑140 (MQ=255)                                                                                         gCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCg                                   <  2:1242527/140‑1 (MQ=255)                                                                                          cAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCgg                                  <  2:917143/140‑1 (MQ=255)                                                                                          cAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCgg                                  <  1:1325751/140‑1 (MQ=255)                                                                                            atatCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTa                                >  1:723844/1‑140 (MQ=255)                                                                                              atCGGTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGc                              <  1:1376939/140‑1 (MQ=255)                                                                                                  gTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCg                                   >  1:623559/1‑131 (MQ=255)                                                                                                  gTCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCg                                   <  2:623559/131‑1 (MQ=255)                                                                                                   tCATGATGCGTTCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGaa                         <  2:1070010/140‑1 (MQ=255)                                                                                                             ttCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATAAc               >  1:1414476/1‑140 (MQ=255)                                                                                                             ttCCGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATAAc               >  2:1209850/1‑140 (MQ=255)                                                                                                               ccGCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATAACCg             >  1:474682/1‑140 (MQ=255)                                                                                                                 gCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATAACCggg           <  2:440710/140‑1 (MQ=255)                                                                                                                 gCTGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATAACCggg           <  1:832943/140‑1 (MQ=255)                                                                                                                   tGGATCATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATAACCGGGaa         >  2:1286258/1‑140 (MQ=255)                                                                                                                        cATGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATAACCGGGAAAAGga    >  1:224077/1‑140 (MQ=255)                                                                                                                          tGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATaa                <  2:826732/126‑1 (MQ=255)                                                                                                                          tGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATaa                >  1:826732/1‑126 (MQ=255)                                                                                                                          tGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATAACCGGGAAAAGgaga  >  2:487035/1‑140 (MQ=255)                                                                                                                          tGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGCCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATaa                >  1:903010/1‑126 (MQ=255)                                                                                                                          tGCAGGGCGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGCCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATaa                <  2:903010/126‑1 (MQ=255)                                                                                                                                |                                                                                                                                      GGCTTTAAGTATCTGCGCGGGATGCACCTTAACGATGCGAAAAGCACCTTTGGCAGCCGCGTTGACCGCCATCATAGCCTCGGTGAAGGCAATATCGGTCATGATGCGTTCCGCTGGATCATGCAGGACGACCGTTTCGACGGCATTCCGCTGATCCTCGAAACCATCAACCCGGATATCTGGGCAGAAGAGATCGCCTGGCTGAAAGCGCAACAAACTGAAAAAGCGGTAGCCTGAAGATGAATAACCGGGAAAAGGAGA  >  NC_000913/2251461‑2251721

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