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

GGCGAGCGGTGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACACGACGCCGCGC  >  NC_000913/30481‑30743                                                                                                                                        |                                                                                                                                ggcgAGCGGTGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCa                                                                                                                              <  1:661117/139‑1 (MQ=255)        ggTGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGc                                                                                                                       <  2:314665/139‑1 (MQ=255)           gCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCgtg                                                                                                                    <  2:1187658/139‑1 (MQ=255)              aaGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCa                                                                                                                 <  2:1604098/139‑1 (MQ=255)                  cTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCa                                                                                                             <  1:1347899/139‑1 (MQ=255)                    gTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATa                                                                                                           >  2:226890/1‑139 (MQ=255)                    gTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATa                                                                                                           >  2:502789/1‑139 (MQ=255)                    gTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGATAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCTCGGTTTTGCGGTGGACGAAGCAACATTACGCGCAAACCTGCGTGTCACGCATa                                                                                                           >  2:1444355/1‑139 (MQ=255)                           tGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCt                                                                                                    >  1:604475/1‑139 (MQ=255)                                 ttGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGa                                                                                              <  2:1367575/139‑1 (MQ=255)                                              cggcAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAg                                                                                 >  1:749955/1‑139 (MQ=255)                                                   aCCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCAt                                                                            >  2:779883/1‑139 (MQ=255)                                                       tCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCAt                                                                        >  1:1406982/1‑139 (MQ=255)                                                       tCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCAt                                                                        >  2:441869/1‑139 (MQ=255)                                                         cGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCg                                                                      <  1:1577457/139‑1 (MQ=255)                                                         cGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCg                                                                      <  2:1315011/139‑1 (MQ=255)                                                          ggTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTATGTTACAGGGCATTCATCGc                                                                     <  2:1172164/139‑1 (MQ=255)                                                          ggTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGc                                                                     <  1:1115436/139‑1 (MQ=255)                                                              aaaGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACg                                                                                       <  1:1485762/117‑1 (MQ=255)                                                              aaaGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACg                                                                                       >  2:1485762/1‑117 (MQ=255)                                                                     tgGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAGCATTACCCGCAAACCTGCGTGTCAAGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAAcc                                                          >  2:1142359/1‑139 (MQ=255)                                                                     tgGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAAcc                                                          >  2:1705843/1‑139 (MQ=255)                                                                         gaAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCa                                                      <  2:338745/139‑1 (MQ=255)                                                                          aaaaaCGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTc                                                                                                                  >  2:123876/1‑78 (MQ=255)                                                                          aaaaaCGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTc                                                                                                                  <  1:123876/78‑1 (MQ=255)                                                                            aaaCGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCa                                                                                                                 >  2:512574/1‑77 (MQ=255)                                                                            aaaCGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCa                                                                                                                 <  1:512574/77‑1 (MQ=255)                                                                            aaaCGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTa                                                                                         >  1:90712/1‑101 (MQ=255)                                                                            aaaCGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTa                                                                                         <  2:90712/101‑1 (MQ=255)                                                                             aaCGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCa                                                  <  2:1406982/139‑1 (MQ=255)                                                                                    tAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTa                                                                                         <  2:1514002/93‑1 (MQ=255)                                                                                    tAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTa                                                                                         >  1:1514002/1‑93 (MQ=255)                                                                                        gATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAgggg                                       <  1:226890/139‑1 (MQ=255)                                                                                           cACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCAc                                    >  2:843338/1‑139 (MQ=255)                                                                                            aCCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATa                                                              >  1:68614/1‑112 (MQ=255)                                                                                            aCCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATa                                                              <  2:68614/112‑1 (MQ=255)                                                                                              cGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGc                                      <  1:1701286/134‑1 (MQ=255)                                                                                              cGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGc                                      >  2:1701286/1‑134 (MQ=255)                                                                                              cGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCt                                 >  2:510651/1‑139 (MQ=255)                                                                                              cGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCt                                 >  1:236904/1‑139 (MQ=255)                                                                                                cccAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGt                                                                                          >  2:1112661/1‑80 (MQ=255)                                                                                                cccAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGt                                                                                          <  1:1112661/80‑1 (MQ=255)                                                                                                   aGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGa                                                                                              >  2:364010/1‑73 (MQ=255)                                                                                                   aGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGa                                                                                              <  1:364010/73‑1 (MQ=255)                                                                                                     aaCCACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAagcca                          >  1:1071087/1‑139 (MQ=255)                                                                                                        cACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAagccagcc                       <  2:749955/139‑1 (MQ=255)                                                                                                        cACGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAagccagcc                       >  1:1365814/1‑139 (MQ=255)                                                                                                          cGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTa                                                                                    >  2:108636/1‑76 (MQ=255)                                                                                                          cGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTa                                                                                    <  1:108636/76‑1 (MQ=255)                                                                                                          cGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCt                     >  2:965588/1‑139 (MQ=255)                                                                                                          cGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCt                     >  2:1320013/1‑139 (MQ=255)                                                                                                          cGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCt                     >  1:360224/1‑139 (MQ=255)                                                                                                          cGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACACATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCg                                                                      >  2:198022/1‑90 (MQ=255)                                                                                                          cGGTTTTGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACACATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCg                                                                      <  1:198022/90‑1 (MQ=255)                                                                                                               ttGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGt                                                                                          <  2:39204/65‑1 (MQ=255)                                                                                                               ttGCGGTGGACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGt                                                                                          >  1:39204/1‑65 (MQ=255)                                                                                                                       gACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTg                    >  1:870053/1‑127 (MQ=255)                                                                                                                       gACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTg                    <  2:870053/127‑1 (MQ=255)                                                                                                                       gACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACacgacg        >  1:525437/1‑139 (MQ=255)                                                                                                                       gACGAAGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACacgacg        <  2:525437/139‑1 (MQ=255)                                                                                                                            aGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGc                                      >  2:116840/1‑104 (MQ=255)                                                                                                                            aGCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGc                                      <  1:116840/104‑1 (MQ=255)                                                                                                                             gCAACATTACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTTGTCCACACGAcgccgcgc  <  1:1142359/139‑1 (MQ=255)                                                                                                                                    tACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAagccagcc                       <  2:1249914/111‑1 (MQ=255)                                                                                                                                    tACCCGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAagccagcc                       >  1:1249914/1‑111 (MQ=255)                                                                                                                                        |                                                                                                                                GGCGAGCGGTGCGAAGACTGTCAAAATGAAATTTGGTCACCACGGCGGCAACCATCCGGTTAAAGATGTGGAGAAAAACGTGGTAATGATCACCGCCCAGAACCACGGTTTTGCGGTGGACGAAGCAACATTACCTGCAAACCTGCGTGTCACGCATAAATCCCTGTTCGACGGTACGTTACAGGGCATTCATCGCACCGATAAACCGGCATTCAGCTTCCAGGGGCACCCTGAAGCCAGCCCTGGTCCACACGACGCCGCGC  >  NC_000913/30481‑30743

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