Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,390,897	A→G	W13R (TGG→CGG)	ymjD ←	uncharacterized protein YmjD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,390,897	0	A	G	100.0%	149.0 / NA	43	W13R (TGG→CGG)	ymjD	uncharacterized protein YmjD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/26);  total (17/26)

AGTTTGTTTAATGGTTTGCTTATCTGCCAGTTGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCAGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGC  >  NC_000913/1390773‑1391021                                                                                                                             |                                                                                                                             aGTTTGTTTAATGGTTTGCTTATCTGCCAGTTGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCTAGATTATTGGCAGCAATATGCCGGTCCatat                                                                                                                      <  2:21110/133‑1 (MQ=255) aGTTTGTTTAATGGTTTGCTTATCTGCCAGTTGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCatat                                                                                                                      >  1:408897/1‑133 (MQ=255) aGTTTGTTTAATGGTTTGCTTATCTGCCAGTTGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCatat                                                                                                                      <  1:94339/133‑1 (MQ=255) aGTTTGTTTAATGGTTTGCTTATCTGCCAGTTGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCatat                                                                                                                      >  2:496218/1‑133 (MQ=255) aGTTTGTTGAATGGTCTGCTTATCTGCCAGCTGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCatat                                                                                                                      <  1:207293/133‑1 (MQ=255)  gtttgtttAATGGTTTGCTTATCTGCCAGTTGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATc                                                                                                                     <  1:282133/133‑1 (MQ=255)     tgtttAATGGTTTGCTTATCTGCCAGTTGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACt                                                                                                                  >  1:80154/1‑133 (MQ=255)       tttAATGGTTTGCTTATCTGCCAGTTGATTAATCACATGGCGGGCTATCTGGCCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTAt                                                                                                                >  1:468273/1‑133 (MQ=255)         tAATGGTTTGCTTATATGCCAGTTGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTg                                                                                                              <  1:168107/133‑1 (MQ=255)           aTGGTTTGCTTATCTGCCATTTTATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCg                                                                                                            <  2:510659/133‑1 (MQ=255)                 tGCTTATCTGCCAGTTGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAAGATGACGGTCCATATAACTATTGCGGATTTg                                                                                                      >  1:168616/1‑133 (MQ=255)                     tATCTGCCAGTTGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAata                                                                                                  <  1:496218/133‑1 (MQ=255)                       tCTGCCAGTTGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATg                                                                                                <  1:354298/133‑1 (MQ=255)                          gCCAGTTGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGttt                                                                                             >  1:74410/1‑133 (MQ=255)                              gTTGATTAATCACATGGCGGGCTATCTGGCCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCAtt                                                                                         <  2:468273/133‑1 (MQ=255)                              gTTGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCAtt                                                                                         >  1:259522/1‑133 (MQ=255)                                tGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTAc                                                                                       <  2:80154/133‑1 (MQ=255)                                         acaTGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTcc                                                                              <  2:408897/133‑1 (MQ=255)                                          aaTGGCGGCCTATCGGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCt                                                                             <  1:433572/132‑1 (MQ=255)                                           aTGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATTTGTATCTTTACAGAGTTTCCTc                                                                            >  1:145051/1‑133 (MQ=255)                                               cGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTATCACTAAAGAGTTTCCTCAAAg                                                                        >  1:228951/1‑133 (MQ=255)                                                  gCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAAc                                                                     <  2:259522/133‑1 (MQ=255)                                                      tCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAg                                                                 <  1:341234/133‑1 (MQ=255)                                                         gACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTa                                                              >  1:420701/1‑133 (MQ=255)                                                             gCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTCCAGAGTTACCTTAAAGAACGGAGCTAATAc                                                          >  1:251102/1‑133 (MQ=255)                                                                aGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATATTTCATTAGAGATTTTCCTCAAAGATCGGATCTAATACTGa                                                       >  2:221993/1‑133 (MQ=255)                                                                  cACCAAGAATAAGTACATTTTTCATGAAAATTTTCCCAGATGATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGaaa                                                     <  2:228951/133‑1 (MQ=255)                                                                        gAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCtt                                               <  1:346049/133‑1 (MQ=255)                                                                          aTCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCa                                             <  2:76968/133‑1 (MQ=255)                                                                             aGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGcc                                          >  2:7220/1‑133 (MQ=255)                                                                               tACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATAATGTAATTGCTTCAGCcac                                        >  1:199738/1‑133 (MQ=255)                                                                                      ttCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCt                                 >  1:397857/1‑133 (MQ=255)                                                                                          tGAAAATTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGc                             <  2:391514/133‑1 (MQ=255)                                                                                           gAAAACTATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTACTGCTGATTTGAATATGTTTCTATACAGAGTTTTCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGcc                            <  2:91009/133‑1 (MQ=255)                                                                                                 tATCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATa                      <  1:495239/133‑1 (MQ=255)                                                                                                   tCCCAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATaaa                    <  2:420701/133‑1 (MQ=255)                                                                                                   tCCCAGATCATTGGCTGCTATAAGCCGGTCTATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATaaa                    <  1:221993/133‑1 (MQ=255)                                                                                                     ccAGATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGa                  <  1:288265/133‑1 (MQ=255)                                                                                                        gATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATaaa                    >  2:374373/1‑128 (MQ=255)                                                                                                        gATTATTGGCAGCAATATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATaaa                    <  1:374373/128‑1 (MQ=255)                                                                                                         attattGGCAACAATATGACGGTCTATAGCACAATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTACTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCa              <  2:168616/133‑1 (MQ=255)                                                                                                          ttattGGCGGCAATATGCTGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCat             <  1:193465/133‑1 (MQ=255)                                                                                                                     aaTATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTgcgc  <  1:7220/133‑1 (MQ=255)                                                                                                                     aaTATGCCGGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTgcgc  >  1:116760/1‑133 (MQ=255)                                                                                                                             |                                                                                                                             AGTTTGTTTAATGGTTTGCTTATCTGCCAGTTGATTAATCACATGGCGGGCTATCTGACCGCCAGCACCAAGAATCAGTACATTTTTCATGAAAATTATCCCAGATTATTGGCAGCAATATGCCAGTCCATATCACTATTGCGGATTTGAATATGTTTCATTACAGAGTTTCCTCAAAGAACGGAGCTAATACTGAAATTGCTTCAGCCACATAGGCTTTGCCGTCATAAAGATCCATATGGTTTGCGC  >  NC_000913/1390773‑1391021

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