Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,936,477	A→G	N67S (AAC→AGC)	rbsB →	ribose ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,936,477	0	A	G	100.0%	184.3 / NA	54	N67S (AAC→AGC)	rbsB	ribose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/24);  total (30/24)

AATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATC  >  NC_000913/3936341‑3936610                                                                                                                                         |                                                                                                                                      aaTGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGGGGTGCTGGACTCCCAGAACAGccc                                                                                                                                    >  2:195066/1‑140 (MQ=255)      gATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGa                                                                                                                               <  2:697470/140‑1 (MQ=255)       aTGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGaa                                                                                                                              <  1:348711/140‑1 (MQ=255)               aGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTgg                                                                                                                      <  2:912742/140‑1 (MQ=255)                  cacCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGa                                                                                                                   <  2:306210/140‑1 (MQ=255)                  cacCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGa                                                                                                                   <  1:198410/140‑1 (MQ=255)                  cacCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGa                                                                                                                   <  1:243687/140‑1 (MQ=255)                   acCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGaa                                                                                                                  >  2:902136/1‑140 (MQ=255)                       tCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTg                                                                                                              <  1:87727/140‑1 (MQ=255)                       tCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTg                                                                                                              <  1:5514/140‑1 (MQ=255)                         gcgcTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCa                                                                                                            >  1:320956/1‑140 (MQ=255)                             tggtggTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGAc                                                                                                        >  2:640478/1‑140 (MQ=255)                                    ctcCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCg                                                                                                 >  1:60620/1‑140 (MQ=255)                                          gCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTcgcg                                                                                           <  2:949617/140‑1 (MQ=255)                                                    ccGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAAt                                                                                 >  1:296317/1‑140 (MQ=255)                                                        tCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTctg                                                                             >  1:291343/1‑140 (MQ=255)                                                             gTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGAt                                                                        <  2:703993/140‑1 (MQ=255)                                                              tATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGAtt                                                                       >  2:110595/1‑140 (MQ=255)                                                              tATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGAtt                                                                       >  1:395832/1‑140 (MQ=255)                                                              tATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAAATTGGCTATAACCTGGCGGCGCACGGCTCCCCGAACAGCCCGGCGTAAGAGCTTGCGTACAGGCATCACTTTACCGATCACGGCACAAAAAATctgctgtgt                                                                       >  2:1014624/1‑137 (MQ=255)                                                                   cTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAAccc                                                                  >  1:583600/1‑140 (MQ=255)                                                                   cTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAAccc                                                                  >  2:104741/1‑140 (MQ=255)                                                                    tGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACccg                                                                 >  1:653474/1‑140 (MQ=255)                                                                      aaaGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACccgac                                                               >  2:978186/1‑140 (MQ=255)                                                                       aaGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACccgacc                                                              <  2:901428/140‑1 (MQ=255)                                                                              cgcgcAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGGTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCg                                                       >  2:966908/1‑140 (MQ=255)                                                                                cgcAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTTGACTCCCAGAACAGCCCGGCGAAAGAACTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGAc                                                     >  1:414351/1‑140 (MQ=255)                                                                                cgcAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGAc                                                     >  2:46442/1‑140 (MQ=255)                                                                                cgcAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGAc                                                     >  1:136952/1‑140 (MQ=255)                                                                                cgcAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGAc                                                     <  2:153484/140‑1 (MQ=255)                                                                                 gcAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACg                                                    <  1:902136/140‑1 (MQ=255)                                                                                   aGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCa                                                  <  1:640478/140‑1 (MQ=255)                                                                                   aGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCa                                                  <  2:653474/140‑1 (MQ=255)                                                                                   aGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCa                                                  <  1:46442/140‑1 (MQ=255)                                                                                   aGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCa                                                  <  2:1005283/140‑1 (MQ=255)                                                                                     aaaGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGt                                                <  1:104741/140‑1 (MQ=255)                                                                                          ggCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTa                                           <  1:394066/140‑1 (MQ=255)                                                                                          ggCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTa                                           <  2:320956/140‑1 (MQ=255)                                                                                           gCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTaa                                          >  1:113016/1‑140 (MQ=255)                                                                                               aTAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTggg                                             <  1:515804/133‑1 (MQ=255)                                                                                               aTAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTggg                                             >  2:515804/1‑133 (MQ=255)                                                                                                              aCCTGGTGGTGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAAc                       <  1:1016678/140‑1 (MQ=255)                                                                                                                  ggtggtGCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAgg                   >  2:1048308/1‑140 (MQ=255)                                                                                                                   gtggtgCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGc                  >  2:807269/1‑140 (MQ=255)                                                                                                                   gtggtgCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGc                  >  1:510083/1‑140 (MQ=255)                                                                                                                      gtgCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGaa               >  1:708567/1‑140 (MQ=255)                                                                                                                      gtgCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGaa               >  2:186573/1‑140 (MQ=255)                                                                                                                        gCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACa             >  1:595606/1‑140 (MQ=255)                                                                                                                        gCTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACa             >  1:297465/1‑140 (MQ=255)                                                                                                                         cTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACATCAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAAc                       <  2:969770/129‑1 (MQ=255)                                                                                                                         cTGGACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAAc                       >  1:969770/1‑129 (MQ=255)                                                                                                                            gACTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATccc         <  2:902535/140‑1 (MQ=255)                                                                                                                              cTCCCAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCgg       >  1:789298/1‑140 (MQ=255)                                                                                                                                 ccAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTa    >  2:439603/1‑140 (MQ=255)                                                                                                                                 ccAGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTa    >  2:872693/1‑140 (MQ=255)                                                                                                                                   aGAACAGCCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATc  >  1:387559/1‑140 (MQ=255)                                                                                                                                     aacagcCCGGCGAAAGAGCTGGTGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGa                <  2:621147/124‑1 (MQ=255)                                                                                                                                     aacagcCCGGCGAAAGAGCTGGTGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGa                >  1:621147/1‑124 (MQ=255)                                                                                                                                         |                                                                                                                                      AATGCGATGGCAAAAGACACCATCGCGCTGGTGGTCTCCACGCTTAACAACCCGTTCTTTGTATCGCTGAAAGATGGCGCGCAGAAAGAGGCGGATAAACTTGGCTATAACCTGGTGGTGCTGGACTCCCAGAACAACCCGGCGAAAGAGCTGGCGAACGTGCAGGACTTAACCGTTCGCGGCACAAAAATTCTGCTGATTAACCCGACCGACTCCGACGCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATC  >  NC_000913/3936341‑3936610

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