Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,878,603	T→C	D84G (GAT→GGT)	ygbF ←	CRISPR adaptation ssRNA endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,878,603	0	T	C	100.0%	158.3 / NA	46	D84G (GAT→GGT)	ygbF	CRISPR adaptation ssRNA endonuclease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/19);  total (27/19)

TATACATTTGTTCTTAAAGCATTTTTTCCCATAAAAACAACCCACCAACCTTAATGTAACATTTCCTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCATCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCATCACTACATTGCCTTCTTCCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGA  >  NC_000913/2878472‑2878737                                                                                                                                    |                                                                                                                                       tataCATTTGTTCTTAAAGCATTTTTTCCCATAAAAACAACCCACCAACCTTAATGTAACATTTCCTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATc                                                                                                                                 >  1:1027616/1‑139 (MQ=255)           ttCTTAAAGCATTTTTTCCCATAAAAACAACCCACCAACCTTAATGTAACATTTCCTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTa                                                                                                                       >  2:671346/1‑139 (MQ=255)           ttCTTAAAGCATTTTTTCCCATAAAAACAACCCACCAACCTTAATGTAACATTTCCTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTa                                                                                                                       >  2:682998/1‑139 (MQ=255)                           tCCCATAAAAACAACCCACCAACCTTAATGTAACATTTCCTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCaa                                                                                                       <  2:1027616/139‑1 (MQ=255)                                   aaaCAACCCACCAACCTTAATGTAACATTTCCTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTgg                                                                                               <  1:682998/139‑1 (MQ=255)                                         cccACCAACCTTAATGTAACATTTCCTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCa                                                                                         <  1:438882/139‑1 (MQ=255)                                         cccACCAACCTTAATGTAACATTTCCTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCa                                                                                         <  2:805861/139‑1 (MQ=255)                                         cccACCAACCTTAATGTAACATTTCCTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCa                                                                                         >  1:805861/1‑139 (MQ=255)                                                           aCATTTCCTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTc                                                                       >  2:380460/1‑139 (MQ=255)                                                                tCCTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAACTCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGc                                                                  >  2:259897/1‑139 (MQ=255)                                                                tCCTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCgggg                                                                                                                         <  1:491756/84‑1 (MQ=255)                                                                tCCTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCgggg                                                                                                                         >  2:491756/1‑84 (MQ=255)                                                                 ccTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTgg                                                                                               <  2:733119/109‑1 (MQ=255)                                                                 ccTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTgg                                                                                               >  1:733119/1‑109 (MQ=255)                                                                 ccTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGcc                                                                 >  1:983646/1‑139 (MQ=255)                                                                 ccTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGcc                                                                 >  1:693987/1‑139 (MQ=255)                                                                 ccTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGcc                                                                 >  1:486361/1‑139 (MQ=255)                                                                     attattaAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATg                                                             >  2:417284/1‑139 (MQ=255)                                                                      ttattaAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGc                                                            <  1:380460/139‑1 (MQ=255)                                                                           aaaGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCATCa                                                       >  1:498974/1‑139 (MQ=255)                                                                             aGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCaa                                                                                                       <  1:836314/89‑1 (MQ=255)                                                                             aGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCaa                                                                                                       >  2:836314/1‑89 (MQ=255)                                                                                tCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCgggg                                                                                                                         <  2:1019322/68‑1 (MQ=255)                                                                                tCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCgggg                                                                                                                         >  1:1019322/1‑68 (MQ=255)                                                                                tCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCATCa                                                       <  2:941292/134‑1 (MQ=255)                                                                                tCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCATCa                                                       >  1:941292/1‑134 (MQ=255)                                                                                 cAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAATTCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGc                                                            <  1:1058657/128‑1 (MQ=255)                                                                                 cAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTAc                                                                                                                              <  1:957075/62‑1 (MQ=255)                                                                                 cAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTAc                                                                                                                              >  2:957075/1‑62 (MQ=255)                                                                                 cAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCATCACTACAt                                                 >  1:82236/1‑139 (MQ=255)                                                                                 cAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTAACGGGGTACGCCTGTTTAACCCAACTGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGc                                                            >  2:1058657/1‑128 (MQ=255)                                                                                         tCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTAc                                                                                                                              >  1:847918/1‑54 (MQ=255)                                                                                         tCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTAc                                                                                                                              <  2:847918/54‑1 (MQ=255)                                                                                            ttGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCATCACTACATTGCCTTCTTcc                                      >  1:360195/1‑139 (MQ=255)                                                                                             tGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCATCACTACATTGCCTTCTTCCg                                     >  2:878065/1‑139 (MQ=255)                                                                                             tGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCATCACTACATTGCCTTCTTCCg                                     >  2:259238/1‑139 (MQ=255)                                                                                              gTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGtt                                                                                <  2:225722/95‑1 (MQ=255)                                                                                              gTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGtt                                                                                >  1:225722/1‑95 (MQ=255)                                                                                                           tAAAAAAGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCAGCACTACATTGCCTTCTTCCGCCAGTCCAGCTAtt                       >  1:689591/1‑139 (MQ=255)                                                                                                            aaaaaaGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCATCACTACATTGCCTTCTTCCGCCAGTCCAGCTAttt                      <  2:82236/139‑1 (MQ=255)                                                                                                              aaaaGACACCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCATCACTACATTGCCTTCTTCCGCCAGTCCAGCTATTTGt                    <  2:486361/139‑1 (MQ=255)                                                                                                                    cacCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCg                                                                      >  1:951401/1‑83 (MQ=255)                                                                                                                    cacCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCg                                                                      <  2:951401/83‑1 (MQ=255)                                                                                                                    cacCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCATCACTACATTGCCTTCTTCCGCCAGTCCAGCTATTTGTTc                  >  1:354321/1‑135 (MQ=255)                                                                                                                    cacCAACCTTAAACCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCATCACTACATTGCCTTCTTCCGCCAGTCCAGCTATTTGTTc                  <  2:354321/135‑1 (MQ=255)                                                                                                                               aaCCACCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCATCACTACATTGCCTTCTTCCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACg   <  2:693987/139‑1 (MQ=255)                                                                                                                                accaccCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCATCACTACATTGCCTTCTTCCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGa  >  2:410111/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                       TATACATTTGTTCTTAAAGCATTTTTTCCCATAAAAACAACCCACCAACCTTAATGTAACATTTCCTTATTATTAAAGATCAGCTAATTCTTTGTTTTCAAACAGGTAAAAAAGACACCAACCTTAAACCATCCAAATCTACCGGGGTACGCCTGTTTAACCCAAATGTCTGGAACTCAAATCCCGTTTCCGTATTCGTTGCCCATGCCATCACTACATTGCCTTCTTCCGCCAGTCCAGCTATTTGTTCCCAGATCATTTCACGA  >  NC_000913/2878472‑2878737

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