Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,334,049	T→C	Q514Q (CAA→CAG)	eptA ←	lipid A phosphoethanolamine transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,334,049	0	T	C	100.0%	186.9 / NA	50	Q514Q (CAA→CAG)	eptA	lipid A phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/22);  total (28/22)

CAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCTTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACAGCAGCATCGGCACCTGTTTTTGGCTATCCGGGGCGATGGCATAAGGCAGACCGTGC  >  NC_000913/4333903‑4334190                                                                                                                                                   |                                                                                                                                              cAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCtgt                                                                                                                                            <  2:40992/150‑1 (MQ=255)      ccTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAg                                                                                                                                      <  1:676277/151‑1 (MQ=255)        tGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAgtg                                                                                                                                    >  1:571639/1‑151 (MQ=255)         gCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAgt                                                                                                                                     <  1:462677/149‑1 (MQ=255)           aaTAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGtt                                                                                                                                  >  2:268216/1‑150 (MQ=255)           aaTAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTg                                                                                                                                 <  2:390495/151‑1 (MQ=255)              aaCAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGt                                                                                                                              >  2:637865/1‑151 (MQ=255)                 aGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGttt                                                                                                                            >  2:87533/1‑150 (MQ=255)                 aGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTg                                                                                                                           >  1:555075/1‑151 (MQ=255)                   cgaaTCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTgcg                                                                                                                         <  2:554105/148‑1 (MQ=255)                        cGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGt                                                                                                                    <  2:71722/151‑1 (MQ=255)                          tCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGCATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGttt                                                                                                                  >  2:595732/1‑151 (MQ=255)                              caacaaTCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTgc                                                                                                              <  2:415500/151‑1 (MQ=255)                               aacaaTCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTgca                                                                                                             <  1:121539/151‑1 (MQ=255)                                     cAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGt                                                                                                       >  1:619254/1‑151 (MQ=255)                                      aGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGtt                                                                                                      >  2:125631/1‑151 (MQ=255)                                           tttCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGt                                                                                                 >  1:155008/1‑151 (MQ=255)                                           tttCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGt                                                                                                 >  2:173292/1‑151 (MQ=255)                                              cattcaCTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCa                                                                                               >  2:9328/1‑150 (MQ=255)                                              cattcaCTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCa                                                                                               <  1:7885/150‑1 (MQ=255)                                                 tcactcacTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAAcc                                                                                            <  1:593228/150‑1 (MQ=255)                                                    ctcactcTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGa                                                                                         <  1:134703/150‑1 (MQ=255)                                                     tcactcTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGAt                                                                                        >  2:28649/1‑150 (MQ=255)                                                      cactcTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTTCAGGCAGTTCTGGTCAACCTGAt                                                                                        <  1:272557/149‑1 (MQ=255)                                                        ctctCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCg                                                                                    >  1:448127/1‑151 (MQ=255)                                                           tcCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGttt                                                                                 >  2:615792/1‑151 (MQ=255)                                                            ccTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGtttt                                                                                <  2:182470/151‑1 (MQ=255)                                                               gCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGa                                                                             <  2:12069/151‑1 (MQ=255)                                                                 aaGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGAt                                                                            >  1:637963/1‑150 (MQ=255)                                                                   gTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAAt                                                                         >  1:558640/1‑151 (MQ=255)                                                                       gCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGCGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAAtcctc                                                                     <  1:197975/151‑1 (MQ=255)                                                                              atcatcCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATGCCGTTTTTGATAATCCTCCGACAg                                                               <  1:515703/150‑1 (MQ=255)                                                                              atcatcCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGc                                                              >  2:132129/1‑151 (MQ=255)                                                                                     gCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCAcagca                                                       <  1:114035/151‑1 (MQ=255)                                                                                       aGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACagcagc                                                     >  2:69832/1‑151 (MQ=255)                                                                                       aGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACagcagc                                                     <  2:38844/151‑1 (MQ=255)                                                                                       aGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACagcagc                                                     >  1:451349/1‑151 (MQ=255)                                                                                          cTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTCCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCAGACAGCCACagcagc                                                     >  2:405907/1‑148 (MQ=255)                                                                                           tGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACAGCAGCATCg                                                 >  1:393885/1‑151 (MQ=255)                                                                                                 tACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACAGCAGCATCGGCACCt                                           >  2:118919/1‑151 (MQ=255)                                                                                                      cGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACAGCAGCATCGGCACCTGtttt                                      >  2:596665/1‑151 (MQ=255)                                                                                                           cAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACAGCAGCATCGGCACCTGTTTTTgg                                   <  2:122054/149‑1 (MQ=255)                                                                                                           cAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACAGCAGCATCGGCACCTGTTTTTGGCt                                 <  1:475481/151‑1 (MQ=255)                                                                                                                       aTCCCAATAGCGTGGACAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACAGCAGCATCGGCACCTGTTTTTGGCTATCCGGGGCGAt                     <  2:417561/151‑1 (MQ=255)                                                                                                                         cccAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACAGCAGCATCGGCACCTGTTTTTGGCTATCCGGTGCGATg                    >  2:344750/1‑150 (MQ=255)                                                                                                                          ccAATAGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACAGCAGCATCGGCACCTGTTTTTGGCTATCCGGGGCGATGGc                  <  2:592075/151‑1 (MQ=255)                                                                                                                               aGCGTGGAGAATAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACAGCAGCATCGGCACCTGTTTTTGGCTATCCGGGGCGATGGCATa               >  1:243039/1‑149 (MQ=255)                                                                                                                                        aaTAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACAGCAGCATCGGCACCTGTTTTTGGCTATCCGGGGCGATGGCATAAGGCAGACCGt    >  1:263769/1‑151 (MQ=255)                                                                                                                                          tAAATTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACAGCAGCATCGGCACCTGTTTTTGGCTATCCGGGGCGATGGCATAAGGCAGACCGTGc  >  2:175834/1‑151 (MQ=255)                                                                                                                                           aaaTTGTCCTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACAGCAGCATCGGCACCTGTTTTTGGCTATCCGGGGCGATGGCATAAGGCAGACCGTGc  >  2:302659/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                              CAGTCCCTGCAATAACAGCGTATCGTCTTCAACAATCAGAATTTTCATTCACTCACTCTCCTGCAAGTTTGCAGAATATCATCCGCAGCCTGGTAATACTTCGTCTCAACGCCAGTTAATCCCAATAGCGTGGAGAATAAATTGTCTTGTGAATAGTGTTGCGTTTGCGCCTGTTTTTGCAGGCAGTTCTGGTCAACCTGATACCGTTTTTGATAATCCTCCGACAGCCACAGCAGCATCGGCACCTGTTTTTGGCTATCCGGGGCGATGGCATAAGGCAGACCGTGC  >  NC_000913/4333903‑4334190

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