Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,154,941	C→A	intergenic (‑94/‑60)	argE ← / → argC	acetylornithine deacetylase/N‑acetylglutamylphosphate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,154,941	0	C	A	100.0%	110.6 / NA	32	intergenic (‑94/‑60)	argE/argC	acetylornithine deacetylase/N‑acetylglutamylphosphate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (21/11);  total (21/11)

TCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTAGTGACCTATGTAAATCGCCATCCGC  >  NC_000913/4154817‑4155076                                                                                                                             |                                                                                                                                         tctcGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGttt                                                                                                                           <  2:463563/140‑1 (MQ=255)                   aTTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAAc                                                                                                         <  2:71529/140‑1 (MQ=255)                                          ccTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATaagaa                                                                                             >  1:232861/1‑129 (MQ=255)                                          ccTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATaagaa                                                                                             <  2:232861/129‑1 (MQ=255)                                               ggTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGAt                                                                            <  2:523929/140‑1 (MQ=255)                                                   aTGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTg                                                                        >  2:473627/1‑140 (MQ=255)                                                   aTGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTg                                                                        >  1:47184/1‑140 (MQ=255)                                                     gATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGaa                                                                      >  2:291401/1‑140 (MQ=255)                                                          tATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGc                                                                 >  1:2109/1‑140 (MQ=255)                                                                  ttCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATtgtg                                                         <  1:291401/140‑1 (MQ=255)                                                                   tCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTgg                                                        <  2:402105/140‑1 (MQ=255)                                                                    cATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTggg                                                       >  1:515813/1‑140 (MQ=255)                                                                    cATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTggg                                                       >  1:453287/1‑140 (MQ=255)                                                                      tGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTg                                                     >  2:350239/1‑140 (MQ=255)                                                                      tGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTg                                                     >  1:272026/1‑140 (MQ=255)                                                                      tGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACC‑GGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGc                                                    >  1:364402/1‑140 (MQ=255)                                                                        cAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTATCCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGcc                                                   <  1:270398/140‑1 (MQ=255)                                                                           tATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGc                                                >  2:6071/1‑140 (MQ=255)                                                                                   aaTAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGc                                        >  2:248703/1‑140 (MQ=255)                                                                                   aaTAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGc                                        >  1:512259/1‑140 (MQ=255)                                                                                     tAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTg                                      >  2:285285/1‑140 (MQ=255)                                                                                     tAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTg                                      >  1:145041/1‑140 (MQ=255)                                                                                       aaaaTACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGgc                                    >  1:21489/1‑140 (MQ=255)                                                                                           tACACTAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCag                                <  1:23278/140‑1 (MQ=255)                                                                                                tAACGTTGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTa                           <  2:512072/140‑1 (MQ=255)                                                                                                     ttGAGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTAGTGAc                      <  2:21489/140‑1 (MQ=255)                                                                                                        aGCGTAATAAAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTAGTGACCTa                   >  1:190230/1‑140 (MQ=255)                                                                                                        aGCGTAATAAAACCCA‑CAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCCGGGGCAAAGCTAGGGACCTaa                  >  1:50172/1‑139 (MQ=255)                                                                                                            taataaAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATAGGCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTAGTGACCTATGt                >  2:247109/1‑140 (MQ=255)                                                                                                               taaAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTa                           <  1:365629/125‑1 (MQ=255)                                                                                                               taaAACCCACCAGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTa                           >  2:365629/1‑125 (MQ=255)                                                                                                                       accaGCAGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTAGCGACCTATGTAAATCGCCATcc    >  1:182804/1‑140 (MQ=255)                                                                                                                         cagcagTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTAGTGACCTATGTAAATCGCCATCCGc  >  2:14982/1‑140 (MQ=255)                                                                                                                             |                                                                                                                                         TCTCGATAAATGGCGGTAATTTGTTTTTCATTGTTGACACACCTCTGGTCATGATAGTATCAATATTCATGCAGTATTTATGAATAAAAATACACTAACGTTGAGCGTAATAAAACCCACCAGCCGTAAGGTGAATGTTTTACGTTTAACCTGGCAACCAGACATAAGAAGGTGAATA‑GCCCCGATGTTGAATACGCTGATTGTGGGTGCCAGCGGCTACGCTGGCGCAGAGCTAGTGACCTATGTAAATCGCCATCCGC  >  NC_000913/4154817‑4155076

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑