Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,122,826	G→A	Q150* (CAA→TAA)	yceB ←	lipoprotein, DUF1439 family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,122,826	0	G	A	100.0%	212.2 / NA	53	Q150* (CAA→TAA)	yceB	lipoprotein, DUF1439 family
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/29);  total (24/29)

ACGAAAAATAAAGGGCGCGAATGCGCCCTCGTGATTAATCAGTAAATGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTGCTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAGATCGCACCTTTTTCTTTATCAAACACGGGCA  >  NC_000913/1122679‑1122972                                                                                                                                                    |                                                                                                                                                   aCGAAAAATAAAGGGCGCGAATGCGCCCTCGTGATTAATCAGTAAATGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTg                                                                                                                                                 >  1:465351/1‑151 (MQ=255)  cGAAAAATAAAGGGCGCGAATGCGCCCTCGTGATTAATCAGTAAATGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTgg                                                                                                                                                >  1:68508/1‑151 (MQ=255)     aaaaTAAAGGGCGCGAATGCGCCCTCGTGATTAATCAGTAAATGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTa                                                                                                                                             <  2:256499/151‑1 (MQ=255)        aTAAAGGGCGCGAATGCGCCCTCGTGATTAATCAGTAAATGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTaa                                                                                                                                            >  1:509543/1‑149 (MQ=255)               gcgcGAATGCGCCCTCGTGATTAATCAGTAAATGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGc                                                                                                                                   <  1:617884/151‑1 (MQ=255)                cgcgAATGCGCCCTCGTGATTAATCAGTAAATGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCg                                                                                                                                  >  2:674284/1‑151 (MQ=255)                         gcCCTCGTGATTAATCAGTAAATGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTg                                                                                                                         >  2:26078/1‑151 (MQ=255)                               gTGATTAATCAGTAAATGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCa                                                                                                                    >  1:586106/1‑150 (MQ=255)                                tGATTAATCAGTAAATGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCaa                                                                                                                   >  1:32168/1‑150 (MQ=255)                                    tAATCAGTAAATGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAg                                                                                                               <  1:341592/150‑1 (MQ=255)                                         aGTAAATGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGc                                                                                                         >  1:135083/1‑151 (MQ=255)                                         aGTAAATGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGc                                                                                                         <  1:542798/151‑1 (MQ=255)                                            aaaTGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAAc                                                                                                      >  2:104981/1‑151 (MQ=255)                                              aTGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGt                                                                                                    >  2:406060/1‑151 (MQ=255)                                                       cAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCAc                                                                                           >  1:385040/1‑151 (MQ=255)                                                         aTTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCg                                                                                         <  2:582482/151‑1 (MQ=255)                                                               ccTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGc                                                                                    >  2:728058/1‑150 (MQ=255)                                                                cTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTg                                                                                                                         <  2:259752/112‑1 (MQ=255)                                                                cTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTg                                                                                                                         >  1:259763/1‑112 (MQ=255)                                                                      tcacttcaATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCAttttt                                                                              <  1:93475/149‑1 (MQ=255)                                                                        acttcaATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCTTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCg                                                                          <  2:211560/151‑1 (MQ=255)                                                                           tcaATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGtt                                                                       >  2:413685/1‑151 (MQ=255)                                                                                       gCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCAt                                                           >  2:678318/1‑151 (MQ=255)                                                                                         cAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCg                                                         <  2:61417/151‑1 (MQ=255)                                                                                           gTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGAc                                                       <  2:255924/151‑1 (MQ=255)                                                                                            ttttttCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGAcc                                                      >  2:753903/1‑151 (MQ=255)                                                                                               tttCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCAc                                                    <  2:135080/150‑1 (MQ=255)                                                                                                ttCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACtt                                                  >  1:526315/1‑151 (MQ=255)                                                                                                ttCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACt                                                   <  1:71462/150‑1 (MQ=255)                                                                                                   gCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCa                                               <  1:748093/151‑1 (MQ=255)                                                                                                    ccATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCAt                                              >  1:655965/1‑151 (MQ=255)                                                                                                      aTTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCAtct                                            >  2:384592/1‑151 (MQ=255)                                                                                                        tGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATctct                                          >  2:311590/1‑151 (MQ=255)                                                                                                         gCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATctct                                          <  2:546977/150‑1 (MQ=255)                                                                                                         gCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATctct                                          <  1:33965/150‑1 (MQ=255)                                                                                                            tCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTc                                       <  1:333709/150‑1 (MQ=255)                                                                                                              gCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAgg                                    <  2:425366/151‑1 (MQ=255)                                                                                                               cccTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGa                                   <  1:285316/151‑1 (MQ=255)                                                                                                                   ggCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAg                                 <  1:246322/149‑1 (MQ=255)                                                                                                                   ggCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAg                                 <  2:488064/149‑1 (MQ=255)                                                                                                                    gCTGCCATCTTCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAGATc                              <  2:480432/151‑1 (MQ=255)                                                                                                                              tCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAGATCGCACCtttt                     <  1:12316/150‑1 (MQ=255)                                                                                                                              tCGCGCAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAGATCGCACCtt                       <  2:172137/148‑1 (MQ=255)                                                                                                                                gcgcAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAGATCGCACCTTTTTCt                  >  1:583842/1‑151 (MQ=255)                                                                                                                                 cgcAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAGATCGCACCTTTTTCtt                 <  1:574900/151‑1 (MQ=255)                                                                                                                                 cgcAGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAGATCGCACCTTTTTCtt                 >  1:311024/1‑151 (MQ=255)                                                                                                                                    aGGACGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAGATCGCACCTTTTTCTTTAt              >  1:684587/1‑151 (MQ=255)                                                                                                                                       aCGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAGATCGCACCTTTTTCTTTATCaa           <  1:172533/151‑1 (MQ=255)                                                                                                                                       aCGTAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAGATCGCACCTTTTTCTTTATCaa           <  1:75251/151‑1 (MQ=255)                                                                                                                                          tAAGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAGATCGCACCTTTTTCTTTATCAAac         >  2:14255/1‑150 (MQ=255)                                                                                                                                           aaGCAGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAGATCGCACCTTTTTCTTTATCAAacac       <  2:385019/151‑1 (MQ=255)                                                                                                                                               aGGTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAGATCGCACCTTTTTCTTTATCAAACACGGGc   <  1:492973/151‑1 (MQ=255)                                                                                                                                                ggTTACTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAGATCGCACCTTTTTCTTTATCAAACACGGGCa  <  1:521486/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   ACGAAAAATAAAGGGCGCGAATGCGCCCTCGTGATTAATCAGTAAATGGAATGACAATTTCGCCTGGCTTCACTTCAATGCCTTTCGCCAGTTTTTTCGCCATTGCTTCGCCCTGGCTGCCATCTTCGCGCAGGACGTAAGCAGGTTGCTGGTTAAAGTAATTGCGTAATGCCTGGTTCAAATAGGGAAGCAACGTTTGCATCACCGTTTGCATTTTTTCCGGTTGTACCGTCGCATCGACCACTTCCATCTCTTTCAGGAAGATCGCACCTTTTTCTTTATCAAACACGGGCA  >  NC_000913/1122679‑1122972

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