Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,517,815	T→C	E224E (GAA→GAG)	aroB ←	3‑dehydroquinate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,517,815	0	T	C	100.0%	213.4 / NA	58	E224E (GAA→GAG)	aroB	3‑dehydroquinate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/27);  total (31/27)

CAGCGACCGCTTCACCATGTAACCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGTTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTCAAGAATAATGCCGTATTTGATGACTTCTGCCAGCCC  >  NC_000913/3517679‑3517946                                                                                                                                         |                                                                                                                                    cAGCGACCGCTTCACCATGTAACCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTca                                                                                                                                  >  2:938474/1‑140 (MQ=255) cAGCGACCGCTTCACCATGTAACCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTca                                                                                                                                  >  1:1410949/1‑140 (MQ=255) cAGCGACCGCTTCACCATGTAACCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTca                                                                                                                                  >  1:342165/1‑140 (MQ=255)  aGCGACCGCTTCACCATGTAACCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTcac                                                                                                                                 >  1:816784/1‑140 (MQ=255)   gCGACCGCTTCACCATGTAACCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTcaca                                                                                                                                >  1:1541606/1‑140 (MQ=255)      aCCGCTTCACCATGTAACCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCacaaca                                                                                                                             >  1:935861/1‑140 (MQ=255)      aCCGCTTCACCATGTAACCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCacaaca                                                                                                                             <  2:1476566/140‑1 (MQ=255)      aCCGCTTCACCATGTAACCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCacaaca                                                                                                                             >  1:167159/1‑140 (MQ=255)         gCTTCACCATGTAACCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACg                                                                                                                          <  2:1387671/140‑1 (MQ=255)              aCCATGTAACCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGaa                                                                                                                     >  2:127593/1‑140 (MQ=255)              aCCATGTAACCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGaa                                                                                                                     <  2:1203549/140‑1 (MQ=255)                cATGTAACCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATa                                                                                                                   <  2:701583/140‑1 (MQ=255)                cATGTAACCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATa                                                                                                                   <  2:154642/140‑1 (MQ=255)                     aaCCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTa                                                                                                              >  1:335123/1‑140 (MQ=255)                     aaCCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTa                                                                                                              >  1:412226/1‑140 (MQ=255)                       ccAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACg                                                                                                            >  1:171213/1‑140 (MQ=255)                       ccAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACg                                                                                                            >  1:1381086/1‑140 (MQ=255)                        cAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCTCGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGc                                                                                                           <  2:1059497/140‑1 (MQ=255)                        cAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGc                                                                                                           <  2:1571401/140‑1 (MQ=255)                         aaTTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGcc                                                                                                          >  2:1024918/1‑140 (MQ=255)                                  aCCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGa                                                                                                 >  1:1547565/1‑140 (MQ=255)                                  aCCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGa                                                                                                 >  1:122715/1‑140 (MQ=255)                                  aCCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGa                                                                                                 <  2:432339/140‑1 (MQ=255)                                  aCCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGa                                                                                                 <  2:1407778/140‑1 (MQ=255)                                      cATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGt                                                                                             >  1:350401/1‑140 (MQ=255)                                         ttCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCa                                                                                          <  2:365187/140‑1 (MQ=255)                                                  aaTGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGcaaca                                                                                 >  1:1443054/1‑140 (MQ=255)                                                  aaTGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGcaaca                                                                                 >  1:1538997/1‑140 (MQ=255)                                                     ggCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACg                                                                              <  2:19485/140‑1 (MQ=255)                                                               aaGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGAtt                                                                    <  2:171213/140‑1 (MQ=255)                                                                     tgtCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTc                                                              >  1:1202332/1‑140 (MQ=255)                                                                     tgtCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTc                                                              >  1:167287/1‑140 (MQ=255)                                                                     tgtCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAACTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTc                                                              >  1:439396/1‑140 (MQ=255)                                                                         ccAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTccagc                                                          <  2:935861/140‑1 (MQ=255)                                                                         ccAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTccagc                                                          <  2:342165/140‑1 (MQ=255)                                                                          cAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTccagcc                                                         >  2:780213/1‑140 (MQ=255)                                                                          cAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTccagcc                                                         >  1:1260569/1‑140 (MQ=255)                                                                                 aGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGTGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTaaa                                                  >  1:249513/1‑140 (MQ=255)                                                                                        cACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCa                                           <  2:1381086/140‑1 (MQ=255)                                                                                        cACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCa                                           >  1:1291769/1‑140 (MQ=255)                                                                                        cACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCa                                           <  2:167159/140‑1 (MQ=255)                                                                                        cACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGGCCAGACGCAACAACGCATCCAGATTCTCTTCCCGCCAGTTAAAAAACGCa                                           >  2:1176879/1‑140 (MQ=255)                                                                                          cGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCAcc                                         <  2:335123/140‑1 (MQ=255)                                                                                             aaCCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTc                                      >  1:1426393/1‑140 (MQ=255)                                                                                             aaCCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTc                                      <  2:367559/140‑1 (MQ=255)                                                                                             aaCCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTc                                      <  2:1431391/140‑1 (MQ=255)                                                                                             aaCCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTc                                      >  1:720128/1‑140 (MQ=255)                                                                                                  ggTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTCAAGaa                                 >  1:320212/1‑140 (MQ=255)                                                                                                  ggTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTCAAGaa                                 >  1:1408657/1‑140 (MQ=255)                                                                                                    tttCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTCAAGaata                               <  2:96705/140‑1 (MQ=255)                                                                                                      tCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTCAAGaataat                             <  1:1024918/140‑1 (MQ=255)                                                                                                        gcgcTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTCAAGAATAATGc                           <  2:618666/140‑1 (MQ=255)                                                                                                         cgcTCGTCGGCGGCGACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTCAAGAATAATGcc                          <  2:1410949/140‑1 (MQ=255)                                                                                                                     gcgACAACTTCTGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTCAAGAATAATGCCGTATTTGATGAc              <  2:1443054/140‑1 (MQ=255)                                                                                                                                tGCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTCAAGAATAATGCCGTATTTGATGACTTCTGCCAGcc   >  2:566770/1‑140 (MQ=255)                                                                                                                                 gCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTCAAGAATAATGCCGTATTTGATGACTTCTGCCAGccc  <  2:1291769/140‑1 (MQ=255)                                                                                                                                 gCCTTCAGCTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTCAAGAATAATGCCGTATTTGATGACTTCTGCCAGccc  <  2:350401/140‑1 (MQ=255)                                                                                                                                 gCCTTCAACTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTCAAGAATAATGCCGTATTTGATGACTTCTGCCAGccc  <  2:439396/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                    CAGCGACCGCTTCACCATGTAACCAATTGCCATACCCCATTTCAGCTTCAATGGCATGACCAAAGGTGTGTCCCAGATTCAGTAAAGCACGTAACCCGGTTTCGCGCTCGTCGGCGGCGACAACTTCTGCCTTCAGTTCACAACAACGGCGAATACAGTACGCCATTGCCGGACCGTCCAGACGCAACAACGCATCCAGATTCTCTTCCAGCCAGTTAAAAAACGCACCGTCAAGAATAATGCCGTATTTGATGACTTCTGCCAGCCC  >  NC_000913/3517679‑3517946

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