Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	30,616	T→C	P322P (CCT→CCC)	carA →	carbamoyl phosphate synthetase small subunit, glutamine amidotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	30,616	0	T	C	100.0%	68.5 / NA	20	P322P (CCT→CCC)	carA	carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/10);  total (10/10)

GGCGAGCGGTGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACACGACGCCGCGCCGT  >  NC_000913/30481‑30746                                                                                                                                        |                                                                                                                                   ggcgAGCGGTGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCa                                                                                                                                 >  2:652203/1‑139 (MQ=255)          tGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCgt                                                                                                                        >  1:715658/1‑139 (MQ=255)                    gTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATa                                                                                                              >  1:230592/1‑139 (MQ=255)                    gTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATa                                                                                                              >  1:386472/1‑139 (MQ=255)                                 ttGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGa                                                                                                 <  1:652203/139‑1 (MQ=255)                                 ttGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGa                                                                                                 <  1:74238/139‑1 (MQ=255)                                 ttGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGa                                                                                                 <  1:778115/139‑1 (MQ=255)                                     tCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGt                                                                                             <  2:715658/139‑1 (MQ=255)                                                      aTCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCa                                                                            <  2:199182/139‑1 (MQ=255)                                                                 gATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATa                                                                 <  1:486660/139‑1 (MQ=255)                                                                     tgGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATccc                                                                                                        <  1:398277/96‑1 (MQ=255)                                                                     tgGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATccc                                                                                                        >  2:398277/1‑96 (MQ=255)                                                                     tgGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAAcc                                                             >  1:243670/1‑139 (MQ=255)                                                                      ggAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCg                                                            >  2:658891/1‑139 (MQ=255)                                                                                                cccAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGa                                  >  2:454734/1‑139 (MQ=255)                                                                                                          cGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCt                                    <  2:167603/127‑1 (MQ=255)                                                                                                          cGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCt                                    >  1:167603/1‑127 (MQ=255)                                                                                                             ttttGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGt                     <  2:230592/139‑1 (MQ=255)                                                                                                                            aGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACACGAcgccgcg      >  2:788476/1‑139 (MQ=255)                                                                                                                                aCATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACACGACGCCGCGCCgt  <  2:386472/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   GGCGAGCGGTGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACACGACGCCGCGCCGT  >  NC_000913/30481‑30746

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