Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,325,822	G→A	intergenic (‑81/+577)	yjdN ← / ← yjdM	conserved protein YjdN/conserved protein YjdM

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,325,822	0	G	A	100.0%	109.6 / NA	32	intergenic (‑81/+577)	yjdN/yjdM	conserved protein YjdN/conserved protein YjdM
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/17);  total (15/17)

AGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGC  >  NC_000913/4325687‑4325919                                                                                                                                        |                                                                                                  aGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCg                                                                                               <  1:23805/140‑1 (MQ=255)     aaTCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGc                                                                                           >  2:729931/1‑140 (MQ=255)           gTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAgg                                                                                     <  1:565798/140‑1 (MQ=255)             cGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCa                                                                                   <  1:729931/140‑1 (MQ=255)                 aCAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCg                                                                               <  2:486074/140‑1 (MQ=255)                     ttACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCAc                                                                           >  1:356474/1‑140 (MQ=255)                     ttACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCAc                                                                           >  2:679204/1‑140 (MQ=255)                          ggCAAAAGAGAGGTACGGACTTAACGGCATGATGGTGACCTCAGCTAAGAGAAGCCAGTTAAGTGTAGTTCGCAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCAccattg                                                                      >  2:685256/1‑135 (MQ=38)                           gCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGc                                                                     <  2:501208/140‑1 (MQ=255)                                   gagGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATgcg                                                             <  2:356474/140‑1 (MQ=255)                                      gTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATgcggcg                                                          >  1:151553/1‑140 (MQ=255)                                       tAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCC‑GTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTa                                                        <  1:824800/140‑1 (MQ=25)                                                aaCGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCtt                                                >  1:442270/1‑140 (MQ=255)                                                  cGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTAt                                              <  2:656736/140‑1 (MQ=255)                                                    gCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCC‑TATCCg                                           >  2:800911/1‑140 (MQ=25)                                                        gatgGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                        >  2:707456/1‑140 (MQ=255)                                                                 cTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAggg                               >  2:876858/1‑140 (MQ=255)                                                                    aGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGt                              >  1:532817/1‑138 (MQ=255)                                                                    aGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGt                              <  2:532817/138‑1 (MQ=255)                                                                    aGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGc                            >  2:612660/1‑140 (MQ=255)                                                                        aaGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                        <  1:876858/140‑1 (MQ=255)                                                                        aaGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                        <  1:864869/140‑1 (MQ=255)                                                                        aaGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                        <  2:442270/140‑1 (MQ=255)                                                                           agaAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAtt                     >  1:231935/1‑140 (MQ=255)                                                                             aaGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCTCCATTGCAGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTg                   <  1:685256/140‑1 (MQ=18)                                                                             aaGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTg                   >  2:691898/1‑140 (MQ=255)                                                                               gCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTa                 <  1:707456/140‑1 (MQ=255)                                                                                  aGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGTTGCGGGAATTTGTAGGc              <  2:151553/140‑1 (MQ=255)                                                                                  aGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGc              <  2:115678/140‑1 (MQ=255)                                                                                  aGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCATGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGc              <  2:694434/140‑1 (MQ=255)                                                                                              ttCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcgc  >  1:97541/1‑140 (MQ=255)                                                                                              ttCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcgc  >  1:83156/1‑140 (MQ=255)                                                                                              ttCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcgc  >  1:734850/1‑140 (MQ=255)                                                                                                                                        |                                                                                                  AGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGC  >  NC_000913/4325687‑4325919

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: ‑