Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,154,549	T→C	E100G (GAG→GGG)	argE ←	acetylornithine deacetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,154,549	0	T	C	100.0%	145.5 / NA	40	E100G (GAG→GGG)	argE	acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/17);  total (23/17)

GCCATACTGGTTTCTTCATCAGCAGTCGCCAGAATGTAGAGCGGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGC  >  NC_000913/4154402‑4154689                                                                                                                                                    |                                                                                                                                             gCCATACTGGTTTCTTCATCAGCAGTCGCCAGAATGTAGAGCGGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCg                                                                                                                                           >  1:393111/1‑151 (MQ=255)     tACTGGTTTCTTCATCAGCAGTCGCCAGCATGTAGAGCGGTTTTTTCAGTTTCGTGACGTCGCCATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTc                                                                                                                                         >  2:454314/1‑149 (MQ=255)     tACTGGTTTCTTCATCAGCAGTCGCCAGAATGTAGAGCGGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTc                                                                                                                                         >  2:530662/1‑149 (MQ=255)            ttcttcATCAGCAGTCGCCAGAATGTAGAGCGGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTaaa                                                                                                                                 <  2:309120/150‑1 (MQ=255)            ttcttcATCAGCAGTCGCCAGAATGTAGAGCGGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAAc                                                                                                                                >  2:300315/1‑151 (MQ=255)                tcatcaGCAGTCGCCAGAATGTAGAGCGGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACg                                                                                                                               <  2:200945/148‑1 (MQ=255)                tcatcaGCAGTCGCCAGAATGTAGAGCGGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGAt                                                                                                                            >  2:3836/1‑151 (MQ=255)                 catcaGCAGTCGCCAGAATGTAGAGCGGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGAt                                                                                                                            <  1:524028/150‑1 (MQ=255)                               aGAATGTAGAGCGGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGAc                                                                                                             <  1:630889/151‑1 (MQ=255)                                  aTGTAGAGCGGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCg                                                                                                           >  2:169574/1‑150 (MQ=255)                                    gTAGAGCGGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGtc                                                                                                         <  2:75868/150‑1 (MQ=255)                                         gCGGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCaaa                                                                                                   <  2:145463/151‑1 (MQ=255)                                                tttCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAGCGGATCGCGCGTCCAGCGACCGTCATCAAATGGCAc                                                                                             <  2:150729/150‑1 (MQ=255)                                                            gACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGc                                                                                >  2:6392/1‑151 (MQ=255)                                                              cGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGccc                                                                              >  1:567057/1‑151 (MQ=255)                                                                      tCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGc                                                                       >  2:301549/1‑150 (MQ=255)                                                                       cgcgTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGca                                                                      <  1:41073/150‑1 (MQ=255)                                                                          gTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGcaacaa                                                                  <  2:325784/151‑1 (MQ=255)                                                                          gTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGcaacaa                                                                  <  1:596006/151‑1 (MQ=255)                                                                             gcgcATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAg                                                                 <  2:552384/149‑1 (MQ=255)                                                                                   cAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGc                                                         >  1:201357/1‑151 (MQ=255)                                                                                     aGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGccc                                                       >  1:259016/1‑151 (MQ=255)                                                                                           aaCGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTc                                                 >  1:190578/1‑151 (MQ=255)                                                                                             cGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCg                                               >  2:283349/1‑151 (MQ=255)                                                                                                 aaaaaGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATAc                                           >  2:29887/1‑151 (MQ=255)                                                                                                  aaaaGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACt                                          >  1:478345/1‑151 (MQ=255)                                                                                                    aaGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCGTGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTg                                        <  2:181283/151‑1 (MQ=255)                                                                                                              tGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATAtt                              >  2:13112/1‑151 (MQ=255)                                                                                                               gTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATAtt                              >  2:277793/1‑150 (MQ=255)                                                                                                                         ccTAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGc                   >  2:35147/1‑151 (MQ=255)                                                                                                                           tAAGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCg                  >  1:570056/1‑150 (MQ=255)                                                                                                                            aaGCCGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAg                <  1:530659/151‑1 (MQ=255)                                                                                                                               ccGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTc             <  2:38944/151‑1 (MQ=255)                                                                                                                                cGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTg                                                   >  2:301355/1‑112 (MQ=255)                                                                                                                                cGTAAAGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTcc            >  1:27136/1‑151 (MQ=255)                                                                                                                                   aaaGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGCCCAGCGACCGTCATCAAAGGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTg                                                   <  1:301353/109‑1 (MQ=255)                                                                                                                                    aaGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATTCTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGa        <  2:393115/151‑1 (MQ=255)                                                                                                                                    aaGCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGa        <  2:77509/151‑1 (MQ=255)                                                                                                                                      gCTTGCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAAc      >  1:71048/1‑151 (MQ=255)                                                                                                                                          gCCGTCATGCCCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGc  >  2:594310/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                             GCCATACTGGTTTCTTCATCAGCAGTCGCCAGAATGTAGAGCGGTTTTTTCAGTTTCGTGACGTCGACATCGCGTAGCGCATCAAGGATAAACGCAAAAAAGCCTTTCATGTCGGCGGTGCCTAAGCCGTAAAGCTTGCCGTCATGCTCCGTCAGTGTAAACGGATCGCGCGTCCAGCGACCGTCATCAAATGGCACCGTATCGGTATGCCCCGCCAGCAACAAGCCGCCAGCCCCCTGTCCGATACTTGCCAGCATATTGAATTTGTTGCGAGTTCCTGGAACAGGC  >  NC_000913/4154402‑4154689

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