Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,153,753	A→G	F365F (TTT→TTC)	argE ←	acetylornithine deacetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,153,753	0	A	G	100.0%	86.1 / NA	25	F365F (TTT→TTC)	argE	acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/13);  total (12/13)

ACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCT  >  NC_000913/4153627‑4153882                                                                                                                               |                                                                                                                                  aCATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATGAACCgt                                                                                                                             <  1:644985/133‑1 (MQ=255)          cGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGa                                                                                                                    <  1:408551/133‑1 (MQ=255)          cGGCACTGCAGTTTGCCACCAGAGCCGGATCCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGa                                                                                                                    <  1:169779/133‑1 (MQ=38)                 ggAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATtcat                                                                                                             >  2:74285/1‑133 (MQ=255)                  gAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATtcatc                                                                                                            <  1:74285/133‑1 (MQ=255)                              aGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGc                                                                                                >  2:569634/1‑133 (MQ=255)                                cgcCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCt                                                                                              <  1:576739/133‑1 (MQ=255)                                  ccGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGa                                                                                            >  2:632964/1‑133 (MQ=255)                                            cgcgAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGc                                                                                  <  1:632964/133‑1 (MQ=255)                                                 gcgcCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGc                                                                             <  1:569634/133‑1 (MQ=255)                                                             gCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGc                                                                 >  2:151910/1‑133 (MQ=255)                                                               cTACGTTTTATTGCCCGCAAAAATGGTGAATTCCCTGGGTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGt                                                               <  2:200666/133‑1 (MQ=255)                                                                            ccAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCa                                                         >  1:563340/1‑126 (MQ=255)                                                                            ccAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCa                                                         <  2:563340/126‑1 (MQ=255)                                                                                         tGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGgcg                                     <  1:151910/133‑1 (MQ=255)                                                                                                    ggTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACag                          <  2:444083/133‑1 (MQ=255)                                                                                                     gTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATGACGTTTGAATAAACGGCGCTTGGGTACagt                         >  1:172278/1‑133 (MQ=255)                                                                                                     gTTATCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAACGTCGGGCTTAACGTTCGAATAAACCGCGCTTCGGTACagt                         >  2:263291/1‑133 (MQ=255)                                                                                                        aTCAGTTCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATGAACAGCGCTTCGGTACagtagt                      >  1:207489/1‑133 (MQ=255)                                                                                                           aGTTCGCGAGTGGGCTTGATGATCCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTca                   <  1:263291/133‑1 (MQ=255)                                                                                                              tCGCGGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTcacca                >  2:352478/1‑133 (MQ=255)                                                                                                               cgcgGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTcaccac               >  2:100272/1‑133 (MQ=255)                                                                                                               cgcgGGTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCAAAACGTTTGAATAAACGACGCTTCGGAACAGCAGTTcaccac               >  2:225774/1‑133 (MQ=255)                                                                                                                  gggTGGGCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTc            <  1:106683/133‑1 (MQ=255)                                                                                                                       ggCTTGATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGTCTGATTAATTGAGACAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGttt       >  2:563610/1‑133 (MQ=255)                                                                                                                            gATGAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCt  <  1:547978/133‑1 (MQ=255)                                                                                                                               |                                                                                                                                  ACATTATTGCGGCACTGGAGTTTGGCAACAGCGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGTTTTAATGCCAGCAAAAATGGTGAATTACCTGGGTTATCAGTTCGCGGGTGGGCTTGATAAACCGTGTTTCCAGATATTCATCAGGTTGATGAGCCTGATTAATTGAGCCAGGCCCCAACACCAGCGTCGGGCATAACGTTTGAATAAACGGCGCTTCGGTACAGTAGTTCACCACTTCGGTTTTTGCT  >  NC_000913/4153627‑4153882

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