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%	178.3 / NA	51	P322P (CCT→CCC)	carA	carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/23);  total (28/23)

TGGCGAGCGGTGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACACGACGCCGCGCCGTT  >  NC_000913/30480‑30747                                                                                                                                         |                                                                                                                                    tggcgAGCGGTGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGc                                                                                                                                   >  1:333192/1‑139 (MQ=255)  ggcgAGCGGTGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCa                                                                                                                                  <  2:481254/139‑1 (MQ=255)       gCGGTGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCt                                                                                                                             <  1:505403/139‑1 (MQ=255)             cGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCgtgt                                                                                                                       <  1:1231032/139‑1 (MQ=255)                aGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCAc                                                                                                                    <  1:1122172/139‑1 (MQ=255)                          aaTGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATcc                                                                                                          <  2:626446/139‑1 (MQ=255)                            tGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGc                                                                                                                           >  1:1204483/1‑120 (MQ=255)                            tGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGc                                                                                                                           <  2:1204483/120‑1 (MQ=255)                              aaaTTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGt                                                                                                      >  2:1215891/1‑139 (MQ=255)                                  ttGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGa                                                                                                  <  2:307362/139‑1 (MQ=255)                                        accacGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTAc                                                                                            >  2:250130/1‑139 (MQ=255)                                                    aCCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGGAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATGCCTGTGGGACGGTACGTTACAGGGCAt                                                                                >  1:349056/1‑139 (MQ=255)                                                           ggTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGc                                                                         <  1:1654550/139‑1 (MQ=255)                                                                      tgGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAAcc                                                              >  2:471773/1‑139 (MQ=255)                                                                            aaaaCGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCAtt                                                        <  1:296326/139‑1 (MQ=255)                                                                             aaaCGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCa                                                                                                                                  >  2:1570335/1‑64 (MQ=255)                                                                             aaaCGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCa                                                                                                                                  <  1:1570335/64‑1 (MQ=255)                                                                                   ggTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATaaa                                                                                                             >  2:816035/1‑79 (MQ=255)                                                                                   ggTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATaaa                                                                                                             <  1:816035/79‑1 (MQ=255)                                                                                   ggTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTc                                                 >  1:816261/1‑139 (MQ=255)                                                                                          aTCACCGCCCAGAACCACGGTTTTGTGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGc                                          <  2:594515/139‑1 (MQ=255)                                                                                           tCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCa                                         >  2:3805/1‑139 (MQ=255)                                                                                              ccGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCAccc                                      >  1:42763/1‑139 (MQ=255)                                                                                              ccGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCAccc                                      >  1:1578398/1‑139 (MQ=255)                                                                                               cGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGa                                                                                                  <  2:738899/78‑1 (MQ=255)                                                                                               cGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCt                                     <  1:3805/139‑1 (MQ=255)                                                                                               cGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCt                                     >  2:1525554/1‑139 (MQ=255)                                                                                               cGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCt                                     >  1:851954/1‑139 (MQ=255)                                                                                               cGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGATCGa                                                                                                  >  1:738899/1‑78 (MQ=255)                                                                                                gCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGt                                                                                              >  2:712619/1‑81 (MQ=255)                                                                                                gCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGt                                                                                              <  1:712619/81‑1 (MQ=255)                                                                                                 cccAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGa                                   >  1:505523/1‑139 (MQ=255)                                                                                                   cAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAag                                 >  2:745952/1‑139 (MQ=255)                                                                                                    aGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAagc                                <  2:1369360/139‑1 (MQ=255)                                                                                                      aaCCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTc                                                                                                                      >  2:367261/1‑51 (MQ=255)                                                                                                      aaCCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTc                                                                                                                      <  1:367261/51‑1 (MQ=255)                                                                                                      aaCCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCg                                                                          <  1:402671/95‑1 (MQ=255)                                                                                                      aaCCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCg                                                                          >  2:402671/1‑95 (MQ=255)                                                                                                         cACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAagccagcc                           >  1:206383/1‑139 (MQ=255)                                                                                                           cGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCt                         >  1:638940/1‑139 (MQ=255)                                                                                                                ttGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTcc                    >  1:703169/1‑139 (MQ=255)                                                                                                                 tGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCAATCATCg                                                                          >  1:1022167/1‑84 (MQ=255)                                                                                                                 tGCGGTGGACGAAGCAACATTACCCGCAAACCAGCGTGTCACTCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCg                                                                          <  2:1022167/84‑1 (MQ=255)                                                                                                                   cGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAgggg                                           >  2:1382207/1‑113 (MQ=255)                                                                                                                   cGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAgggg                                           <  1:1382207/113‑1 (MQ=255)                                                                                                                        gACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACa                                                                                      <  1:1111656/65‑1 (MQ=255)                                                                                                                        gACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACa                                                                                      >  2:1111656/1‑65 (MQ=255)                                                                                                                              gCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACACGAcgccgcgc      <  2:333192/139‑1 (MQ=255)                                                                                                                                  cATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCacac                <  1:981348/125‑1 (MQ=255)                                                                                                                                  cATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCacac                >  2:981348/1‑125 (MQ=255)                                                                                                                                  cATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACACGACGCCGCGCCgtt  >  1:1604430/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                    TGGCGAGCGGTGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACACGACGCCGCGCCGTT  >  NC_000913/30480‑30747

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